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Abstract 



The express purpose of these Lecture Notes is to go through some aspects of 



the simphcial quantum gravity model known as the Dynamical Triangulations 
approach. Emphasis has been on lying the foundations of the theory and on 



X 

H 

I illustrating its subtle and often unexplored connections with many distinct 

mathematical fields ranging from global riemannian geometry, moduli theory, 
number theory, and topology. Our exposition will concentrate on these points 
so that graduate students may find in these notes a useful exposition of some of 
the rigorous results one can establish in this field and hopefully a source of in- 
spiration for new exciting problems. We also illustrate the deep and beautiful 
interplay between the analytical aspects of dynamical triangulations and the 
results of MonteCarlo simulations. The techniques described here are rather 
novel and allow us to address successfully many high points of great current 
interest in the subject of simplicial quantum gravity while requiring very lit- 



tie in the way of fancy field theoretical arguments. As a consequence, these 
notes contains mostly original and yet unpublished material of great potential 
interest both to the expert practitioner and to graduate students entering in 
the field. Among the topic addressed here in considerable details there are: 
(i) An analytical discussion of the geometry of dynamical triangulations in 
dimension n = 3 and n = 4; (ii) A constructive characterization of entropy 
estimates for dynamical triangulations in dimension n = 3, n = 4, and a com- 
parision of the analytical results we obtain with tha data coming from Monte 
Carlo simulations for the 3-sphere and the 4-sphere (Hi) A proof that in 
the four-dimensional case the analytical data and the numerical data provide 
the same critical line ^4(^2) characterizing the infinite- volume limit of simpli- 
cial quantum gravity; (iv) An analytical characterization of the critical point 
/c^*"** which is in an extremely good agreement with the location of the critical 
point obtained by MonteCarlo simulations; (v) A proof that corresponding to 
such a critical point, 4-dimensional (simply connected) manifolds undergo a 
higher order phase transition possibly concealed by a bistable behavior: this 
settles down in the affirmative a long standing debate as to the nature of the 
phase transition in 4-dimensional simplicial quantum gravity; (vi) A geomet- 
rical analysis of the branched polymer phase of 4D-gravity; (vii) We also show 
that in the three- dimensional case the confrontation between the analytical 
and numerical data is very satisfactory; (viii) A proof that in dimension n = 3 
there are pronounced hysteresis effects which are possibly responsible for the 
observed first-order nature of the phase transition at criticality. 
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I. INTRODUCTION 



Recent years witnessed a massive introduction of methods of statistical field theory in 
attempts to quantize gravity along the lines of a conventional field theory fill? 0- Such an 
emphasis is the direct outspring of the possibility of discretizing gravity in a way consistent 
with the underlying reparametrization invariance of general relativity [Q. The basic 
idea is to use a variant of the standard Regge calculus [^], first suggested by Weingarten 

in an attempt to make quantum gravity well defined as a lattice statistical field theory. 
This discretization is known as the theory of Dynamical Triangulations (DT) or simplicial 
quantum gravity. It consists of a replacement of the (Euclidean) functional integral over 
all equivalence classes of metrics with a summation over all abstract triangulations of the 
given manifold. The fixed edge-length of the links of the triangulations plays here the 
role of lattice spacing of the more familiar lattice gauge theories. The basic idea is to 
search, in the parameter space of these discretized quantum gravity models, for critical 
points where the lattice spacing can be taken to zero and contact can be made to continuum 
physics. Credit to this picture is lend by noticing that within this framework, Euclidean 
two-dimensional quantum gravity can be formulated as the scaling limit of a dynamically 
triangulated collection of surfaces thought of as an ordinary statistical system. Contact 
with the continuum, as described by Liouville theory can successfully be made, and the 
robust alliance of the theory of critical phenomena and Monte Carlo simulations even sheds 
light on many aspects of the theory not liable of an analytical approach. 

If we extend the theory to higher dimensions we have to face the fact that we have no 
continuum theory of Euclidean Quantum Gravity with which we can compare, nonetheless 
many aspects of the theory can be investigated by computer-assisted simulations with a good 
degree of precision. These numerical methods allow us to get a rather complete picture of 
the phase diagram of the discretized theory and they serve as an important inspiration for 
analytical studies of higher dimensional models of quantum gravity. These lecture notes 
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attempts to partially fill the gap between such numerical studies and a fully- fiedged ana- 
lytical approach to higher dimensional dynamical triangulations. Our emphasis here is on 
the applications of elementary ideas of Piecewise-Linear (PL) geometry, global riemannian 
geometry, and number theory in order to understand dynamically triangulated models of 
quantum gravity in dimension n = 3,4. The occasion for such an analytical approach comes 
by exploiting the techniques developed for the recent proof of the conjecture concerning 
the existence of an exponential bound to the number of distinct triangulations on a PL- 
manifold M", {n > 3), of given volume and topology. This exponential bound is germane 
to the definition of dynamical triangulations as a lattice statistical field theory model. An- 
alytical proofs are well known in two dimensions 0, and computer simulations supported 
very strongly the existence of such bound also in higher dimensions [|l|. The actual proof 
for higher dimensions draws from techniques of controlled topology and geometry initiated 



by Cheeger, Gromov, Grove, and Petersen 0, |10|, and it suggests that such mathematical 
methods may bear relevance to the whole program of dynamical triangulations. A check 
with Monte Carlo simulations (in dimension n = 4) shows an excellent agreement between 
the numerical data and the analytical results obtained by pursuing this geometrical approach 
to higher dimensional dynamically triangulated models. It is not yet possible to reach the 
level of sophistication of the two-dimensional case, but as we show in these notes, the an- 
aljd;ical methods we develop in dimensions n = 3 and n = 4, provide a good start for an 
analytical understanding of dynamical triangulations in such dimensions. In particular, we 
characterize with greater precision the entropy estimates 0, and provide the corresponding 
generating functions. Among other properties common to both the 3- dimensional and 4- 
dimensional case, we prove that in 4- dimensions dynamically triangulated models it makes 
sense to sum over all simply connected manifolds. We also discuss the existence of a polymer 
phase in these models [0. We analytically prove that the system admits a critical point 
^2 = kf^^ corrisponding to which a dynamically triangulated 4-manifold undergoes a higher 
order phase transition possibly concealed by a bystable behavior. This is an important re- 
sult suggesting that dynamical triangulations provide a reliable candidate to a consistent 



quantum gravity theory, and shows that the importance of the interplay between analytical 
approaches and Monte Carlo simulations continues unabated. 

The arguments presented here should be amenable to considerable improvements, and 
hopefully these notes can be used as a working tool. For this reason, we have presented 
a detailed summary of the more important and useful results in section III. In section II 
we review the definition of the model. In section IV we recall few elementary aspects of 
the PL-geometry of dynamical triangulations. In section V we establish the connection 
with Gromov's spaces of bounded geometry. This connection provides mathematical foun- 
dation to the heuristic arguments which motivate the approximation of distinct riemannian 
structures with distinct dynamical triangulations. In section VI we introduce the Moduli 
space associated with dynamical triangulations. These moduli spaces parametrize the set 
of inequivalent deformations of dynamical triangulations needed to approximate riemannian 
structures with large symmetries. In sections VII and VIII we provide all entropy esti- 
mates needed for the higher dimensional DT models. Here we exploit elementary number 
theory in order to provide such estimates and construct the associated generating functions. 
We discuss also the 2-dimensional case which fully agrees with the known results. The 4- 
dimensional case, and the 3-dimensional case are discussed in section IX. The agreement 
between the analytical results we obtain and the Monte Carlo simulations is impressive in 
dimension n = 4. We discuss in detail an analytical proof of the existence of the phase 
transition whose presence is strongly suggested by the computer simulations. As recalled we 
prove that definitively it is not a first order transition as hinted by some recent simulations. 

Finally, we show that also in dimension n = 3 the analysis of the data coming from the 
two approaches is quite satisfactory. In particular, we prove that the correct asymptotics 
for the canonical partition function has a bistable behavior yielding for the hysteresis effects 
which are believed to signal a first order nature of the transition point. 
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II. THE MODEL: SIMPLICIAL QUANTUM GRAVITY 



Let M be an n-dimensional, (n > 2), manifold of given topology, and with a finite num- 
ber of fixed (n — l)-dimensional boundaries S^, k = 1,2,.... Let Riem{M) and Diff{M) 
respectively denote the space of riemannian metrics g on M, and the group of diffeomor- 
phisms on M. In the continuum formulation of quantum gravity the task is to perform a 
path integral over equivalence classes of metrics: 

Z(A,G,S,,/.)= Y: I P[^?(M)]e-^«[^>^'^l (1) 

Top{M) JRiemiM)/Diff{M) 

weighted with the Einstein- Hilbert action associated with the riemannian manifold {M,g), 
viz., 

Sg[A,G,^ = Af (PiJ7j-^^=l (riJ7jR+ boundary terms (2) 

The boundary terms depend on the metric h and on the extrinsic curvature induced on the 
boundaries S^, and they are such that for the action obtained by gluing any two manifolds 
(Ml, and (M2, g2) along a common boundary (S, h), we get Sg^^g,^ [A, G] = Sg-^ [A, G, S] + 
5*32 [A, G, S]. In this way, the partition function (P associated with a manifold M with two 
boundaries and {T,2,h2), satisfies the basic composition law 

hi] S2, = Yu ^(^1' hi)Z{i:i, hi] S2, h2) (3) 

which describes how M can interpolate between its fixed boundaries {lli,hi) and (S2,/i2) 
by summing over all possible intermediate states (Sj, hi). 

Even without pretense of rigor, it is clear that the formal a priori measure P[5f(M)] char- 
acterizing such sort of path integration should satisfy some basic properties. In particular: 
(i) 'D[g{My\ should be defined on the space of riemannian structures Riem{M)/Diff{M) 
so to avoid counting as distinct any two riemannian metrics gi{M) and g2{M) which differ 
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one from the other simply by the action of a diffeomorphism of 0: M — > M, viz., such that 
g2{M) = <P*g^{M)- 

(a) The measure Vlgi^M)] should play a kinematical role and not a dynamical one. More 
explicitly, for a given r G M"*", let [Riem{M) / Dif f{M)]x^^x2;r denote the set of all riemannian 
structures on a manifold M with two marked points Xi and X2 with preassigned distance 
dg(M){xi,X2) = r. When restricted to this set of riemannian structures, Vlgi^M)] should 
factorize: fluctuations in the geometry of M which are localized in widely separated regions, 
should be statistically independent. This implies that 'D[g{My\ should not describe the a 
priori existence of long range correlations on the set of riemannian manifolds considered. 
Such correlations should be generated by the spectrum of fluctuations of (|l]) by means 
of a statistical suppression-enhancement mechanism similar to the energy versus entropy 
argument familiar in statistical mechanics. In a field theoretic formalism, this property is 
translated in the familiar requirement of ultra-locality, namely in the absence of derivatives 
in the formal LF' norm on the space of metrics, (the De Witt supermetric) . 
(Hi) Finally, V[g{M)] should be so constructed as to allow for the introduction of a Diff{M)- 
invariant short distance cut-off representing the shortest wavelength allowed in discussing 
fluctuations of the geometry of M. This cut-off should be removable under appropriate 
circumstances, in particular when, in a critical phase, long range correlations are generated. 

Summing over topologies. Long before addressing the physical and mathematical charac- 
terizations of such requirements on V[g{M)], both in the continuum field-theoretic formalism 
and in its possible discretized versions, (in this connection see the recent work of Menotti 
and Peirano one must also discuss the proper meaning to attribute to the formal sum 
over topologies appearing in (|l]): are we summing over homotopy types, homeomorphism 
types, or over smooth types of manifolds? 

The homotopy type of a manifold is a topological notion which (with the exception of 
dimension two) is too weak to appear of some immediate utility. The homeomorphism 
type of a manifold M is the notion which is more directly related to a natural summation 
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over topologies, however it is difficult to handle in dimension larger than two, (see below). 
Finally, the smooth type yields for a summation over the possible distinct differentiable 
structures that a manifold M can carry. This latter summation appears more relevant to a 
field- theoretic formalism since in an expression such as (|I]) the diffeomorphism group of the 
underlying manifold Diff{M) plays a basic role, but as we shall see momentarily, a sum- 
mation over the smooth types arises more naturally in simplicial gravity. Thus, this latter 
interpretation for J2top{m) is perhaps the most appeahng since it suggests an unexpected 
bearing of simplicial quantum gravity on the field-theoretic formalism, and in this sense 
the two approaches are not alternative to each other. Simplicial quantum gravity may be 
necessary in addressing the important issue of summation over distinct topologies, whereas 
the field-theoretic formalism is necessary in order to understand the nature of the continuum 
limit (if any exists) of the discretized models. 

In dimension two there is equivalence between smooth structures and the homeomor- 
phism types, and the summation over topologies can be given the unambigous meaning 
of a summation over the Euler characteristic x{M) of the manifolds, since this invariant 
characterizes surfaces. However, even in such a simple case it is very difficult to obtain 
clear-cut results for what concerns a reliable procedure for summing over topologies. The 
use of matrix models of 2Z}-gravity and the associated double scaling limit have shed some 
light on this issue, but we are still far from a clear understanding [|I|. 



In dimension three topological manifolds are uniquely smoothable |T2|, and the ques- 
tion of summing over topologies is again synonimous of summation over smooth structures. 
Thus, as long as we confine our attention to the topological category or to the smooth cate- 
gory, summation over topologies is reduced to the yet unsolved problem of enumerating the 
homeomorphism types of three-manifolds. This enumeration cannot be realized as long as 
the Poincare conjecture is not proved. For instance, if there were a fake three- sphere then 
one could prove |[T^ that there cannot be finitely many homeomorphism types of three- 
manifolds even under bounds on curvatures of the manifold, (typically, under curvature, 
volume and diameter bounds one gets topological finiteness theorems for homeomorphism 
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and diffeomorphism types ||T^ for dimension n 7^ 3,4). 



In dimension four the situation is even more complex since smoothing theory is not yet 
completely known. In open contrast to the theory of 3-manifolds, there is no equivalence be- 
tween topological manifolds and smooth structures, and Donaldson-Freedman's theory shows 
that there are topological manifolds which are not smoothable as well as manifolds admitting 
uncountably many inequivalent smooth structures |T^. The situation is further worsened 
by the fact that the topological classification of all (compact, orientable) 4-manifolds is log- 
ically impossible, (again, the fundamental group should be blamed for this). As stressed by 
Frohlich this circumstance has even be used to foster the credence that only simply- 



connected, spinable 4-manifolds should contribute to the gravitational path integral. But 
clearly this point of view, if not better substained, cannot be advocated. Suppression of a 
class of manifolds from path integration can be justified only on a dynamical ground, as a 
form of statistical suppression driven by the spectrum of fluctuations of the theory. 

This situation does not improve even if we confine our attention to a manifold M of fixed 
topology and with a given smooth structure, for, in that case we have the additional problems 
associated with: (i) The unboundedness of the Euclidean action; (ii) The mathematical 
difficulties in defining a proper path integration over the stratified manifold of riemannian 
structures; (in) The Einstein-Hilbert action is not renormalizable. 

Simplicial quantum gravity. The hope behind simplicial quantum gravity is that some 
of the above problems, concerning both the characterization of the measure T>[g{My\ or the 
issues related to T^TopiM)^ can be properly addressed, in a non-perturbative setting, by ap- 
proximating the path integration over inequivalent riemannian structures with a summation 
over combinatorially equivalent piecewise linear manifolds. The first attempt of using PL 
geometry in relativity dates back to the pioneering work of Regge . His proposal was to 
approximate Riemannian (Lorentzian) structures by PL-manifolds in such a way as to obtain 
a coordinate-free formulation of general relativity. The basic observation in this approach is 
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that parallel transport and the (integrated) scalar curvature have natural counterparts on 
PL manifolds once one gives consistently the lengths of the links of the triangulation defining 
the PL structure. The link length is the dynamical variable in Regge calculus, and classically 
the PL version of the Einstein field equations is obtained by fixing a suitable triangulation 
and by varying the length of the links so as to find the extremum of the Regge action. If 
the original triangulation is sufficiently fine, this procedure consistently provides a good ap- 
proximation to the smooth spacetime manifold which is the corresponding smooth solution 
of the Einstein equations. This approach can also be successfully extended so as to provide 
a quantum analogue to Regge calculus 0, in which we replace the formal path integration 
over the space of Riemannian structures with an integration over the link variables. 

Dynamical triangulations are a variant of Regge calculus in the sense that in this for- 
mulation the summation over the length of the links is replaced by a direct summation over 
abstract triangulations where the length of the links is fixed to a given value a. In this 
way the elementary simplices of the triangulation provide a Diff-invariant cut-off and each 
triangulation is a representative of a whole equivalence class of metrics. Regge calculus 
still mantain its validity and it provides both the metric assignement for the PL manifold, 
obtained by gluing the simplices, and the corresponding action. Since all simplices now are 
identical, the action will only depend on the numbers, respectively and Nn^2i of n- and 
(n — 2)-dimensional simplices of the n-dimensional PL manifold. In this way we get that 
the Einstein- Hilbert action for n-dimensional (Euclidean) gravity formally goes into the 
combinatorial action 

where |A^„_2(<9T) is the boundary term, and A;„, and /c„_2 are (bare) coupling constants 
related to the cosmological constant A and to the gravitational coupling G, respectively. In 
particular, we can view l/kn-2 as a bare gravitational coupling constant. 
The partition function associated with such discretized action is 
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TeT 

where the summation is over distinct triangulations T in a suitable class of triangulations T. 
Roughly speaking, we consider any two triangulations (say with the some number of vertices) 
distinct if there is no map between the vertices which is compatible with the assignements of 
links, triangles, etc., while T restricts the class of triangulations to those satisfying suitable 
topological constraints. The proper choice of the class T is strictly connected to the difficult 
question of how to sum over topologies in quantum gravity mentioned above. 

The interpretation of (^) as providing also a sum over topologies is based on the observa- 
tion that a PL manifold is uniquely smoothable in low dimension (in particular for n = 4, by 
Cerf theorem, see e.g., fl^)) (it is also worth stressing that every smooth manifold admits 
a natural C^-smooth PL-structure). In particular, in dimension n = A there is a bijective 
correspondence between (isotopy) classes of smooth structures and PL-structures. Thus (U), 
if no restrictions are imposed on T, subsumes in a rather natural way a summation over 
topologies if we decide that J^Top actually means summation over smooth structures. Even 
if not ambiguous in its definition, (^ blends summation over metric structures and summa- 
tion over smooth structures only in a formal way. Already in dimension two, the sum (^ 
is divergent since if the topology is not fixed, the number of distinct triangulations grows 
factorially with the volume of the manifold {i.e., with the number of simplices). However, if 
one fixes the topology, the number of distinct triangulations of a 2-dimensional PL manifold 
is exponentially bounded |^ according to 

(A)^^iV2('^)(^--2)-^ (6) 

where A is a suitable constant and 'jstr, the string exponent, is a topological suscettivity 
generated by the quantum fiuctuations of the metric. In this way, Y^ti^t where T is 

a given PL surface, is well defined. Topology is then allowed to fiuctuate {i.e., one attempts 
to extend the summation to all PL surfaces), by a delicate limiting procedure {Double 
Scaling Limit). In higher dimension is not yet known how to perform the summation (^, 
but numerical simulation as well as the results of the two-dimensional theory, suggests 

13 



that makes sense if one fixes the topology. This issue is clearly related to a rigorous 
characterization of the discretized counterpart of the formal a priori measure T>[g{M)]. As a 
matter of fact, a necessary condition for attributing a meaning to (|^), for fixed topology, is to 
require that the number of distinct triangulations of a given PL manifolds is exponentially 
bounded as a function of the number of (top-dimensional) simplices Recently |]^, the 
existence of the exponential bound in all dimensions and for all (fixed) topologies has been 
proved. The existence of the exponential bound implies that for a fixed kn-2 there is a 
critical lower value /c^'"**(A;„_2) of kn such that the partition function is well defined for 
kn > k^^*{kn-2) and divergent for fc„ < k'^''^ {kn-2) ■ To be more specific this implies that 
we can introduce the (canonical) partition function for fixed (lattice) volume: The effective 
Entropy 

ly (AT, fc„_2)e// = E e^"-^"-(^) (7) 
r6r(Af„) 

where T{Nn) denotes the class of distinct triangulations of fixed volume (A^n), fixed topology 
and boundary conditions. This effective entropy will characterize the infinite (lattice) volume 
limit of the theory, defined by the approach to the critical line fc„ — > kn^^{kn-2) in the 
{kn-2, kn) coupling constant plane. The existence of the infinite volume limit is a necessary 
but not sufficient condition for the existence of a physically significant continuum limit of 
the theory, and we only expect interesting critical behavior, i.e., the onset of long range 
correlations, at certain critical values of kn-2- In the rest of these notes we will be mainly 
interested in the analytical characterization of the canonical partition function (the effective 
entropy) W{N, kn-2)eff, and in comparing the obtained results with the existing numerical 
data coming from Monte Carlo simulations. 

It is not difficult to see the origin of the difficulties in dealing with the set of dynamically 
triangulated manifolds (of given volume and topology) considered as a statistical system 
described by (0). Roughly speaking we have a collection of identical simplices whose only 
interaction is basically associated with unpenetrability and gluings according to certain rules. 
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The energetic term is very simple and controls the average curvature. Thus the free energy 
of the system is basically characterized by the entropic factor enumerating the number of 
distinct dynamical triangulations of given topology, volume, and average curvature. Such 
an estimate is deeply non-local, {e.g., see (^) for n = 2), and as such its characterization 
is conceptually different from the standard approaches used to enumerate distinct configu- 
rations, say of spins, on a rigid lattice. This situation is not new in statistical mechanics 
since it is reminiscent of what happens when dealing with the hard sphere gas: a set of 
identical spherical particles whose only interaction is associated with unpenetrability. Here 
too, the free energy is of entropic origin, and the characterization of such an entropy is a 
highly non local problem since it is related to the characterization of the densest packings 
in spheres: the insertion of a sphere may change the packing up to a distance proportional 
to the inverse of the average separation between the surfaces of the spheres. For dynamical 
triangulations, as we shall see, the source of non- locality is subtler since there is a delicate 
feedback with the topology, {e.g., see (|]), where the critical exponent 7^ depends from the 
Euler characteristic of the surface obtained by gluing the simplices). But these topologi- 
cal difficulties would be present also in the hard sphere gas had we considered the spheres 
evolving in a topological non-trivial ambient space. Thus, it is perhaps not so surprising 
to note that the first proof of the existence of an exponential bound for the entropy of dy- 
namical triangulations was carried out by exploiting the geometry of sphere packings in 
riemannian manifolds. This interplay between the geometry of packings, topology, and the 
metric properties of riemannian manifolds, (or more general metric spaces), can be traced 
back to M. Gromov's ideas. These ideas have sprung a renaissance in recent developments 
in Riemannian geometry, and we will exploit them here by showing that they are a good 
source of inspiration also in simplicial quantum gravity and in the fascinating field of the 
statistical mechanics of extended objects. 
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III. SUMMARY OF RESULTS 



For the convenience of the reader we present here a summary of the constructive re- 
sults that we prove in these notes and that can be useful in applications to higher dimen- 
sional dynamical triangulations. If not otherwise stated, we refer to dynamically triangu- 
lated manifolds, M, in dimension n > 2. The dynamical triangulations in question have 
Nn-2{T) = A + 1 bones cr""^, Nn(T)=N top- dimensional simplices u", a given average num- 
ber, b{n,n — 2) = ^n{n + 1){N/X + 1), of n-simplices incident on a bone, and a minimum 
number q of n-dimensional simplices incident on a bone, (typically q — 3). Sometimes we 
need to mark n of these bones and a top-dimensional simplex ctq incident on one of them. 
We refer collectively to triangulations with such markings as rooted triangulations. The 
topology of the triangulations explicitly enters our results through the representation variety 
Hom{TvijM),G) . ^j^g q£ conjugacy classes of representations of the fundamental group of the 
manifold, 7ri(M) into a (compact) Lie group G. This representation variety parametrizes 
the set of inequivalent deformations of dynamical triangulations needed to approximate rie- 
mannian structures endowed with G-structures, (for instance, this parameter space tells us 
how to deform a 3-dimensional dynamical triangulation so as to approximate unambiguosly 
inequivalent flat 3- tori). Finally, in what follows we also introduce a normalizing factor 
< c„ < 1, of the form 

c„ = CoA^;(") • e""^" (8) 

where the constants < Co < 1, r(n) > 0, an < 0, depend only on the dimension n. 
These normalizing constants are necessary since our enumeration procedure (in establishing 
entropy estimates) overcounts the number of distinct dynamical triangulations admitted by 
a manifold of given topology. In dimension n = 2, c„ reduces to a constant, whereas in 
dimension n = 4, r(n = 4) = 5/2 and Q;n=4 = — | ln(cos~^(l/4)) and similarly in dimension 
n — 3. With these notational remarks along the way, we have: 
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The generating function. Let < t < 1 be a generic indeterminate, and let p\{h) 
denote the number of partitions of the generic integer h into (at most) X + n — 1 parts, each 
< {b{n,n — 2) — q){X + 1). In a given holonomy represention G:7ri(M; ctq) — >■ G, and for 
a given value of the parameter b — b{n,n — 2), the generating function for the number of 
distinct rooted dynamical triangulations with Nn-2{Tj^'') = A + 1 bones and given number 
((X h) of n- dimensional simplices incident on the n marked bones is given by 



{b{n, n - 2) - g)(A + 1) + (A + 1) - n 
(A + 1) - n 

where < c„ < 1 is the above-mentioned normalizing factor, and 



g[W{e,X,b;t)]^Cr,-J2Pxih)t''^c„ 

h>0 



(9) 



n 
m 



;i - r)(i - 1 



...il-f 



m+1 ^ 



(1 _ tn-m)^l - tn-m-l) . . . (1 _ t) 

is the Gauss polynomial in the variable t. 



(10) 



The Rooted Entropy. In a given represention ©:7ri(M; a^ ) G, and for a given value 
of the parameter b = b{n, n — 2), the number W{Q, A, b) of distinct rooted dynamical trian- 
gulations with Nn_2{T^^^) = A + 1 bones, is given by 



w{e,x,b)^cn 



^ {b{n, n - 2) - g)(A + 1) + (A + 1) - n ^ 



V 



(A + 1) -n 



(11) 



Asymptotics for the Entropy. The number of distinct dynamical triangulations, with 
A + 1 bones, and with an average number, b = b{n, n — 2), of n- simplices incident on a bone, 
on an n- dimensional, (n < i), PL- manifold M of given fundamental group 7ri(M), can be 
asymptotically estimated according to 



W{X,b) ~ 



ip-q 



\l-2n 



^ ^/2^\ {b-qf 



{b-q + 1)^-9+1 



n A+i 



b-q 



b{n,n — 2) ^.^/o.^ 2n+3 



n{n + 1) 



(12) 
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where W^^ is a topology dependent parameter, and D^dim[Hom{Tri{M),G)]. 



The Canonical Pcirtition Function. Let us consider the set of all simply-connected n- 
dimensional, (n = 3, dynamically triangulated manifolds. Let us set 



27 / 27 

gfcn-2 + 1 + W [ — ekn-2 + 1)2 _ 1 

^ V ^ 



1/3 



+ 



27 



,27 



«„_2 + 1 _ w e'^"-^ + 12-1 



1/3 



- 1, (13) 



and let k^_2 denote the unique solution of the equation 

1 1 

r]*{kn-2) = -(1 - -TTi t) = Vn 



(14) 



3' A{K^2)' 

where rjmax = 1/4, (for n = 4) and rjmax = 2/9, (for n = 3). Let < e < 1 small enough, 
then for all values of the inverse gravitational coupling kn-2 such that 



K-2 - e < K-2 < +00, 



(15) 



the large N -behavior of the canonical partition function (effective entropy) is given by the 
uniform asymptotics 



27r 



^-n-l/2g(iV+l)/(r,„,,) 



WQ{lpmax{kn-2)\f^) + ^W_i{lpmax{kn-2)\l N) 



N 



(16) 



where Wr{z) = r(l — r)e^^^'^Dr-i{z) , (r < 1), Dr-i{z) and T{1 — r) respectively denote the 
parabolic cylinder functions and the Gamma function, and where the constants and are 
given by 

V^l/xi'^ - 2r/max)"" i^max{kn-2) 



6 = 



^(1 - 3r]max){i - '^Vmax) fviVmax) 



ri*-y^ (1 -2?7*)-" 



ipmax{kn-2) ^(1 - 3ri*){l - 2r]*) 



'^max{kn-2) \J ~fm {V* 



(17) 



(18) 
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Where, for notational convenience we have set 



f{r]) = -77log?7 + (1 - 2r]) log(l - 2r)) - (1 - Stj) log(l - Stj) + 



(19) 



with fjj = df/dr], and /^^=d^//d77^,and 

i^max{kn-2) = Sgn{r]* - r]max)\/'2[f(V*) " fiVmax)]- 



(20) 



(Note that for the 3-dimensional case the above expression for W{N, kn-2) is shghtly more 
general if one uses the variables (A^o, -/V3)-see below for details). 



The infinite volume limit. The critical value of the coupling k^'^^ corresponding to the 
infinite-volume limit for simply- connected n- dimensional dynamically triangulated manifolds 
is given by 



■ A{k2) + 2 In On 
m h lim 

3 N„^oo N„ 



(21) 



which, for < e < 1 small enough, holds for all /c^"! — e < kn-2 < k^l\. Whereas, for 

kn-2 < kn-2 < +00 wc get 



kT\kn-2) 



-n{n + 1) I log 



;i - 2r/, 



\{l-2rimax) 



r7te^(l-3w)^'"^'''"''^^ 



+ kn-2Vmax + ,lim 



In On 



(22) 



The four-dimensional case. The critical value of the coupling corresponding to the 
infinite-volume limit for a simply- connected A-manifold, (e.g., the 4-sphere S"^), is given by 



kr\k2)^lQ{log 



2r]„ 



\(l-2?j„ 



r/fc^(l-3r/™a,,)(i-3''— ) 



log 



2.066 



(23) 
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(with /cg"* < k2 < +00), and rjmax = 1/4, and 

kf'^ ~ 1.387 (24) 

(note that this is a theoretical upper bound; for instance for rjmax = 1/4.1 one would get 
j^cnt ^ 1.2), whereas one obtains 

^-* = 101n^lM±| (25) 
^ 3e/2.066 ^ ^ 

when kf^* — e < k2 < k^'^* . Corresponding to the critical value k2 — > kf^^ the system 
undergoes a higher order phase transition. 



The three-dimensional case. Let us consider the set of all simply- connected 3- 
dimensional dynamically triangulated manifolds described in terms of the variables N^, Nq. 
Let k^Q^ , and k^"^ respectively denote the unique solutions of the equations 



3^ A{ko) 



1 1 , 

3^^~A{koj' 



1 _ 1 

~l~ Vmin 

6 6 



1 

+ 6 



9 



(26) 



(27) 



Let < e < 1 small enough, then for all values of the inverse gravitational coupling ko such 
that 



k^^ + €<ko< A;^"* - e, 



(28) 



the large N -behavior of the canonical partition function W{N, ko)eff for 3-dimensional sim- 
plicial quantum gravity on a simply connected manifold is given by the uniform asymptotics 



^jV-7/2g(iV+l)/(»7^i„) 



Wo(V'mm(^o)V N) + -j^W^i{'^rnvn{ko)\J N) 



N 



+ 



+ 



27r 



^-7/2g(JV+l)/(r;^a,) 



Woi'tpmaxiko)]/^) + ^«^-l(V'max(^o) V-^) 



N 



(29) 
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where the constants ao, a\, and ^q, are given by 



\J{1- 3r7mi„)(l - 2rjrain) fviVmin) 



(30) 



^^-1/2 (1 _2ry* 



i^minih) .^(1 - 3r/*)(l - 2r7*) 



(31) 



^0 = 



max / Ymax 
(1 - 3?7„aa;)(l - 2r7„aa;) fviVmax) 



(32) 



6 = 



Co 



(1 -2r7*) 



V'max(fco) ^(1 - 3r/*)(l - 2r7*) 



(33) 



r/ie critical value of the coupling k^"^ corresponding to the infinite-volume limit for the 
3-sphere §^ is given by 



(34) 



and corresponding to the critical point 



ko k^"' ~ 3.845 



(35) 



the system described by ^2^) exhibits a bistable behavior leading to hysteresis effects. 
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IV. TRIANGULATIONS 



In this section we recall few basic notions of Piecewise-Linear (PL) topology in order 
to provide a self-contained set up for the some of the mathematical problems we wish to 
discuss in simplicial quantum gravity. Together with standard material, (for which we refer 
freely to [T^), we also collect few definitions and results which are not readily accessible in 
textbooks. 

A. Preliminaries: simplicial manifolds and pseudo-manifolds 

We first list some well known, but necessary, preliminaries. 

By an n-simplex cr" = (xq, . . . , x„) with vertices Xq, . . . ,Xn we mean the following sub- 
space of R'^, (with d > n), 

n 

cr" = ^A,x, (36) 

i=0 

where xq, . . . , are n + 1 points in general position in i?*^, and YJi=o Aj = 1 with Aj > 0. 

A face of an n-simplex a" is any simplex whose vertices are a subset of those of a", and 
a simplicial complex i^' is a finite collection of simplices in such that if a" G -fC then so 
are all of its faces, and if a", G K then cr" fl a™ is either a face of a" or is empty. The 
/i-skeleton of K is the subcomplex C K consisting of all simplices of K of dimension 
< h. 

Let i^' be a (finite) simplicial complex. Consider the set theoretic union \K\ G of all 
simplices from K 

\K\ = U^^K(T (37) 

Introduce on the set \K\ a topology that is the strongest of all topologies in which the 
embedding of each simplex into \K\ is continuos, (the set A C l-f^l is closed iS A H a'' is 
closed in for any cr^ G K). The space |K| is the underlying polyhedron, geometric carrier 
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of the simplicial complex it provides the topological space underlying the simplicial 
complex. The topology of 1/^1 can be more conveniently described in terms of the star of 
a simplex cr, star (a), the union of all simplices of which cr is a face. The open subset of 
the underlying polyhedron \K\ provided by the interior of the carrier of the star of a is the 
open star of a. Notice that the open star is a subset of the polyhedron \K\, while the star 
is a sub-collection of simplices in the simplicial complex K. It is immediate to verify that 
the open stars can be used to define the topology of \K\. The polyhedron \K\ is said to be 
triangulated by the simplicial complex K. More generally, a triangulation of a topological 
space M is a simplicial complex K together with a homeomorphism |i^| — >^ M. 

Simplicial maps. A simphcial map f: K ^ L between two simplicial complexes K and L 
is a continuous map f'\K\ — > \L\ between the corresponding underlying polyhedrons which 
takes n-simplices to n- simplices for all n, (piecewise straight- line segments are mapped to 
piecewise straight-line segments). The map / is a simplicial isomorphism if f~^: L ^ K is 
also a simplicial map. Such maps preserve the natural combinatorial structure of M". Note 
that a simphcial map is determined by its values on vertices. In other words, ii f: ^ 
carries the vertices of each simplex of K into some simplex of L, then / is the restriction of 
a unique simplicial map. 

Combinatorial properties. Sometimes, one refers to a simphcial complex X as a simphcial 

division of A subdivision K' of is a simplicial complex such that \K'\ = \K\ and each 
n-simplex of K' is contained in an n- simplex of K, for every n. A property of simplicial 
complex K which is invariant under subdivision is a combinatorial property or Piecewise- 
Linear (PL) property of K, and a Piecewise-Linear homeomorphism f:K^L between two 
simplicial complexes is a map which is a simplicial isomorphism for some subdivisions K' 
and L' of K and L. 

Piecewise-Linear manifolds. A PL manifold of dimension n is a polyhedron M — \K\ 
each point of which has a neighborhood, in M, PL homeomorphic to an open set in i?". 
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PL manifolds are realized by simplicial manifolds under the equivalence relation generated 
by PL homeomorphism. Any piecewise linear manifold can be triangulated, however notice 
that any particular triangulation of a piecewise linear manifold is not well defined up to 
piecewise linear homeomorphism, (because such homeomorphisms involve subdivisions of 
the given triangulations). A triangulated space can be characterized as a PL- manifold 
according to the |[16|| 

Theorem 1 A simplicial complex K is a simplicial manifold of dimension n if for all r- 
simplices G K, the link of , link{a^) has the topology of the boundary of the standard 
[n — r)-simplex, viz. if link {a^) ~ §>"'~^'^'^. 

Recall that the link of a simplex o" in a simplicial complex K is the union of all faces ct/ of 
all simplices in star{a) satisfying o"/ fl a = 0, also recall that the Cone on the link link{a^), 
C{link{a^)), is the product link{a^) x [0, 1] with link{a^) x {1} identified to a point. The 
above theorem follows [|T^ by noticing that a point in the interior of an r-simplex o"^ has a 
neighborhood homeomorphic to B"^ x C{link{a^)), where denotes the ball in R". Since 
link{a'') ^ and C{^^'^'''~^) ~ B"-^^, we get that \K\ is covered by neighborhoods 

homeomorphic to B^ x B^^^ ~ B^ and thus it is a manifold. Note that the theorem holds 
whenever the links of vertices are (n — l)-spheres. As long as the dimension n < 4, the 
converse of this theorem is also true. But this is not the case in larger dimensions, and 
there are examples of triangulated manifolds where the link of a simplex is not a sphere. 
In general, necessary and sufficient conditions for having a manifolds out of a simplicial 
complex require that the link of each cell has the homology of a sphere, and that the link of 
every vertex is simply-connected, (see e.g., Thurston's notes [|r§|). Since we are interested 
in dimension 2 < n < 4, we can simply disregard such subtler characterizations. However 
(see the next comments on gluings), this issue is actually bypassed by the particular way 
triangulations are generated in simplicial quantum gravity. 



Gluings. In dynamical triangulations theory one is not given from the outset a triangulation 
of a manifold, rather one generates an n-dimensional PL-manifold by gluing a finite set of 

24 



n-simplices {cr"}. A rather detailed analysis of such gluing procedures is given in Thurston's 
notes ||T^, (clearly for different purposes), and here we simply recall the most relevant facts. 

Given a finite set of simplices and the associated collection of faces, a gluing is a choice of 
pairs of faces together with simplicial identifications maps between faces such that each face 
appears in exactly one of the pairs. The identification space, K, resulting from the quotient 
of the union of the simplices by the equivalence relation generated by the identification maps, 
is homeomorphic to the polyhedron of a simplicial complex. The gluing maps are linear, and 
consequently the simplicial complex K obtained by gluing face-by-face n-simplices has the 
structure of a PL- manifold in the complement of the {n — 2)-skeleton. Since the link of an 
(n — 2)-simplex is a circle, it is not difficult to prove that the PL-structure actually extend to 
the complement of the {n — 3)- skeleton, and that the identification space of a gluing among 
finite ra-simplices is a PL-manifold if and only if the link of each cell is PL-homeomorphic 
to the standard PL-sphere. 

Pseudo-manifolds. It is important to stress that in dimension n > 2, not every simplicial 
complex M = \T\ obtained by gluing simplices along faces is a simplicial manifold, and in 
general one speaks of pseudo-manifolds. 

Definition 1 M = \T\ is an n- dimensional pseudomanifold if: (i) every simplex of T is 
either an n-simplex or a face of an n- simplex; (ii) each {n — l)-simplex is a face of at most 
two n- simplices; (Hi) for any two simplices a", r" ofT, there exists a finite sequence of 
n-simplices cr" = (Tq , cr", . . . , cr" = r" such that cr" and crJYi have an {n — l)-face in common, 
(i.e., there is a simplicial path connecting cr" and r"/ 

Recall that a regular point p of a polyhedron |T| is a point having a neighborhood in |T| 
homeomorphic to an n-dimensional simplex, otherwise p is called a singular point. Absence 
of singular points in a pseudo-manifolds characterizes triangulated manifolds. Moreover, an 
n-dimensional polyhedron |T| is a pseudo- manifold if and only if the set of regular point 
in |T| is dense and connected and the set of all singular points is of dimension less than 
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n — 1. Thus, in the apphcations to gravitational physics, the occurrence of pseudo-manifolds 
is associated with the presence of local irregularities that should not have a sensible effect in 
the continuum limit of the theory, since the set of regular points has the density properties 
mentioned above. 

Dehn-Sommerville relations. In order to construct a simplicial manifold T by gluing 
simplices cr", through their n—1- dimensional faces, the following constraints must be satisfied 

n 

5:(-l)W,(T) = xiT) (38) 

i=0 

if n is even, and 1 < k < n/2. Whereas if n is odd 

with 1 < k < {n — l)/2. These relations are known as the Dehn-Sommerville equations. The 
first, (|38|), is just the Euler-Poincare equation for the triangulation T of which Ni{T) denotes 
the number of z-dimensional simplices, the f-vector of the triangulation T. The conditions 
(PPP or (HDf), are a consequence of the fact that in a simplicial manifold, constructed by 



gluings, the link of every {2k — l)-simplex (if n is odd) or 2A;-simplex (if n is even), is an 
odd-dimensional sphere, and hence it has Euler number zero. 

Note also that in order to generate an n-dimensional polyhedron there is a minimum 
number q{n) of simplices cr" that must join together at a bone, (for instance, in dimension 
2, we have q{2) = 3, otherwise we have no polyhedral surface). In general, since we will 
usually consider only the class of regular simplicial manifolds we assume that the number, 
g(cr*), of ra-dimensional simplices a" which share the subsimplex a* is such that 

q{a') >n-i + l (41) 
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with i < n — 2. 



There is a simple but important relation which must be satisfied by the numbers, 
of n-dimensional simplices sharing the (n — 2)-dimensional subsimplices, the bones of T, 



cr" ^, viz., 



E ^(^ ) = ^^^n (42) 

which follows by noticing that the number of (n — 2)-dimensional sub- simplices in an n- 
dimensional simplex is |n(n + 1). We shall exploit this relation quite intensively, and for 
later convenience, let us introduce the average number of simplices a^, in T, incident on the 
(n — 2)-dimensional subsimplices, cr"~^, namely 

It is easily verified that the relations (|38D , (|39|) , ( ^OD leave — 1, {n even) or |(n — 1) 
unknown quantities among the n ratios Ni/Nq, . . . , Nn/N^, p[. Thus, in dimension n = 
2,3,4, the datum of b{n,n — 2), (trivial for n = 2), and of the number of bones A^„_2, 
fixes, through the Dehn-Sommerville relations all the /-vectors Ni{T) of the dynamical 
triangulation considered. 



Together with the Dehn- Sommerville relations, equation (^2]) implies the following 
bounds for the average incidence number b{n, n — 2), when Nn-2 » 1: 
(i) For n = 2: 

6(2,0) = 6; (44) 

(a) for n = 3: 

3 < 6(3, 1) < 6; (45) 

(Hi) for n = 4: 
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3 < 6(4, 2) < 5. (46) 



Notice however that the lower bounds in (^) and (^6|) are not optimal. This fact was 
noticed during computer simulations (see e.g. |^5[). The actual characterization of the 
optimal lower bounds for 6(n, n — 2) in dimension n = 3 and n = 4 it is rather important 
since it bears relevance to the location of the critical point in simplicial quantum gravity. It 
is not trivial to characterize such bounds since they are a consequence of the following result 
from the combinatorial topology of 3- dimensional PL-manifolds proved by D. Walkup ||55|| : 



Theorem 2 For any combinatorial 3-manifold the inequality 

Ni > ANo - 10 (47) 
holds with equality if and only if it is a stacked sphere 

(a stacked sphere is by definition a triangulation of a sphere which can be constructed from 
the boundary of a simplex by successive adding of pyramids over some facets). 
By exploiting this result we can prove the following 

Lemma 1 For any triangulation T — > M" of a closed n-dimensional PL- manifold M , the 
actual lower hound for the average incidence number b{n, n — 2) is given for n = 3 by 

6(n,n-2)U=3 = ^ (48) 

and for n = 4 by 

6(n,n-2)|„=4 = 4 (49) 



Proof. Let us consider first the three-dimensional case. The Dehn-Sommerville relations 
for 3-manifolds are 

No- Ni + N2- N3 = (50) 
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N2 = 2N3 (51) 

from which we immediately get = Ni — Nq. On the other hand, Walkup's theorem 
imphes Nq < ^Ni + |. Thus, 

iVs > ^iVi - ^ (52) 
which yields the stated result for n = 3. 

In dimension n = 4, the Dehn-Sommerville relations read 

No-Ni + N2-N3 + N4 = x (53) 

2Ni - 3N2 + 4:N3 -5N4 = (54) 

5iV4 = 2Ns (55) 

Let us apply Walkup's theorem to to the link of every vertex v in our 4-dimensional trian- 
gulation. In this way we get 

Ni[link{v)] > 4:o[link{v)] - 10 (56) 

and by summing over all vertices one obtains 

3N2 > 8iVi - lOA^o (57) 

which is the basic inequality that we are going to exploit. From it, (multiplying by 3), we 
obtain 

9N2 > 24A^i - 30No (58) 

There is another independent relation of this sort which can be obtained if we multiply ( |55[ ) 
by 30, and exploit viz., 

30No - 30Ni + 3ON2 - 18Ns = 30x (59) 
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From ([5^) one gets 2Ni — 3N2 + = which, (upon multiphcation by 9), yields I8N3 = 
27 N2 — I8N1. Inserting this latter expression in (|59|) we eventually get 

3N2 = 30x + l2Ni - 30No (60) 

which together with ( pBD provides a direct proof of the bound for b{n, n — 2)\n=4- Explicitly, 
let us subtract (^) to (^8|). We get 

6N2 > l2Ni - 30x, (61) 

namely 

N2 > 2Ni - 5x (62) 
Since 2Ni - ?>N2 + 2N^ = 0, we get 2N3 - 2N2 = N2- 2Ni, thus, according to (|6|), we have 

2N3 - 2N2 > -5x (63) 

And from (^3]) we obtain 

5N^ > 2N2 - 5x (64) 
which, for N2 » x, yields for h{n,n — 2)|„=4 > 4. □ 



It is also interesting to note that according to relation (|4^) the incidence numbers 
{g((T"^^)} of a simplicial manifold cannot be provided by a monotonically increasing se- 
quence of integers, for otherwise Eo-n-2g(cr"^^) would be 0{N'^) in plain contrast with (^2|). 
The actual distribution of the incidence numbers g(cr"~^) of a triangulation has typically a 
peak around a rather small value of g(cr""^) (> 3) and only a few of the bones cr"~^ can 
support a large number {i.e., 0{Nn)) of incident n-simplices. In particular, from ( ^21) its is 
immediate to prove that if we let qmax = ^'^Pt{n,b) suPo-n-2{g((T"^^)} denote the maximum 
value of g(cr"~^), as cr"~^ varies in the set of triangulations with N simplices a" and given 
value of b{n, n — 2), then, for Nn-2 >> 1, 
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_ n(n + l) ( b{n,n-2)-q \ 

2F{n) [ b{n,n-2) j ^^^^ 

where is the number of bones over which q^^x can be attained, and q is the average 
incidence over the remaining bones, i.e., 

E (66) 

''^"-2 F(n)+1 

An elementary example of such a configuration is easily constructed in dimension n — 2 
by considering two distinct vertices Vi and V2 on each one of which N/2 triangles are incident. 
The link of both vertices is an A^/2-gon, and by gluing the stars of Vi and V2 through such 
links we get a triangulation of the sphere with triangles, with q{vi) — q{v2) — y- 

Actually, the existence of this sort of triangulation is a rather general fact of some 
relevance to us. Indeed it provides, for large N^, a triangulation minimizing the average q, 
{q = q+1 = 4 in the above example), with few vertices carrying a large incidence. In general, 
the existence of such triangulations follows from the construction of a triangulation 
5"+^ of the (n+l)-dimensional ball 5"+^ inducing on S'^ ~ 95"+^ a triangulation S" of 

the n-dimensional sphere S" with A^„ simplices and Nn-2 = -^{n — l){n + l)Nn + n bones. ^ 
of the n-dimensional simplices are incident on n distinguished bones {(7*^~^}i=i,...,„ spanning 
an embedded S""^ in T". Thus, for n — 2, {(T*7~^}i=i,...,n consists of two vertices {viz., S° 
which is indeed topologically identified with a double point). For n = 3, {cr*"~^}j=i ... „ 
consists of 3 edges bounding a triangle (~ §^), whereas for n = 4 the supporting bones 
{cr*"~^}i=i^..._„ yields for 4 triangles bounding a tetrahedron (viz., an S^). 

In order to construct such triangulation of consider the (n + l)-dimensional ball 5"+^ 
triangulated by a sequence of A„+i, [n + l)-dimensional simplices {a"'~^^} all sharing a 
common (n — l)-dimensional simplex a^"^ ^ o^, i.e., 

<yT\<y*). <+M = <'(^*)- (67) 

and such that 
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F[a^^\a.), ajtii^*)] ^ ^[^j^i ^Ji' (^*)] = ^* (68) 

where F[aJ~^^ (a^:) , cr]*+i (o"*)] = cr""'"^((j^,) fl cr]^_^i (cr*) is the n- dimensional face through which 
cr^^^(o"*) and aj^i{a^:) are glued together. The boundary of such T"+^ j^n+i ^ 
angulation of the sphere T" with Nn = nNn+i top-dimensional simplices a", and 

Nn-2 = li ~ + l)^n + n bones. The intersection fl T" defines an embedded S"~^ 
in T" generated by n bones, {o"*"~^}j=i_...^„, on each of which ^ simplices cr" are incident. 

These results can be naturally extended to PL-manifolds of arbitrary topology by gluing, 
the above triangulation of the n-sphere into a small triangulation of the given PL-manifod 
M, (say a triangulation of M with a minimal number of vertices). Explicitly, let Tsrnaii{M) 
be a triangulation with few simplices, Nn{Tsmaii) , of an n-dimensional Pl-manifold M of 
given topology, (for a characterization of such type of triangulations and for examples see 
the paper by W. Kiihnel in [l^). Consider the large Nn triangulation T(§") of the n-sphere 
S" defined above, let (t"(M) be a simplex in M and (§") a simplex in S". Set 

T(M) = (TwK^) - ^i{M)) Uf (r(§") - a2"(§")) (69) 

where a{-) denotes the interior of the given cr's, and / is an homeomorphism from the 
boundary of <(M) to the boundary of a^(S"). If iV„(§") » iV„(T,„„;;), such T(M) 
provides a large triangulation of M whose distribution of incidence numbers {q'(o""~^)} 
has the properties of the collection of incidence numbers of T(§"). 

Anticipating a little bit, when this sort of triangulations entropically dominate in the 
canonical partition function (0) we should expect that the resulting statistical system as 
described by (0) exhibits a very large effective dimensionality as ^ 00. A useful clue to 
this is to note that an effective dimensionality on a triangulated manifold can be defined by 
counting the number of vertices visited in a walk of k steps, starting from a typical vertex. If 
the triangulation is sufficiently regular then, for k large, the number of visited vertices grows 
as k^, n being the topological dimension of the underlying PL- manifold. On the other hand, 
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triangulations such as the ones described above are such that all vertices can be visited in a 
few steps: as increases the system effectively behaves as if we were dealing with a regular 
infinite- dimensional lattice. 

Such remarks are elementary, nonetheless they play an important role in what follows 
since the incidence numbers {^'(o""^^} are directly related with curvature when discussing 
dynamical triangulations. Thus there is an intimate relations between the distributions of 
the incidence numbers metric properties, (curvature), and effective dimension of 

the resulting statistical system. 

B. Distinct triangulations of the same PL manifold 

By the very definition of PL manifolds, it follows that there exist different triangulations, 
T«, of the some PL manifold M. For later convenience, it is better to formalize this remark, 
and recall the following standard characterization by W. Tutte [Q: 

Definition 2 Two triangulations, T'^^^ and T*^^^ of the some underlying PL manifold M are 
identified if there is a one-to-one mapping of vertices, edges, faces, and higher dimensional 
simplices ofT^^^ onto vertices, edges, faces, and higher dimensional simplices ofT^"^^ which 
preserves incidence relations. If no such mapping exists the corresponding triangulations are 
said to be distinct. 

Notice that sometimes, {e.g., W. Thurston [|l^,p.l05), such triangulations are said to be 
combinatorially equivalent. However, we shall avoid this terminology since, in simplicial 
quantum gravity, combinatorial equivalence is used as synonimous of PL- equivalence. 

Not surprisingly, in connection to the simplicial quantum gravity program natural ques- 
tions to ask are among the classics of PL- Topology: (i) When can a topological manifold be 
triangulated? (ii) Do any two triangulations of a manifold admit isomorphic subdivisions? 



33 



The answer to the first question is, for general manifolds, negative : Kirby, Siebenmann 



and Wall showed that there exist topological manifolds which do not admit combinatorial 
structures, and that there are non-standard structures for PL-manifolds of dimension at 
least 5. What Kirby and Siebenmann proved |T6| is that for n > 5 there is a single obstruc- 



tion k{M) G H'^{M, Z2), (where if*(M, Z2) denotes the i-th cohomology group of M with Z2 
coefficients field), to the existence of a piecewise linear structure on a topological manifold 
M, and that the vanishing of this obstruction {k{M) = 0) allows for the existence of a 
worth of \H'^{M, 1,2)1 distinct PL structures. Moreover, if n < 7, every PL manifold admits 
a compatible differentiable structure which is unique, (up to diffeomorphisms) , if n < 6. 
The Hauptvermutung is the subject of question (ii), and it is false in general. For instance, 
as hinted before, the 5-sphere admits PL-inequivalent triangulations. In the dimensions 
of interest to simplicial quantum gravity, the situation is more encouraging. Indeed, for 
dimension n < 3 both questions (i), and (ii), have an affermative answer. However, it is not 
yet known if every n- dimensional topological manifold is homeomorphic to a polyhedron 
for n > 4. Actually, (see the paper by T. Januszkiewicz fS^), there is a closed topological 
four- manifold M with the following properties: (i) M not homotopy equivalent to a PL 
manifold; (ii) M is not homemorphic to a simplicial complex. The construction of this ex- 
ample is rather involved and exotic (the basic ingredient being the homology manifold), 
but shows once more, if necessary, that the interplay between the topology of 4-dimensional 
manifolds and quantum gravity cannot be easily tackled. For this reason when discussing 
the 4-dimensional case, we shall not consider explicitly topological 4-manifolds, and work 
directly in the PL-category. When discussing deformations of dynamical triangulations we 
shall need to consider smooth structures, but this is simply a technical issue since as already 
recalled the category of smooth 4-manifolds is equivalent to the category of PL 4-manifolds. 

Another natural question to ask is under what conditions a given set of integers Ni is 
the /-vector of a triangulated polyhedron. This problem is of some relevance to simpli- 
cial quantum gravity, since it can be related to the counting of distinct triangulations, {the 
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Entropy problem). Indeed, a possible way of stating the entropy problem is to enumerate 
all distinct triangulations, on a manifold of given topology, for a given choice of possible 
/-vectors. However, for general triangulated n-manifolds the characterization of all possible 
/-vectors is beyond reach. Partial answers can be obtained only under rather severe con- 
straints. Typically, one has to assume that the polyhedron in question is realized as a convex 
polyhedron in Euclidean d-dimensional space, M'^. In this case, the possible /-vectors are 
characterized by Stanley's theorem , according to which the Dehn-Sommerville relations 
are necessary but not sufficient (they suffice in dimension 2 if we add the obvious condition 
A^o > 4 expressing the fact that to bound a solid takes more than three vertices). In general, 
for a sequence of integers Nq, Ni, . . . , N^-i, set A^_i = 1 and m = [|]. For < p < d, define 



hp by 



and for, 1 < p < m, set 



hp = Ei-^y-' 

i=p yp J 



9p hp hp— I 



i-l 



(70) 



(71) 



then we have 17 



Theorem 3 A sequence of integers Nq, Ni, . . . , N^-i occurs as the number of vertices, edges, 
. . ., {d — 1)- faces of a d-dimensional convex simplicial polytope if and only if the following 
conditions hold: hp = h^-p, for < p < m, (these are just a rewriting of Dehn-Sommerville 
relations); gp> for 1 < p < m; and, if one writes 
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(72) 



with rip > rip-i > 
way), then 



> rij. > r > 1, (there is a unique way of decomposing an integer in this 



9p+i < 
for 1 < p < m — 1. 



^ rip + 1 ^ ( r2p_i + 1 ^ ^ 

+ . . . + 
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(73) 
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This result was conjectured by McMullen |]T8[, sufficiency and later necessity was proven re- 



spectively by Billera and Lee |]T9[, and Stanley and McMullen [|T^, |^0[. Finally, it is worth 
noticing that for the class of all spheres the determination of all /-vectors has been conjec- 
tured by Stanley [^. Owing to the strong limitations imposed by the convexity constraint, 
it is not very convenient to formulate the entropy problem in terms of the characterization 
of the possible /-vectors of simplicial manifolds, and in enumerating distinct triangulations 
we shall follow a different strategy. However, it must be stressed that progress in Stanley's 
theory, (and its unexpected connections with algebraic geometry) may shed light on this 
question. 

C. Dynamical Triangulations 

Let T — ^ M = |T| be a simplicial manifold, (henceforth, if there is no danger of confusion, 
we shall use the shorthand notation M = \T\ for denoting such a triangulated manifold). A 
triangulation T can be used to produce a metric on the underlying PL manifold M = \T\ by 
declaring that all simplices in the triangulation are isometric to the standard simplex of the 
given dimension. This procedure provides a collection of compatible metrics on M, which can 
be extended to a piecewise-flat metric on M in an obvious way, here we follow the elegant 



construction provided by Frohlich |jT5|. Let us choose an orthonormal frame {ea(o"")}a=i 
in the interior of each simplex G T, (this frame can be conveniently located at the 
barycentre of cr'^), and introduce another frame by choosing n directed edges {-E'/^(o"")}j]^]^ 
in the boundary of cr". We can write E^{a^) = J2a=i ^^('^")^a(o""), where the n x n matrix, 
the n-bein, (tj^(cr")) is, for each cr", a regular matrix, a" being a non-degenerate simplex. 

Fixed edge-length metric. In terms of the chosen edges £'^(cr"), we can write the re- 
maining edges in the boundary of cr" as E^y{a^) = i?^(cr") — E,y{a^), with 1 < fi < u < n. 
The metric inside each simplex cr" is assumed to be the usual Euclidean metric, 6. It can 
be expressed in terms of the matrix 5(cr") = g^u^cr^) = {E^,Ei,)s. We assume that each 
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edge E^{a"-) or Ei^,^{a'^) has a, fixed length a — |^^((t")| = |^^,^((7")|, this characterizes fixed 
edge-length triangulations as compared to the variable edge-length triangulations typical of 
Regge calculus. In this way we get a metric in each simplex cr" e M which is extended to 
M by specifying that the metric is continuous when one crosses the faces of the cr^'s, and 
that each face of a simplex is a linear flat subspace of the simplex in question. Thus the 
metric on a fixed edge-length PL manifold is entirely specified by the fixed distance elements 
a between nearest neighbor vertices of the corresponding simplicial manifold, viz., 

Definition 3 For any a e a Dynamical Triangulation Ta of a polyhedron M — \T\, 
with distance cut-off a, and with Nn{Ta) simplices cr", is a division of a simplicial manifold 
M = \T\ generated by gluing, along their adjacent faces, iV„(T(j) equilateral simplices {cr*^} 
with edge-length |^^((7")| = |^^..((7")| = a, V(t" e T. 

Note that in order to avoid overcounting dynamical triangulations differing only by a 
rescaling of the edge-lengths, one always refers to a given value of the cut-off, (typically 
a — 1), when characterizing distinct triangulations. 

D. Dynamical triangulations as length-spaces 

The metric g^i, generated by extending over M — \Ta\ the flat metrics g^u{o-'") of the 
simplices {cr"}, is not the object of primary interest for unravelling the metric geometry of a 
dynamically triangulated manifold. For one thing, this metric is manifestly singular around 
the {n — 2)-dimensional simplices, (which have conical neighborhoods), for the other there is 
nothing sacred in g^^,^ and one can capture more explicitly the geometry of M by looking at 
all curves between given points pi and p2 in M and by setting the distance d{pi,p2) equal to 
the inflmum of the length of these curves. Incidentally one can reconstruct g at any (regular) 
point p from the distance function. Setting f{p) = [d{p,p)]^, it is easily verifled that g is 
provided, to leading order, by the second differential D^f{p)\p=p- 
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In order to discuss some aspects of the interior geometry of a dynamical triangulation 
thought of as a metric space, we introduce few prehminary notions hinting to a characteri- 
zation of geodesies on a dynamical triangulation without using calculus pO . 



A curve in a dynamically triangulated manifold M is a continuous map c: I ^ M = \Ta\, 
where I is an interval. The length L{c) of the curve c: [a, 6] — M is then defined according 
to 

i 

L(c) = sup J2 d{cith^i), c{th)) (74) 

h=l 

where the supremum is taken over all subdivisions a = < ti < ... < tj = b of [a,b], 
and d{, ) denotes the distance function associated with the given dynamical triangulation 
Ta ^ M = \Ta\. We assume that c is rectifiable, {i.e., L{c) < oo), and introduce the arc 
length parametrization as the non- decreasing continuous surjective map 

s:[a,b]^[0,Lic)] (75) 

defined by s(t) = L{c\ [a, t]). In terms of the arc-length, the curve c is travelled at unit speed, 
viz., 

c: [0, L(c)] ^ M (76) 

with c(s(t)) = c{t), and 

L(c|[si, 82]) = \si - S2I (77) 

Similarly we can characterize a curve c: I ^ M travelled at speed f > if the curve is such 
that 

L{c\[ti,t2]) =v\t^-t2\ (78) 

for all ti,t2 G /. If a curve c: I ^ M has constant speed f > and if for any t G / the curve 
c has a neighborhood U in I such that 
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d{c{ti), C(t2)) = V\ti - t2 



(79) 



for all ^1,^2 £ U, then such a c is called a geodesic in the djTiamically triangulated manifold 
M = \Ta\. Such a geodesic is said to be minimizing if ([79| ) holds for all ti,t2 G 

The interior metric dint associated to a dynamical triangulation is defined according 



where inf is taken over all curves from x to y. It is easily verified that for a dynamical 
triangulation such distance function coincides with the original distance function d, dint = d, 
and in such a sense a dynamical triangulation is a Length space It is also a geodesic 



At first sight these facts may seem trivialities, (at least for a dynamically triangulated 
manifold), however they have important consequences when discussing rigorously the way 
a dynamically triangulated manifold approximate smooth riemannian manifolds. A first 
application will be encountered when discussing curvature in the PL-setting. 

PL connections and the Incidence Matrix. In general, associated with a PL metric 
fi'/ii^(c"") = {E^,Ei^)s, o"^ G T, there is also a unique connection (the Levi-Civita connection) 
which can be characterized by the set of matrices, r(c'"", c^i), describing the change of bases 
{-E'^((j")}J]^j^ — ^ {E y,{a'^j^^Yl=\i (regarded as different bases of R"), in passing between two 
adjacent simplices in M, sharing a common face F(cr", o"[^^). In Regge calculus [|T^], the 
Levi-Civita connection so defined is uniquely determined by the length of the edges and by 
the incidence matrix of T, the Nn{T) x A^„(T) matrix I{a^,(j'^), whose entries are 1 if a" 
and 0"" are glued along a common face, otherwise. It follows that when the metric g is 
generated by the fixed edge length prescription, the corresponding Levi-Civita connection is 
uniquely determined by /(cr",(T"). 



to H 



dint{x,y) = inf{L(c)} 



(80) 



space ||5^ since for any pair of points x,y in M there is a minimizing geodesic from x to y. 
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Definition 4 The family ag = {r((Tf , crJYi)}j=i,... with cr" and af'^i adjacent in Ta, is the 
Levi-Civita connection associated with the dynamical triangulation Ta of the PL manifold 
M = \T\. 

Notice that if in place of tlie cliange of bases {E^{a'^)}^^i — > {-^^(cr[Vi)}/i=i consider the 
corresponding change of orthonormal bases {ea(cr")}"^j^ {^a('^IVi)}a=i' passing from 
one simplex to another, then the connection ag can be equivalently expressed in terms of 
orthogonal matrices 0((t", cr^'Yi) given by 

r«, alVi) = t-'«)0(ar, alVi)t(crlVi) (81) 



Simplicial paths and parallel transporters. A trasparent definition of curvature in a 
simplicial manifold can be provided if together with the standard notion of curve on M, 
we introduce the combinatorial notion of simplicial path. A simplicial path in a simplicial 
manifold M is defined as a sequence of simphces {cr"}^^^, / = 2, 3, . . . with a" G M and such 
that and cr"_^^ ^ share a common face. If = then we speak of a simplicial loop 
in the manifold M. With any simplicial loop uj, we can associate the parallel transporter 

K= n r(a;,a;Vi) (82) 

where the product is path-ordered. Let B denote the generic bone in T^, and let uj{B) that 
unique loop winding around the bone B. If <y'j{B) denotes the generic simplex containing 
B, then 



u:{B) = a-iB), . . . , = (83) 

with (Ji{B) n . . . n (T^(-B) = B. The corresponding rotation matrix 

^-(B) =n^(^^(^)'^m(5)) (84) 

i 

is such that all vectors in the {n — 2)-dimensional hyperplane spanned by the bone B are 
eigenvectors with eigenvalue 1, and in the 2-dimensional plane orthogonal to the bone B, 
Ru)(B) reduces to a rotation, where the rotation angle 0(-B) is given by 
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9(B) 

0(B) = E 5,iB) (85) 
i=i 

g(i?) being the number of simplices cr" incident on the bone B, and Sj{B) is the angle 
between the unique two faces of a" containing B. 

Curvature. In Regge calculus, the quantity r{B) = 2tt — 4>{B) is the deficit angle at the 
bone B. Since the angle between two faces sharing a common bone B, in an equilateral 
simplex a", is given by cos~^ l/ra, we can write 

0(5) = q{B) cos-^ ^ (86) 

where q{B) denotes the number of simplices a" incident on the bone B. In defining cur- 
vature for dynamically triangulated manifolds one has to change the perspective on Regge 



calculus somewhat (as advocated by Hamber and Williams ||22|) and view the geometry of 
the triangulation as representing an approximation to some smooth geometry and the local 
curvature as some average curvature for a small volume, (in the next section we formalize 
somewhat this point of view by adopting the Gromov-Hausdorff topology on the space of 
dynamically triangulated manifolds). Thus, to each {n — 2) -dimensional bone B we assign 
an n-dimensional volume density given by 

volniB)^—^ E -^K) (87) 

where {Cj^} is the set of ra-dimensional simplices in the unique simplicial loop uj{B) winding 
around the given bone B. The volume voln{B) can be considered as the natural share of 
the volumes of the n-dimensional simplices to which the bone belongs. It follows that every 
bone 5 in a dynamical triangulation Ta carries a curvature given by 



K{B) 



2tx — q{B) cos 

n 



V0ln{B) 

where voln-2{B) is the usual n — 2-dimensional (Euclidean) measure of the bone B. Recall 
that the Euclidean volume of a j- dimensional equilateral simplex of edge-length a is 
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vol {a- 



(89) 



thus, we get 



K{B) 



n^{n^ — 1) In — I 
a2 i n + l 



2tx - q{B) cos-i ^ 



q{B) 



(90) 



For instance, in dimension four we have voln-2{B) = ^-\/3a^, voln{(y"') = ^"\/5a^, and (88) 



reduces to 



240 /3 
5 



27r - q{B) cos^^ ^ 



(91) 



Curvature assignements. Up to the choice of the cut- off edge- length, a, and of a nu- 
merical factor depending only from the dimension n, curvature in a dynamical triangulation 
T with Nn-2{T) bones -B(a), a = 1, 2, ... , A^„_2(T), is directly provided by the sequence of 
integers {g(a)}, where q{a) denotes the number of simplices a" incident on the bone B{q). 
Henceforth, when speaking of curvature assignements to the bones {B{q)} of a dynamical 
triangulation T^, we shall explicitly refer to the sequence of integers {g(a)}. 



The Einstein-Hilbert action for dynamical triangulations. The Regge version of the 
Einstein-Hilbert action (with cosmological term) for a PL-manifold M, (without boundary), 
is given by §, pj, [pi 



1 



5^e,,e(A,G) = kY.^ol{a-) - ^Y.^(B)vol{B) 



(92) 



When the PL-manifold M is dynamically triangulated — > M, S'/{egge(A, G) simplifies 
considerably reducing to 



S[kn-2, K] = KNn{Ta) - A;„_2 A^„-2 (^a) 



(93) 



where we have set according to standard usage 
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1 COS ^ ~" 

kn = Avol{a^) + -n{n + l)-^Lyol{a''-') (94) 

- ^ (95) 

Thus, for dynamical triangulations the gravitational action reduces to a combinatorial object 
with two (running) couplings /c„_2 and A;„. The former proportional to the inverse gravi- 
tational coupling, while the latter is a linear combination of l/lQnG with the cosmological 
constant A. Accordingly, the formal path integration (|I]) is replaced (see (^) on a dynam- 
ically triangulated PI manifold M, (of fixed topology), by the (grand-canonical) partition 
function 

Z[k^^2,kn]= E 6-'="^"+'="-^^"-^ (96) 
Ter(Af) 

The action S[kn-2, kn] and the associated partition function are so deceiptively simple 
that one may wonder under which rug we are dumping all the notorious problems affecting 
the Einstein-Hilbert action when coming to a path-integral quantization. Obviously most of 
these problems are now shifted in characterizing the summation over distinct triangulations 
J2TeT{M) which replaces the formal path integration over distinct riemannian structures, 
namely in characterizing the dynamical triangulation counterpart of the formal measure 
D[g{M)]. In general, these sort of deals hardly are true bargains. However, as we shall see, 
this trade off in difficulties will pay off since a rather complex problem (Path Integration 
over Riem{M)/Diff{M)), the mathematical formulation of which is rather ambiguous, has 
been transformed into a well-defined counting problem: enumerate all distinct dynamical 
triangulations T G T{M) admitted by a PL manifold of given topology. However, before 
addressing this problem, we have to formalize in more precise terms in which sense dynamical 
triangulations do approximate a riemannian structure on a smooth manifold. 
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V. APPROXIMATING RIEMANNIAN MANIFOLDS WITH DYNAMICALLY 
TRIANGULATIONS OF BOUNDED GEOMETRY 



The fact that curvature for dynamical triangulations is a combinatorial object has a 
number of deep consequences, mainly of topological nature, having an important bearing 
on simplicial quantum gravity. These properties follow from the fact that a dynamically 
triangulated manifold M = \Ta\, with Nn{Ta) ra-simplices and with a given average number 
b{n, — 2) of n- simplices incident on the bones o"""^, (see (0)), is a piecewise-flat singular 
metric space of bounded geometry, a concept which goes back to A.D. Alexandrov |]^. There 
is a definite way in which such singular metric spaces approximate (or are approximated by 
suitably smooth) manifolds, and this will provide the rationale for approximating riemannian 
manifolds with dynamical triangulations. 

A. Dynamical triangulations as singular metric spaces 

To start with an intuitive description, let us consider two geodesic rays starting from 
the some point p in a Riemannian manifold V. Thus, suppose that h: [0, 1] — ^ V and 
l2'- [0, 1] ^ \^ are (minimal) geodesies with /i(0) = ^2(0) = p. Fix t, s G (0, 1) and construct 
the plane Euclidean triangle with (linear-parametrized) lines Li: [0, 1] —>■ M^, L2: [0, 1] — > 
such that Li(0) = ^2(0) = (0,0), ||i:i(s)|| = d{h{0) Ji{s)) , \\L2{t)\\ = d{h{0) , kit)) , and 

\\L,{s)-L2{t)\\ = d{h{s),k{t)). (97) 

If the sectional curvatures of V are everywhere non-negative, then the geodesies h and I2 
tend to come together compared with the corresponding lines Li, L2; i.e., 

||Li(l)-L2(l)|| >f/(/i(l),/2(l)) (98) 

whereas non-positive sectional curvatures force the goedesics to diverge faster than in the 
Euchdean situation: ||Ivi(l) — L2(l)|| < d{li{l) , ^{l)) . 
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Once one has a notion of geodesies, the above characterization of positive or negative 
curvature makes sense for any length-spaces, and in that case we speak of a space of positive 
(or negative) Alexandrov curvature. Actually, similar characterizations can be carried out 
in order to define spaces of Alexandrov curvature bounded above of below by a constant 
A; G M, by comparing the geodesies li, I2 with the corresponding geodesies on the appropriate 
simply connected surface of constant Gauss curvature k. These comparision techniques are 
a standard tool in global riemannian geometry, (see e.g., Chavel's book ||50[| ). 



Geodesic triangles. For A; G R, let S^, denote the model surface oi constant Gauss curvature 
k. Motivated by the above remarks, we characterize dynamical triangulations as singular 
metric spaces with bounded Alexandrov curvature, by comparing suitable triangles in M = 
\Ta\ with triangles in model surfaces S^.. Since in a dynamical triangulation curvature is 
localized around bones a""^, we have to consider geodesic triangles which encircles bones. 
This can be done according to the following procedure. 

Let a = o"""^ be the generic bone in a dynamical triangulation M. Consider its associated 
link, link{cr^^^) ~ and three distinct vertices, labelled 1, 2, 3 in link{a^^^). If q{a) 
simplices a" are incident on the bone a, then link{a^~'^) has a length 

\\link{a''-^)\\ = q{a) ■ a (99) 

a being the edge- length of the given dynamical triangulation. The vertices 1,2,3 are chosen 
(up to a rotation) in such a way that their pairwise distance, d{i,k), i ^ k = 1,2,3,, 
measured along the link link{a^~'^) , satisfies 

d{t,k)>[^] (100) 

where the bracket function [x] stands for the greatest integer not exceeding x. Consider now 
(in the closure of the open star of the bone a) the three geodesic segments pairwise joining 
the vertices so chosen, viz., ci = /(2, 3), C2 = /(I, 3), and C3 = /(I, 2), then 
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Definition 5 The geodesic triangle A(a) = (01,02,03) consisting of the three geodesic seg- 
ments C\, Ci, C3 in M , (called the sides of the triangle), is the geodesic triangle encircling 
the bone a G M. A geodesic triangle A = {ci, 02,03) in is called an Alexandrov (or 
comparision) triangle for A{a) if length{ci) = length{ci), i = 1,2,3. 

It is not difficult to verify tliat sucli a comparision triangle exists, (and it is unique up to a 
congruence), and we have the following 

Lemma 2 A dynamically triangulated manifold M = \Ta\, with Nn{Ta) n-simplices and with 
a given average number b{n,n — 2) ofn- simplices incident on the bones o""~^, is a piecewise- 
flat singular metric space of bounded Alexandrov curvature k, with kmm ^ k < k^ax, where 
kmin, kmax o,re characterized by 



1 ^ -I/-. / N sin(avfiW2;) /ini\ 
— a-gcos [l/n) = ■j== — (101) 



and 



where, (see (^)) 



1 -1/1/ N sinh(av/-fc^i„) 

T^a ■ Qmax cos (1/n) = . (102) 



n{n + l) f b{n,n-2)-q \ 

Qmax = x —T/ ^TT— ^'n- (103) 



2F{n) \ b{n,n-2) 



Proof. We start by showing that comparision triangle exists for any bone of the given 
dynamical triangulation. In order to prove this, let us ffist remark that the geodesic triangle 
A (a) has sides satisfying the triangle inequality 

length{ci) + length{ci+i) > length{ci+2) (104) 

indices taken modulo 3. However, this is not sufficient to insure the existence of a comparision 
triangle in since, if A; > 0, the perimeter of A(q;) must satisfy the bound 
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27r 

length{ci) + length{c2) + length^c-s) < —j=^ (105) 

(this bound is trivially extended to the case < by setting = oo for A; < 0). To prove 
that ( |105| ) holds, let x be any point in the relative interior of the given bone a, and consider 
the set of all tangent vectors to the geodesic rays (lines) emanating from x, of length a, 
pointing into cr", (the generic simplex incident on a), and orthogonal to a. This set of 
vectors forms a spherical (S^)-simplex. Since this set is, up to an isometry, independent 
of X G a, we shall denote this set as S^link{a,a^): the spherical link of a in cr" [pO| . 
Geometrically, the spherical link S^link{a, a"') can be thought of as the arc of circumference 
circumscribed to link{a, a"), the link of the bone a in cr". 

If the Regge curvature around the given bone is positive, then 

q{a) < (106) 

cos ^(1/n) 

Since cos~^(l/?2) > | for n > 2, we get that q{a) < 6. The union of the spherical links 
S^link{a, o"") over all q{a) simplices a"' incident on a is then a circumference of length 
a ■ q{a) cos^^{l/n) < 27ia. Consider now the sphere, S^(r), of radius r, such that a great 
circle of radius a on S^(r) has a circumference of length a ■ q{a) cos^^(l/'n,). An elementary 
computation characterizes the radius r of such §^(r) in terms of q{a) as 

— a ■ g(a) cos~"'^(l/n) = r sin(-) (107) 
27r r 

or, by introducing the curvature A; = of §^(r) 



1 / N 1 / , N sinfavA;) 

—a-q{a)cos-\l/n) = — ^—=^ (108) 
27r y/k 



Since the geodesic triangle A(a) is inscribed in the circumference of length a-q{a) cos ^{1/n), 
we get that (|105|) holds with k given in terms of q{a) by (|108D . By developing ( |108| ) to leading 



order in a, we get, not unexpectedly, that k is given by 



A; ~ — (27r-g(a)cos"^(l/n)) (109) 



a2 
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namely, up to a numerical factor, by the Regge curvature associated with the bone a. This 
elementary analysis can be extended to the case of negative Regge curvature obtaining 

1 , X -X,. I X sinh(av/^) 

— a ■ cos (I/tt.) = -j= (HO) 

27r V — 



thus proving that a dynamically triangulated manifold M = \Ta\, with Nn{Ta) n-simplices 
and with a given average number b{n,n — 2) of n- simplices incident on the bones cr"""^, 
admits comparision triangles on model surfaces, Sjt, of constant curvature k, where k is 
provided by (|108| ) and ( |110| ). Since on such set of djTiamical triangulations the possible 
incidence numbers are bounded between 

q < < Qmax, (111) 

(see (|65|)), the above argument shows that a dynamical triangulation M = |T(j|, as a singular 
metric space, is an Alexandrov space with bounded curvature k^in < ^ < k^ax where kmin 
and kmax are respectively provided by 



1 - -1/1 / \ ^^^{(^vKnax) /110^ 

— a-gcos {l/n) = y== — (112) 



and 



1 _i sinh(aV-fcmm) 

—a-QmaxCOS {l/n) = , (113) 



as stated. □ 



The notion of bounded geometry associated with the characterization of dynamical tri- 
angulations as Alexandrov spaces becomes especially interesting if one considers the set of 
all dynamical triangulations satisfying the hypotheses of lemma |^. Guided by the intuition 
of what happens in similar circumstances for ordinary riemannian manifolds ^ , , we ex- 
pect that the metric space structures of such dynamical triangulations is controlled by such 
curvature bounds. To formulate this kind of controlled behavior, it is useful to introduce a 
topology in the set of dynamical triangulations which arises directly from the classical idea 
of Hausdorff distance between compact sets in a metric space P] . The basic rationale is that 
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given a length cut-off e, two metric spaces are to be considered near in this topology, (one is 
the e-Gromov- Hausdorff approximation of the other), if their metric properties are similar 
at length scales L > e, namely if there is a way of fitting them in a larger matric space so 
that they are close to each other. To formalize this intuition, we introduce the notion of 
rough isometry. 

e-rough isometries. Consider two compact (finite dimensional) metric spaces Mi and M2, 
let dMi{-,-) and dMzi-,-) respectively denote the corresponding distance functions, and let 
0: Ml M2 be a map between Mi and M2, (this map is not required to be continuous ). If 
is such that: (i), the e-neighborhood of 0(Mi) in M2 is equal to M2, and (ii), for each x, 
y in Ml we have 

\dMiix,y) - dM2i<Pix),(f){y))\ < e (114) 

then (J) is said to be an e-Hausdorff approximation. The Gromov-Hausdorff distance between 
the two metric spaces Mi and M2, da {Mi, M2), is then defined according to 

Definition 6 dciMi, M2) is the lower bound of the positive numbers e such that there exist 
e-Hausdorff approximations from Mi to M2 and from M2 to Mi. 



Alternatively [^], ddMi, M2) < e if there is a metric, d{-, ■), on the disjoint union Mi 
extendings the metrics on Mi and M2 and such that Mi C {x G M1UM2: d{x, M2) < e}, 
M2 C {a; G M1UM2: d{x, Ml) < e}. Namely, two compact metric spaces Mi and M2 are 
Gromov-Hausdorff nearby if we can find metrics on the disjoint union of Mi and M2 such 
that Ml and M2 look like the same when they are close to each other. 

The notion of e-Gromov-Hausdorff approximation is the weakest large-scale equivalence 
relation between metric spaces of use in geometry, and is manifestly adapted to the needs 
of simplicial quantum gravity, (think of a manifold and of a simplicial approximation to it). 

Notice that do is not, properly speaking, a distance since it does not satisfy the triangle 
inequality, but it rather gives rise to a metrizable uniform structure in which the set of 
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isometry classes of all compact metric spaces, is Hausdorff and complete. This enlarged 



space does naturally contain riemannian manifolds and metric polyedra. As stressed in [23 



the importance of this notion lies not so much in the fact that we have a distance function, 
but in that we have a way of measuring when metric spaces look alike. 



B. Approximating Riemannian manifolds through dynamical triangulations 

In this section we prove that, dynamical triangulations T^, with a given number A^„(Ta) 
of n-simplices and a given ratio b{n,n — 2), are in the Gromov-Hausdorff closure of the 
cxD-dimensional compact metric space characterized by the following 

Definition 7 For r a real number, D a positive real number, and n a natural number, let us 
define the associated space of Bounded Geometries, TZ{n,r, D,V) , as the Gromov-Hausdorff' 
closure of the space of isometry classes of closed connected n-dimensional riemannian man- 
ifolds {M,g) with Ricci curvature bounded below by r, viz., 

inf {mf{Ric{u,u):u G T^M, \ux\ = 1}} > r 

and diameter bounded above by D, viz., 

diam{M) = sup duiPi^) ^ D 

(p,q)eMxM 



It is rather immediate to recast this characterization in our current setting so as to prove the 
following approximation property (actually a statement of the Lebesgue covering property) 
which has far reaching consequences for simplicial gravity: 

Lemma 3 Let VTn{a,b, N) denote the set of n-dimensional, (n = 2, 3, 4), dynamically tri- 
angulated manifolds M = \Ta\, with Nn{Ta) = top- dimensional simplices and given aver- 
age incidence b{n,n — 2), and let kmin be the corresponding lower bound on the Alexandrov 
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curvature, (see lemma \^). 
geometries with 



Correspondingly, let 7l{n,r, D,V) denote the space of bouded 



V = N ■ valid'') 



(116) 



D < Ci[N ■ vol{a'')]^/' 



n 



(117) 



where vol{a"') denotes the (Euclidean) volume of the equilateral simplex cr" of edge-length 
a, and C\, is a constant (independent from the triangulation considered). Then for any 
n = 2,3,4, any finite b{n,n — 2) G Q^, and G N, there is an e = e(a) > such that for 
any riemannian manifold V G TZ{n, r, D, V) we have 



for some dynamical triangulation M = \Ta\ G VTn{a,b, N) . 

In other words, the set of dynamical triangulations VTn{a,b, N) uniformly approximates 
riemannian manifolds of bounded geometry. 

Note that we restrict the dimension n to the values n = 2, 3, 4, because in such a case, the 
Dehn-Sommerville relations allow us to control the /-vector of the triangulation in terms of 
the number of top-dimensional simplices Nn{Ta) = N and of the average incidence number 
b{n, n — 2). By giving the appropriate combinatorial data, the extension to n > 5 does not 
present particular difficulties. In order to prove lemma ^, we need to recall some standard 
results related to the compactness properties of 7l{n, r, D, V). 

A sequence of compact metric spaces {Mj} converges to a compact metric space M 
provided that dciMi, M) asi oo. Thus we have to prove that there is an e = e(a) > 
such that the Gromov- Hausdorff distance between the generic riemannian manifold M G 
TZ{n,r, D,V) and a suitably chosen dynamical triangulation T^, is lesser than or equal to 



dG{V,M=\Ta\)<e{a) 



(118) 
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e. Note that the e involved is a function of the cut-off a entering the definition of T^. The 
proof of the above lemma is a rather immediate consequence of the following well-known 



properties |T0[: 



Theorem 4 Let M G TZ{n, r, D, V) be the generic riemannian manifold of hounded geome- 
try, and let dM{-,p) denote the corresponding distance function from a chosen point p G M. 
Then, for any given e > 0, it is always possible to find an ordered set of points {pi, . . . ,Pm} in 
M, so that: (i) the open metric halls, (the geodesic halls), BMiPi,^) = {x & M\d{x,pi) < e}, 
i = l,...,m, cover M; (ii) the open halls Bni{pi,^), i = l,...,m, are disjoint, i.e., 
{pi, . . . ,Pm} is a minimal e-net in M. 

We recall that the filling function of a minimal e-net is the number m of points 
{pi, . . . ,Pm}, while the first order intersection pattern of a minimal e-net in M is the set of 
pairs {{i,i)\i,i = 1, . . . ,m; B{pi, e) nB{pj,e) ^ 0}. 

It is important to remark that on TZ{n, r, D, V) neither the filling function nor the inter- 
section pattern can be arbitrary. The filling function is always bounded above for each given 
e, and the best filling, with geodesic balls of radius e, of a riemannian manifold of diameter 
diam{M), and Ricci curvature Ric{M) > {n — l)H, is controlled by the corresponding filling 
of the geodesic ball of radius diam{M) on the space form of constant curvature given by H, 
the bound being of the form |To| < N{n, H{diam{M)f, {diam{M)) / e) . 



The multiplicity of the first intersection pattern is similarly controlled through the ge- 
ometry of the manifold to the effect that its average degree, {i.e., the average number of 
mutually intersecting balls), is bounded above by a constant as the radius of the balls defin- 
ing the covering tend to zero, (i.e., as e ^ ). Such constant is independent from e and can 



be estimated in terms of the parameters n, and H (diam{M)) . 



Definition 8 Two metric spaces (e.g., two riemannian manifolds) Mi and M2 in 
TZ{n, r, D, V) are called e-equivalent ( or e-rough isometric) if and only if they can he equipped 
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with minimal e-nets {pi, . . . ,Pm}i o-nd {j9i, . . . ,Ps}2; respectively, such that s = m and with 
the same intersection pattern. 

We can introduce a natural minimal geodesic ball covering associated with a dynamical 
triangulation according to the following 

Lemma 4 Let (M = |Ta|,Ta) denote a dynamically triangulated manifold with fixed edge 
legth a. With each vertex pi = cr° belonging to the triangulation, we associate the largest 
open metric ball contained in the open star ofpi. Then the metric ball covering of M = \Ta\ 
generated by such balls {Bi} is a minimal geodesic ball covering. It defines the geodesic ball 
covering associated with the dynamically triangulated manifold {M,Ta). 

It is immediate to see that the set of balls considered defines indeed a minimal geodesic 
ball covering. The open balls obtained from {Bi} by halting their radius are disjoint being 
contained in the open stars of {pi} in the baricentric subdivision of the triangulation. The 
balls with doubled radius cover {M,T), since they are the largest open balls contained in 
the stars of the vertices {pi} of T^. 

We are now in a position for proving lemma ^. 

Proof. Given a riemannian manifold M G TZ{n, r, D, V), endowed with with a minimal 
e-net {gi}j=i,...,m we can associated with it the dynamical triangulations T^, with a = a(e), 
having the same combinatorial intersection pattern {pi = c'"°}j=i^...^m and such that \d{qi, qj) — 
d{pi,pj)\ < e for all i, j = 1, . . . ,m. On the disjoint union MjJTa we can define the metric 
d{x, y) = min{(i(g, qi) + e + d{pi,p): i = 1, . . . , m} for q G M, p G T^. More exphcitly, this 
metric is constructed by declaring d{pi,qj) = e, and by setting, for p G M and q G Tq, the 
distance d{p, q) equal to the infimum of the sums d{p, xi)+c/(xi, ^2) + . . . + d{xi, q) where each 
consecutive pair (xj, Xj+i), (or (p, xi), or (x^, q)) either has both elements in one of M or or 
has Xi = Pm, Xi+i = qm- The metric so defined is such that M C {x G M ]J T^: d{x^ To) < e}, 
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Ta C {x & M YiTa- d{x, M) < e}, and thus dG{M,Ta) < 2e, which imphes the stated 
approximation result. □ 

Since a given intersection pattern may be common to distinct dynamical triangulations, 
there can be many distinct dynamical triangulations that are 2e-rough isometric to a given 
riemannian manifold. This is not surprising, the equivalence between riemannian manifolds 
associated with minimal e-nets is very crude, and there are distinct riemannian manifolds 
sharing the same e- net, {i.e., which are 2e-rough isometric). The point is that all such 
manifolds as well as all such dynamical triangulations have a Gromov-Hausdorff distance 
between each other which is smaller than or equal to 2e, and as e ^ 0, viz., as the dynam- 
ical triangulations become finer and finer such dynamical triangulations converge to a well 
defined metric space in 7l{n, r, D, V). 



It is worth noticing that it is possible to obtain a converse of this result, indeed ||T0[, any 
polyhedral compact subset of a Euclidean space is the Gromov-Hausdorff limit of a sequence 
of closed Riemannian manifolds. 

The density result just obtained provides the mathematical rationale for using dynamical 
triangulations in simplicial quantum gravity. The original idea was related to the intuition 
that fixed-edge length triangulations, T^, (of fixed topology M = \Ta\) provide a grid of 
reference manifolds laid upon the cxo-dimensional space of riemannian structures allowed on 
such M. The compactness results associated with the use of Gromov-Hausdorff convergence 
support this picture, and the uniform approximation result of lemma ^ makes precise this 
intuition by proving that dynamically triangulated manifolds do provide a reference grid as 
long as we look at riemannian manifolds and dynamical triangulations of bounded geometry. 
At first sight it may appear surprising that no strict topology control is required on such 
triangulations, but this is actually a bonus coming automatically along with the requirement 
of bounded geometry, {i.e., fixed A^„ (Volume) and fixed average incidence b{n,n — 2) (cur- 
vature bounds)). As we shall see in the next section these natural bounds severely control 
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topology, by avoiding such patologies as infinite genus surfaces (in dimension n = 2), or wild 
homotopy types, etc. One may question the naturality of such bounds in the sense that 
dynamical triangulations is an approach to quantizing gravity in the sense of a statistical 
field theory, and in seeking the associated critical phenomena we eventually need to go to 
the infinite-volume limit. In this limit the bounds on Nn and b{n,n — 2) are removed and 
one consequently loses control on everything. The point we wish to make is exactly that 
the above results tell us how to carry out such sort of limits. One must rigorously establish 
estimates for dynamical triangulations, the associated partition functions and m- points cor- 
relation functions, in a given space of bounded geometry PT„(a, b, N) C TZ{n, r, D, V), and 
then consider a nested sequence X'T„(ai, bi,Ni) C VTn{a2, &2, ^2) C ... C VTn{ai, bi,Ni) . . . 
in which the convergence properties of the partition function and the associated infinite- 
volume limit can be safely discussed. As we shall see, the actual characterization of the 
infinite volume limit does not call for such a sophisticated machinery. However the discus- 
sion of a possible continuum limit may do require the introduction of a suitable projective 
limit procedure for DT„(a, b, N). 



C. Topological finiteness theorems for Dynamical Triangulations of bounded 

geometry 

Since dynamical triangulations, (for given and b{n,n — 2)), are in a Gromov space of 
bounded geometries, TZ{n, r, D, V), we can control to a considerable extent their topology in 
terms of the parameters n, r, D, and V, (namely in terms of the number of top-dimensional 
simplices and b{n,n — 2), in the range of dimension considered). We start with a result 
expressing finiteness of homotopy types of dynamically triangulated manifolds of bounded 



geometry |10 



Theorem 5 For any dimension 2 < n < 4, the number of different homotopy-types of 
dynamically triangulated manifolds realized in VTn{a, b, N) C TZ{n, r, D, V) is finite and is 
a function of n, V^^D"-, andrD^. 
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(Two polyhedra Mi = and M2 = \Ta{2)\ are said to have the same homotopy type if 

there exists a continuous map (p of Mi into M2 and / of M2 into Mi, such that both / ■ </> 
and (p ■ f are homotopic to the respective identity mappings, Imi and ■ Obviously, two 
homeomorphic polyhedra are of the same homotopy type, but the converse is not true). 

This homotopy finiteness property of dynamically triangulated manifolds is an obvious 
adaptation of a more general result proved in a truly remarkable series of papers by P. 



Petersen, K. Grove, and J. Y. Wu [rU|. Actually, for dynamical triangulations, the above 
theorem can be sharpened so as to provide (rather weak) conditions under which two distinct 
dynamically triangulated manifolds share a common homotopy type, and we can estimate 
the number of distinct homotopy types. We have the following 

Theorem 6 For any dimension 2 < n < 4, and for a given value a of the cut-off edge- 
length, the number of different homotopy-types of dynamically triangulated manifolds realized 
in VTn{a, b, N) C 7l{n, r, D, V) is bounded above by 

iHomotopy[VTr,{a, b, N)] < [C{n)-^ ■ Nvol{a'') ■ (119) 

where C{n) is some universal constant depending only on the dimension n, and vol{a^) is 
the euclidean volume of the simplex a" (of edge-length a). 



Proof. Recall |T^, ||T0| that a continuous function ip: [0, a) H^, a > 0, with ip{0) = 0, 



and tp{e) > e, for all e G [0,q;), is a local geometric contractibility function for a (finite 
dimensional) compact metric space M if, for each x E M and e G (0,a), the open ball 
B{x,e), of radius e centered at x G M, is contractible in the larger ball B{x,ip{t)), (which 
says that in M small metric balls are contractible relative to bigger balls) . On a dynamical 
triangulation Ta G VTn{a,b, N), consider the smallest open ball containing the star of the 
generic vertex. This is a ball or radius a which is contractible within itself. Thus, for any 
such Ta we may consider the local geometric contractibility function given by 

V'(e) = e (120) 
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on [0, a]. We can now exploit in our setting a result of P.Petersen to the effect that any 



two dynamical triangulations and in VTn{a,b, N) with Gromov-Hausdorff distance 
dcirl^Tl) < e* are homotopy equivalent if e* ~ a/(32n^), (this bound is not optimal, and 
similarly to what happens for riemannian manifolds with criticality radius bounded below, 
one should get a sharper bound of the form a/[25(n + 1)] [|r^). Since \Ta\ can be covered by 
the open metric balls of radius a centered on the vertices of T^, it also admits a refinenement 
of this covering generated by open metric balls of radius e*/2. If N{e* /2) denotes the number 
of such balls, then an argument due to Yamaguchi |]10| can be used to show that 'DT„(a, b, N) 



contains less than [A^(e*/2)]^ distinct homotopy types. In order to estimate N{e* /2), notice 
that the volume of a metric ball, B{e* /2), of radius e*/2 in \Ta\ is bounded below by 

vol[B{e*/2)] > q{(T^) ■ vol{a''),*/2 = (n + 1) ■ c{n){e* /2Y (121) 

where q{(J^) = (n+ 1) is the minimum number of n- dimensional simplices sharing a vertex, 
and fo/((7")e./2 = c(ri)(e*/2)" is the euclidean volume of a standard euclidean simplex a" 



with edge-length e*/2, (this characterizes the constant c{n) as c{n) = + l/(n!v2")). 
Since the volume of \Ta\ is given by ■ vol{a"'), we get 

Nie* 12) < ^-^"^(^"^ (122) 

^^^^/^^^ (n + l)-c(n)(eV2)" ^^^^^ 



Since e* ~ (a/32n^), by introducing (|122D in the Yamaguchi estimate we get the stated 

;(n)(n+: 
(64n2)" 



result with C{n) = □ 



A similar result, obtained by adapting to dynamical triangulations results obtained by 
Yamaguchi, Grove, Petersen and Wu |jTO|, allows also to control the Betti numbers, (with 
generic field coeffiecients F), of any M = \Ta\ G VTn{a, b, N). Explicitly we have 

Theorem 7 Let bi{M,F) denote the i-th Betti number with field coejficients F, then 

n 

Y,bi{M,F) < (n + l)("+^) [Cin)-^ ■ Nvol{a") ■ a-^'Y (123) 

i=l 
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Proof. The proof is an obvious rewriting of a result of P.Petersen, (see the corollary at p. 
393 of P.Petersen's paper |T0|). 



The above results clearly show that when going to the infinite-volume limit, 
limAr^oo,a-^o < N > vol{(T^) ~ coust., (where < N > denotes the statistical average of N 
with respect to a canonical ensemble of dynamical triangulations), more and more distinct 
topological types of manifolds come into play. However the resulting topological complexity 
is not too wild since the number of distinct homotopy types (and/or the Betti numbers) 
grows polynomially with the inverse cut-off' a~^. 

By considering separately the dimension n = 2, n = 3, n = 4, we can obtain a more 
specific topological finiteness theorem. This result follows again from the local geomet- 
ric contractibility which by construction characterizes dynamically triangulated manifolds. 
Since sufficiently small distance balls on a dynamically triangulated manifold of bounded 
geometry are always contractible, we get, by specializing a theorem of Grove, Petersen, Wu 

m 



Theorem 8 The set of dynamically triangulated manifolds in PT„(a, b, N) C TZ{n, r, D, V) 
contains 

(i) finitely many simple homotopy types (all n), 
(a) finitely many homeomorphism types ifn = 4, 

(Hi) finitely many diffeomorphism types if n = 2 (orn>5, in case we remove the constraint 
2<n<A). 

Note that quantitative estimates associated with these finiteness results have the same struc- 
ture of the ones associated with theorems |^ and |^, and can be easily worked out along the 
same lines. Note also that finiteness of the homeomorphysm types cannot be proved in 
dimension n = 3 as long as the Poincare conjecture is not proved. If there were a fake 
three-sphere then one could prove |]r3| that a statement such as (ii) above is false for n = 3. 



Finally, the statement on finiteness of simple homotopy types, (which actually holds in any 

58 



dimension), is particularly important for the applications in quantum gravity we discuss in 
the sequel. Roughly speaking the notion of simple homotopy is a refinement of the notion 
of homotopy equivalence, relating (in our case) triangulated manifolds (more in general CW 
complexes) which can be obtained one from the other by a sequence of expansion or col- 
lapses of simplices, (for details see [T^). Roughly speaking, it may be thought of as an 



intermediate step between homotopy equivalence and homeomorphism. 

Topology and enumeration of triangulations. The fact that for a dynamical triangu- 
lations with a given number, N, of n-dimensional simplices a", and a given average bone 
incidence b{n, n—2) we have a good a priori control of topology is the basic reason that allows 
us to enumerate distinct dynamical triangulations. An important step in this enumeration 
is to understand how we can reconstruct the triangulation by the knowledge of easy acces- 
sible geometrical data. At this stage it is worth recalling that the Levi-Civita connection 
for a dynamical triangulation is uniquely determined by the incidence matrix /(cr", a") of 
the triangulation itself, thus in order to enumerate distinct dynamical triangulations we can 
equivalently characterize the set of distinct Levi-Civita connections Ug = {r((Tj", crj'Yi)}«=iv- 
The problem we have to face is thus reduced to a rather standard question, familiar in 
gauge theories, where one needs to reconstruct the holonomy of a gauge-connection from 
data coming from the Wilson loop functionals. Here, we need to reconstruct all discretized 
Levi-Civita connection from Wilson-loop datas which reduce to the possible set of curvature 
assignements to the bones, and by taking due care of the non-trivial topological information 
coming from the moduli of locally homogeneous manifolds, {e.g., constant curvature met- 
rics), which need to be described in terms of deformations of dynamical triangulations. In 
order to address this problem we exploit the holonomy representation associated with the 
Levi-Civita connection ag. 
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VI. MODULI SPACES FOR DYNAMICALLY TRIANGULATED MANIFOLDS 



Let {Ta, M — \Ta\) be a dynamically triangulated manifold. If we denote by Qa-{M) the 
family of all simplicial loops in M, starting at a given simplex (Jq, it easily follows that 
fla{M), modulo ctq -based homotopic equivalence, is isomorphic to the fundamental group 
of M, 7ri((TQ,M), based at the given simplex iTq. Moreover, by factoring out the effect of 
loops homotopic to the trivial Uo-based loop, the mapping 

R.u eQaiM) ^ (124) 

yields a representation of the fundamental group 7ri((7o,M). The following definition spe- 
cilizes to the PL setting some of the well known properties of the holonomy of riemannian 
manifolds. 

Definition 9 The holonomy group of the (fixed edge length) PL manifold (M, T^) at the 
simplex aQ, HoI{(t'^) is the subgroup of the orthogonal group 0{n) generated by the set of all 
parallel transporters along loops of the simplicial manifold T based at a^. The subgroup 
obtained by the loops which are homotopic to the identity, is the restricted holonomy group 
of{M,T) ata^, Hol%a^). 

If we change the base simplex ctq to the simplex ctq and fix some path p from ctq to ctq, 
then Hol{aQ) = ApHol{aQ)A~^, (and similarly for the restricted holonomy group). As a 
consequence the holonomy groups at the various simplices of are all isomorphic. And 
we can speak of the holonomy and restricted holonomy group of the (fixed edge length) 
PL manifold M — \Ta\: Hol{M), and Hof{M), respectively. Moreover, if uq and ui are 
simplicial loops based at ctq which represent the same element, [a;], of 7ri(M, ctq), then 
R^Q and R^^ belong to the same arc component of Hol{M). Prom this observation and 
from the definition of restricted holonomy, it naturally follows that there exists a canonical 
homomorphism 
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which defines the Holonomy Representation of 7ri(M). 

Holonomy representations of this sort are of relevance to dynamical triangulations. This 
connection is very subtle and it is related to the well-known fact that generally we cannot 
triangulate locally homogeneous spaces with equilateral flat simplices. For instance, in order 
to model a space of constant curvature with a dynamical triangulation we have to slightly 
deform the triangulation (following a suggestion anticipated many years ago by Romer and 
Zahringer This deformation procedure is apparently trivial in the sense that a small 

alteration of the chosen edge-length allows one to flt the triangulation on a constant curvature 
background. However subtle problems arise if there are non- trivial deformations {i.e., 
moduli) of the underlying constant curvature metric, {e.g, distinct flat tori). Formally, these 
non-trivial deformations are described by the cohomology of the manifold with coefficients 
in the sheaf of Killing vector fields of the manifold, namely with local coefficients in the 
(adjoint representation of the) Lie algebra g of the group G of isometrics considered. This 
cohomology is isomorphic to the cohomology of the fundamental group 7ri(M) with value 
in the holonomy representation of vri(M) in q, and this is one of the basic mechanisms that 
calls into play topology in dynamical triangulations. Such issues can be very effectively dealt 
within the framework of deformation theory in algebraic geometry, (as hinted in the above 
formal sheaf-theoretic remarks), however, we think appropriate here a more pedagogical 
approach, also because this issue has not really been considered previously in dynamical 
triangulation theory. 



A. Romer- Zahringer deformations of dynamically triangulated manifolds 

As a warm up after the above introductory remarks, let us consider in particular the flat 
Levi-Civita connections. Since the connection is flat, the parallel transport along a closed 
simplicial loop u based at ctq depends only on the homotopy class of cu. In this case the 
parallel transport, through its associated holonomy, gives rise to a representation of the 
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fundamental group 

6: 7Ti{M, a^) Hol{M) C 0{n) (126) 

the image of which is the holonomy group of the flat Levi-Civita connection considered. Such 
^ is a homomorphism of '7ri(M, cTq) onto Hol{M), and since 7ri(M) is countable, Hol{M) is 
totally disconnected. In this case one speaks of 6 as the holonomy homomorphism. 

It is a well-known fact that connected compact flat n-dimensional riemannian manifolds 
are covered by a flat n-torus p: T" M, and that the associated group of deck transfor- 
mations is isomorphic to Hol{M), viz., M = T'"'/Hol{M). Quite surprisingly, any finite 
group G can be realized as the holonomy group of some flat compact connected riemannian 
manifold, (this is the Auslander-Kuranishi theorem, see e.g. [0]). Accordingly, flat rieman- 
nian manifolds are completely determined by homomorphisms of their fundamental group 
7ri(M) onto the (finite) group G, up to conjugation. The kind of difficulties encountered in 
dynamically triangulating flat riemannian manifolds are particularly clear when discussing 
the geometry of their covering space. 

Flat tori. It is well known that at least in dimension n > 3, it is not possible to build up 
a flat space {e.g., a flat three-torus) by gluing flat equilateral simplices {o""}, this simply 
follows from the observation that the angle cos~^ ^, > 3, is not an integer fraction of 2tt. 
Thus, for n > 3, dynamical triangulations cannot directly describe flat tori. Superficially, 
this does not seem to be a serious problem, since we can deform a few simplices in such a 
way to match with the 27r- flatness constraint, (see next paragraph for the details). However, 
this is hardly a sound prescription since, as follows from their definition, there are distinct 
fiat tori. More formally, let us fix a basis (ei, . . . , e„) of R*^. Then with each n-ple of integers 
{a}, . . . , a") G Z" we can associate a transformation. A, of M" given by A{x) = a; + X]"=i cl-^Gj. 
This results in a transformation group isomorphic to Z" acting properly discontinuously on 
M". Upon quotienting by this action we get a manifold diffeomorphic to the n-dimensional 
torus T*^. Since R"' is thought of as endowed with its standard fiat metric, the resulting T" 
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is likewise flat. However, by changing base to we change the group of transformation 
associated with A, and the resulting flat torus is not necessarily isometric to the original one. 
Thus, the naive deformation procedure suggested above is ambiguous as long as it does not 
specify which particular flat torus one is modelling. In particular, in dimension n — 2, let 
Oi and 02 be the vectors generating the lattice, and assume that Oi is the first basis vector 
of M^, 0,1 = ei, and that 0,2 = {x,y) lies in the first quadrant of M^. Then the resulting flat 
tori are parametrized by the points in the plane region 

{{x, y)\x^ + y^>l,0<x<^,y>Q} (127) 

to the effect that flat tori generated by lattices corresponding to two distinct points of this 
region are distinct. Thus, even in dimension n — 2, where the flatness 27r-constraint can 
be met, it follows that the lattice generated by equilateral triangles, corresponding to a 
dynamically triangulated torus, can only discretize (up to dilation) the flat torus coming 
from exagonal lattices, (x — y — ^): we cannot dynamically triangulate with equilat- 
eral triangles all remaining flat tori, {e.g., all the rectangular tori 0,2 = {0,y)), (obviously, 
this obstruction persists only if we want the triangulation and the lattice to be consistent, 
otherwise we can dynamically triangulate any 2- torus; in higher dimension we completely 
loose this freedom) . 

The difficulties in associating dynamical triangulations with manifolds of large symme- 
tries is not restricted to the case of flat tori just described. Similar problems are encountered 
in associating dynamically triangulations to any riemannian manifold endowed with a ge- 
ometric structure, namely manifolds locally modelled on homogeneous spaces, for instance 
manifolds of constant curvature, basically because such manifolds do not generally admit 
regular tilings by euclidean equilateral simplices. 

Deforming a triangulation. Even if there is no regular way to generate locally homoge- 
neous PL manifolds by gluing equilateral simplices, {e.g., fiat tori), we can always assume 
that in trying to model a locally homogeneous riemannian manifold M with a dynamical tri- 
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angulation T^, this can be done in such a way that the Gromov- Hausdorff distance between 
M and is minimal. As a matter of fact, we can, at least in line of principle, construct 
an approximating sequence {T^}^=i,2,... of triangulated models of M, interpolating between 
M and T^, by using triangulations of M with simplices which are nearly equilateral, (the 
essence of this remark was very clearly suggested by Romer and Zahringer |TB|). With 
each triangulated manifold of this sequence we can associate the edge-length fluctuation 
functional 

where a(T^) = '^"^ is the average edge-length in T^. The functional s(T^) has 

reasonable continuity properties in the Gromov-Hausdorff topology, and since the Gromov 
space of bounded geometries is compact 0, we can always choose the sequence {T^}/x=i,2,... 
in such a way that it minimizes s(T^). 

However, if there are non-trivial (infinitesimal) deformations of the riemannian structure 
M, preserving the local homogeneous geometry, {e.g., inequivalent flat tori, or moduli of 
Riemann surfaces), the Romer-Zahringer procedure may associate distinct moduli to the 
same approximating dynamical triangulation. Thus, the question remains of how to keep 
track of which locally homogeneous manifold is Gromov-Hausdorff approximated by the 
dynamical triangulation considered. This can be done by exploiting the Gromov- Hausdorff 
continuity of the homotopy type of metric spaces of bounded geometry and by generalizing 
the previous remark relating the holonomy of flat connections to the representations of the 
fundamental group, (see eqn. ( |126|) ). 

B. Dynamical triangulations and locally homogeneous geometries 

In order to formalize the above remarks, we need some preparatory material on rieman- 
nian metrics on homogeneous spaces and the way a dynamical triangulation can approximate 
such particular riemannian structures. From the mathematical point of view, most of this 
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material is standard and as a guide, the reader may found profitable to consult the papers by 
W.Thurston |T6|] and by W.M.Goldman |^ where the classical geometry of deformations of 



locally homogeneous riemannian manifolds is beautifully discussed. As far as applications to 
dynamical triangulation are concerned, the point of view developed here and in the following 
paragraphs is quite new and full of deep implications for unravelling the subtle connection 
between topology and simplicial quantum gravity. 

Locally homogeneous geometries. Let G be a Lie group, H C G a closed (compact) 
subgroup, we denote by X = G/H = {gH: g G G} the set of cosets. To each g E G we can 
associate the left translation Tg-.G/H — >■ G/H given by Tg{g'H) = {gg')H which makes G 
act as a transitive Lie transformation group on the homogeneous space X. Such a space X is 
riemannian homogeneous if G/H is endowed with a riemannian metric relative to which Tg is 
an isometry for every g E G. A locally homogeneous riemannian manifold M is a riemannian 
manifold locally modelled on the geometry of a homogeneous space (X, G) . Constant curva- 
ture metrics are a typical example. They are locally isometric either to euclidean space, the 
sphere, or hyperbolic space, depending on wether the curvature is zero, positive, or negative 
respectively. More generally, locally homogeneous riemannian manifolds include all geomet- 
ric structures used to uniformize surfaces or to discuss Thurston's geometrization program 



for 3-dimensional manifolds IE 



To be more precise, let X be a real analytic manifold and let G be a finite dimensional 
Lie group acting analytically and faithfully {i.e. gx = x for all x G X implies g = id) on X. 
We say that a manifold M has a (G, X) structure if we have an open covering {Ui} of M, 
and coordinate charts (pf. Ui ^ X with transition functions 

7,, = (Pi ■ 0-1; n U,) ^ n U,) (129) 

where each 'jij agrees locally with an element of the group G. 

Let M be a (X, G) manifold and M its universal covering. If we fix a base point Xq in 
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M, M can be identified with the quotient of the set of pairs (x, where x E M and /3 is a 
path from xq to x, with respect to the equivalence relation 

{x,P)c^{x*,P*) (130) 

if and only ii x = x* and (3*13^ ~ 1 in 7ri(A'/, xq). The projection p: M ^ M being defined 
by = x. We wish to recall the characterization of the Developing map D: M ^ X 

associated with a locally homogeneous manifold M. The basic idea is that given a {X, G) 
structure on a manifold M it is natural to try to make the charts {(f/i, (pi)} match up exactly 
and compatibly with the topology of M so as to tesselate M (or portion of M) for suitably 
chosen {(C/j, 0^)}, (in our case the {(C/j, will be associated with equilateral simplices or 
suitable deformations of them) . 

The developing map is defined through analytic continuation along a path a. Fix a 
basepoint G M and a chart (C/q, 0o) the domain of which contains the basepoint xq. 
Given a path a: [0,1] M starting at the basepoint subdivide it at points = Q;(io), 
Xi — a{ti), . . . ,Xn = cx{tn), {to = and t„ = 1), in such a way that each subpath is 
contained in the domain of a single coordinate chart (C/j, 0j). Starting with the given 0o we 
adjust 4>i at xi in such a way that it agrees with the value of cpo at Xg. Since = (pi ■ (pj^ is a 
locally constant map into the group G, the adjustement of (pi is easily obtained by replacing 
(f)i{xi) with ipi = 7oi0i(a^i)- By proceeding in this way, the generic k-th adjusted chart is 
obtained by replacing (pk{xk) with 

V'fe = 7oi (3^1)712(2^2) ■ ■ ■lk-i,k{xk)4>k{xk)- (131) 
The adjusted charts ipi,ip2, ■ ■ ■ j'ipn form the analytic continuation 

# = V'n (132) 

of the chosen 0o along the path a. It can be verified that this analytic continuation only 

depends on the homotopy class [a] of the path, in particular it does not depend from the 
open sets {Ui}, but only from the choice of the first open set {Uq, 0o) containing the basepoint 
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Xq. We get in this way a well-defined application D: M — > X. Notice that different choices 
of (f/o, (/)o) characterize applications D differing by a G- translation in X. Formally we have 
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Definition 10 For a given choice of the initial chart {Uq, 0o) one? of the basepoint Xq G Uq, 
the developing map of a locally homogeneous manifold M is the map D: M ^ X defined by 
D = 00°' ■ p in a neighborhood of (x, a) G M. 

Note that in the definition of developing map it is not strictly necessary that the action of 
G on X is transitive. 

Holonomy homomorphism. We can now recall the characterization of the Holonomy of 
a locally homogeneous manifold M, (the holonomy representation ( [I26| ) is a particular case 
of such a construction). Let a be an element of the fundamental group of M, tti{M,xo)- 
To such a loop there is associated a covering transformation M ^ M in the universal 
cover defined by Trj[{x,T)] = [{x,ar)]. Analytic continuation along the loop a provides 
0{r' = ipn which is well-defined at the basepoint xq, (since x„ = xq), and depends only on 
the homotopy class of the loop a. We call Holonomy of M the map H: tti{M) G defined 
by o" 1-^ h{a) where h^r is that element of G such that 

0fl = h^<Po (133) 

Since D ■ = h„ ■ D, the holonomy map if is a group homomorphism from tti{M) into 
the structure group of the geometry G, the holonomy homomorphism. The holonomy is 
determined up to conjugacy in G according to the fact that the developing map is determined 
only up to left composition with elements of G. 

Note that if the developing map is a covering map, M is a complete (G, X) manifold 
and its holonomy group characterizes M, in the sense that M can be reconstructed as the 



quotient space X/Hol{M), (see Thurston's notes |jT6[ for more details) 
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More generally p6[, the developing map D: M X can be used to pull back the (X, G) 
structure from X to M which thus inherits a natural structure of locally homogeneous 
riemannian manifold. Moreover, the holonomy homomorphism H:7Ti{M) —>■ G determines 
the action of tti (M) on the universal cover M by {X, G) automorphism. This implies that the 
pair (Z), H) associated with the developing map D and a holonomy representation H in the 
structure group of the geometry, the so- called Development pair, is a useful globalization of 
a locally homogeneous structure (X, G) defined on a manifold M by local coordinate charts 
the pair {D,H) completely determines the locally homogeneous geometry on M. 

Moduli of locally homogeneous geometries. For our purposes we need to consider 
not just a particular locally homogeneous structure {M,g*) on a compact manifold M, but 
rather a space of such structures. To this end ||2^, mark a basepoint xq E M and let 
7ri(M, Xo) and M ^ M be the corresponding fundamental group and universal covering 
space. We consider the set of homomorphisms of 7ri(M, Xq) into the structure group of the 
geometry G, Hom{TTi, G) endowed with the pointwise convergence topology, i.e., a sequence 
9n G Hom{ni{M), G) is said to converge to a homomorphism 6 if lim^^oo Qnil^) = 6'(ti;) for 
each uj G 7ri(M, xq). Let us consider, as before, the coordinate map (pQ. f/o C M — ^ X, in 
particular its value at xq G Uq. This defines the germ at xq of the locally homogeneous 
geometry of M. 

The space T(x,g)(^) of isotopy classes of based locally homogeneous structures on 



(M, xo) having the same germ (f)o:UoCM^X at xq can be defined as the set 



J(x,G)(M) = {{D,9)\9 G Homi7r,iM),G)} (134) 

and where D: M ^ X is a. ^-equivariant nonsingular smooth map, viz., a developing map. 
If we topologize X(x,g)(^) using the G°° topology on the developing map D, then the map 

hol:I(x,G){M) Hom{7ri{M), G) (135) 

which associates to the developing pair [D, 6) the corresponding holonomy representation 
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6 is continuous. Notice that if Dif fo{M,xo) is the identity component in the group of all 
diffeomorphisms M —>■ M fixing xq, then Di//o(M, xq) acts properly and freely on 
I(x,G){M) and the map hoi is invariant under this action. It is also easily verified that 
hoi is G- equivariant also respect to the natural G-actions of G on the developing map 
and on the representation space Hom{TTi{M), G). In the quoted Goldman's paper the 



intersted reader can find the statement and the proof of an equivariance slice theorem which 
allows to parametrize the space 1(x,g) {^) of based geometric structures on M in terms of 
a representation 9 G Hom{7Ti{M) , M) and of a diffeomorphism h G Diffo{M,xo)), (see the 
deformation theorem at p. 178 in [|26|). We state this result in a somehow less technical 
fashion 

Theorem 9 The map hol:I(x.G){M) H 0171(71 i{M),G) is an open map, and for any 
given locally homogeneous structure u G I(x.G){M) there is a corresponding neighborhood 
W C X(x,G)(^) such that for any two (based) locally inhomogeneous structures Mi,M2 G W 
such that hol{ui) = hol{u2) = 9 E Hom(7ri(M),G) there is a (based) diffeomorphism h G 
Dif fQ{M, xo) and 7 G G with {h, 7) ■ Mi = 'U2- 

In other words, in a sufficiently small neighborhood of a given locally homogeneous struc- 
ture, any two locally homogeneous structures with the same holonomy representations are 
equivalent, and we can use the holonomy homomorphism as a sort of coordinate map for 
the space of inequivalent locally homogeneous structures that we can obtain from the given 
one by local deformation. 

It also follows from the above theorem that the map 

hohi^xMM) = n^S-T/lf ^ ^ ^ Hom{n,{M),G) (136) 
Dtffo{M,xo) 



is a local homeomorphism. If we define |^ the Deformation Space of locally homogeneous 
structures on M as the quotient 

Def.xMM) = ^^^^ (137) 
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then one would expect that the map hoi extends to a local homeomorphism from 
Def(^x,G){M) and the orbit space Hom{TTi{M),G)/G (endowed with the quotient topol- 
ogy). In general this is the case only if we have more information on the structure of the 
orbit space under the G-action. Typically hoi is well-behaved if we restrict our attention to 
the subspace in Hom{TTi{M), G)/G generated by the conjugacy classes of stable representa- 
tions. In any case we shall work mainly with hoi: I(^x,g){M) — > Hom{7ri{M), G), and simply 
assume, at least at this stage, that the representations considered are such that the local 
properties of the representation variety Hom{TTi{M),G)/G describe, through the holonomy 
map hoi, the local properties of the deformation space Def(^x,G)iM) of locally homogeneous 
structures on M. This is an important point for it will imply that the representations in 
Hom{ni{M),G) p^Qyj^jg ^j^g missiug piccc of information needed for describing the inequiva- 
lent deformations of dynamical triangulations approximating a given locally homogeneous 
riemannian manifold. 

By exploiting the preparatory material so introduced, we now describe in details the 
procedure for parametrizing the deformation space of dynamical triangulations. 



C. Moduli of dynamical triangulations 

Consider on a dynamical triangulation (M, r^) the set of metric balls of radius a, 
{Ui{M)}, centered at the baricenters, Pi, of the iV„_2(M) bones {o"""^} of M. Each such 
Ui{M) is homeomorphic to 5""^ x C {link{a^~'^)) , where C {link{a'^~'^)) is the cone on the 
link of cr"""^, and denotes the topological /c-dimensional ball. Since C{link{a'^~^)) ~ B^, 
we get that Ui{M) ~ x ~ B'^, viz. each metric ball Ui{M) is indeed a topological 
ball. Since their radius is a, the collection of such Ui{M), i = 1, . . . , A^„_2(M), generates an 
open covering of M. 

Assume that the dynamically triangulated manifold (M, T^), endowed with the covering 
{Ui{M)}, is Gromov-Hausdorff near to a locally homogeneous {X,G) riemannian manifold 
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(of bounded geometry) (S, g). This assumption can be translated into the requirement that 
the riemannian manifold (S,(yf) admits a geodesic ball covering {gj}i=i,...,m £ S generated by 
metric balls f/j(S), of radius e = a, such that 



since, by a standard argument already exploited various times, this bound implies that the 
Gromov-Hausdorff distance between M and S is such that dciMjE) < 2a. By itself, the 
dynamically triangulated manifold M cannot be endowed, for any fixed value of the cut-off 
a, with a locally homogeneous metric geometry, (with the exception 

of some particular cases in dimension n = 2). As already stressed this may be a source of 
serious ambiguities. In particular, let us assume that the given locally homogeneous rieman- 
nian manifold (S,(yf) admits non-trivial infinitesimal deformations (moduli) Def(^x,G){^)i 
then, up to the cut-off length scale a, the dynamically triangulated manifold (M, To) is 
roughly-isometric to any of the inequivalent locally homogeneous manifolds associated with 
the moduli of (S,5f), i.e., according to ( |138D , 



for any S* G Def(^x,G){M). Thus, (M, Ta) can be isometric to any locally homogeneous 
structure in Def(^x,G){M). In order to avoid this ambiguity, we need to declare which locally 
homogeneous geometry (M, Ta) is actually approximating. 

To enforce such control, let us consider (M, T^) endowed with the covering {f/j(M)}. 
Choose a basepoint po among the baricenters {pi} of the bones of M, and let Uq{M) denote 
the corresponding metric ball in the cover {Ui{M)}. Finally, mark an ra-dimensional simplex 
(To among those incident on the bone o"o~^ associated with the given choice of po- Notice that 
with these markings we can associate to Uq{M) a well- defined element, i?o('^o~^)) of the 
orthogonal group 0{n): _Ro(o"o~^) is the rotation around the bone ctq"^ of the vectors defining 
the marked simplex a^. Together with these elements, let us also choose a representation 



\dM{Pi,Pj) - d^{qi,qj)\ < a. 



(138) 



dG{M,E*) < 2a 



(139) 
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9: 7[i{M,po) — G 3 0{n) of 7ii{M,po) into the structure group G 3 0(n) of the locally 
homogeneous geometry approximated by (M, T^). 

Since (M, T^) is Gromov-Hausdorff near to a locally homogeneous manifold (H,5f), there 
is by hypothesis an open set Uo(E) on S, corresponding to Uo{M), centered around go ^ 
and a corresponding map 0o: f^o(S) {X, G) associated with the element, i?o(c"o~^)5 of the 
orthogonal group describing the rotation around the bone ag^^. If we think of Ro{aQ^^) as 
an element of G ^ 0{n) associated with the chosen marked simplex ctq, we can generate 
00 by analytical continuation of 0(go) = -Ro(cTo~^) in Uo(E) along the unique simplicial loop 
winding around ctq"^- Thus, associated with the markings of (f/o(M),(TQ) on M there is a 
local chart (t/o(H), 0o) for the locally homogeneous geometry of S. Since Gromov-Hausdorff 
nearby metric spaces share a common homotopy type, it also follows that tti{M,pq) ~ 
7ri(H,go); (provided that a is small enough so that M and H are sufficiently near in the 
Gromov-Hausdorff distance), and the representation of iti{M,pq) induces a corresponding 
representation that we still denote by 9, 9: go) — > G. 

The relevance of such interplay between the markings of (f/o(M), ctq ) on M and the cor- 
responding local chart (?7o(H), ^o) on S, is twofold and follows from the properties of locally 
homegeneous geometries recalled in the previous paragraph. First, the local chart (f/o(S), 0o) 
can be used as a basepoint for the development pair {D, 9: vri(H, go) — G) defining the locally 
homogeneous structure (H, g) considered. Moreover, according to the Goldman-Thurston de- 
formation theorem ^, if the geometry (^,(7) admits non-trivial deformations Defi^x,G){^)^ 
such deformations are locally parametrized by the representation variety ifom(7ri(H), G)/G. 
Thus, the data needed in order to specify unambiguosly which locally homogeneous geometry 
we are approximating with (M, T^) is basically contained in the association of a holonomy 
homomorphism 9: 7ri(M, po) — > G to the dynamical triangulation (M, T^). 

Summing up, we have the following 
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Theorem 10 Let (D, 6': 7ri(S, go) G) be the developing pair of a locally homogeneous 
geometry (H,5f) with a non-trivial deformation space Def(^x,G){'^), o-nd let {M,Ta) G 
VTn{a, b, N) denote a dynamical triangulation approximating (H, g). Then the distinct locally 
homogeneous riemannian manifolds in Defi^x,G){^) ^^i"^ Gromov-Hausdorff approximated by 
the pairs [{M,Ta)r ootid] where: (i) {M,Ta)root is the rooting of the dynamical triangulation 
{M,Ta) generated by the marking of a bone o"o ^^ of {M,Ta), of a simplicial loop Uq winding 
around ctq"'^ , and of a base simplex G oJq, while (ii) is a holonomy homomorphism 

e:7ri{M,po) ^GDO{n) (140) 

varying in a suitable neighborhood of 9 in the representation variety ^""^i'^^i^)''^ ^ 



Proof. The markings in (i) are simply the markings characterizing the data (t/o(M),(To) 
on M inducing a local chart {Uo{E), (po) for the locally homogeneous geometry of S, through 
which S can be unrolled by means of the developing map D. These data are necessary for 
providing a well-defined holonomy representation. Henceforth we shall refer to dynamical 
triangulations with such marked simplices as rooted triangulations for short. 



As we have recalled |26], the map 



hoi: Def(x,G){^) ^ Hom{ni{E), G) (141) 

which associates to the developing pair {D,6), {i.e., to the locally homogeneous manifold 
(H,5f)), the holonomy representation ^ is a continuous and open map. Two locally homoge- 
neous manifolds Si and S2 in De/(x,G)(S) are near in the G°° topology (on the respective 
developing maps Di and D2) and the corresponding geodesic ball coverings {f/j(Si)} and 
{Ui{E2)} are such that ^^(21,22) < e where e depends on the C°°-norm on the respective 
developing pairs \\{D,9i) — (-D, ^^2)||oo- It follows that (M, T^) is Gromov- Hausdorff near 
any locally homogeneous manifolds in De/(x,G)(S). 

By choosing, if necessary, a sufficiently small neighborhood of (2,(7) in Def(^x,G)i'^), 
can assume that (M, T^) shares a common homotopy type with all {^*,g) in De/(x,G)(S)- 
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Thus an immediate application of the Goldman-Thurston deformation theorem |] imphes 
that we can estabhsh a one-to-one correspondence between each pair [{M, To) root, and 
a corresponding locally homogeneous manifold in De/(x,G)(2), simply by declaring that a 
given [{M,Ta)root,0] is associated with the locally homogeneous manifold in De/(x,G)(S) 
defined by the holonomy homomorphism 6, and this correspondence is well-defined for 6 
varying in a suitable neighborhood of 9 in -^°™-(^(*^)''^ ^ □ 

Thus, with a slight abuse of language, it seems appropriate to consider the representation 
variety ^j^g Moduli Space for a dynamically triangulated manifold M E 

VT^{a,h,N). 

Since we are interested in closed manifolds M, the study of this moduli space is basically 
equivalent to the study of discrete cocompact subgroups V of the maximal isometry group G 
that can be ammitted by M, (recall that a discrete subgroup F C G is called cocompact if the 
quotient space G /T is compact). The underlying geometric rationale is that the Lie group in 
question is the isometry group of the simply connected homogeneous riemannian manifold 
X which can be approximated by the (universal cover of the) dynamical triangulation M. 
The discrete subgroup F appears (if it has no torsion) as the deck transformation group of 
the covering X S where S = X/T is the smooth riemannian manifold approximated by 
M, thus F ~ 7ri(M), (if F has torsion we will be dealing with orbifolds). It is intersting 
to remark that even if not every semisimple Lie group G may be realized as the isometry 
group of some symmetric space X, every such G is locally isomorphic, {viz. the isomorphism 
is at the level of the universal coverings of the groups in question), to the isometry group 
of some symmetric space X. The existence of a non-trivial moduli space for dynamical 
triangulations (of closed manifolds) is thus reduced to the study of the rigidity of lattices in 
such groups G [^. 

For the convenience of the reader, we have collected in an Appendix some basic material 
on the geometry of the representation variety Here we stress that since vri(M) is 
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a finitely generated group, (say with m generators), Hom{'Ki{M)^ G) is an analytic subvariety 
p6| of G"", (as a variety is defined by the m-tuples (9(1), . . . , 9(m)) G G"^ such that 
pj(0(l), . . . , 0(m)) = 1 for each relation pi associated with the given presentation of the 
group 7ri(M)). Note that on applying a theorem due to X. Rong |Q, one can actually provide 
a upper bound to the minimal number of generators for 7ri(M) for M G VTn{a, b, N), (the 
bound depends on the parameters a,b, and A^). 

D. Gauge-Fixing of the moduli of a Dynamical Triangulation 

Far most of our efforts will be bent toward an enumeration of the distinct dynamical tri- 
angulations a PL manifold can carry. Since distinct dynamical triangulations (of bounded 
geometry) are meant to approximate distinct riemannian structures, we are actually address- 
ing the characterization of a measure on the space of riemannian structures. As we have seen 
in the previous sections, the correspondence between dynamical triangulations and rieman- 
nian structure ceases to be one-to-one nearby locally homogeneous manifolds, no matter 
how fine is the triangulation. Distinct moduli of locally homogeneous riemannian mani- 
folds may indeed correspond, through a Romer-Zahringer deformation, to a same dynamical 
triangulation. In order to carry over an enumeration of distinct dynamical triangulations 
while preserving a one-to-one correspondence between riemannian structures and dynamical 
triangulations, we have have to impose a suitable gauge fixing so as to remove this possible 
source of ambiguity. The situation is akin to the one familiar in gauge theories where, in 
order to introduce a suitable path integration on the configuration space of the theory, we 
have to remove gauge freedom, (by using a slice theorem for the action of the gauge group), 
to obtain a local (formal) measure on the orbit space labelling the gauge equivalent config- 
urations. Eventually, as in any gauge fixing procedure, in order to obtain the correct count 
of the number of distinct dynamical triangulations, we have to divide the number of distinct 
holonomies by the local volume of the moduli space parametrizing the conjugacy classes of 
representation Q: ttiIM; a^) — G in a neighborhood of the chosen representation. From a 
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combinatorial point of view, the necessity of this factorization comes from the fact that the 
number of distinct holonomies provides the number of distinct dynamical triangulation in 
a given holonomy representation of the fundamental group, namely it is a labelled or rooted 
enumeration. The unrooted enumeration can be obtained from the rooted one by dividing 
by all possible labellings of the inequivalent representations. 

A heuristic remark. In connection with this counting problem for dynamical triangula- 
tion, it is worth remarking a subtle point. Had we addressed the counting of just distinct 
(in the sense of Tutte) dynamical triangulations on a PL manifold of given volume and 
topological type, the above gauge fixing problem would have not explicitly appeared. The 
fact is that in the definition of equivalence according to Tutte, no mention is made of the 
metric properties of the triangulation considered, and thus in the process of counting, we 
enumerate also triangulations potentially fitting on locally homogeneous spaces. The re- 
quirement of triangulating with equilateral simplices, is a restriction of no consequence as 
long as we consider generic riemannian manifold, but quite effective when coming to tri- 
angulating riemannian manifolds of large symmetry. We enumerate distinct triangulations 
basically by considering distinct curvature assignements, viz., namely using as labels for 
distinct triangulations the distinct Levi-Civita connections generated by distinc curvature 
assignements. This strategy can be considered as the discretized counterpart of the tech- 
niques underlying uniformization theory for 2-dimensional surfaces or the geometrization 
program for 3-dimensional manifolds. The simplest example is provided by 2-dimensional 
surfaces. In such a case, a given riemannian metric can be conformally deformed to a metric 
of constant curvature. The conformal factor is related (via the coformal Laplacian) to the 
gaussian curvature of the surface, and in this sense the datum of the local curvature seems 
sufficient to reconstruct the original surface. This fails since there are non- equivalent in- 
finitesimal deformations of the base constant curvature metrics (moduli), and in order to 
reconstruct the metric, we need curvature assignements and the datum of which constant 
curvature metric, {e.g., which flat metric, if the surface in question is a 2-torus), is used. 
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Actually, there is a further subtle problem related to the fact that the action of the semi- 
direct product of the group of diffeomorphisms and of the conformal group makes difficult 
an actual reconstruction of the metric from the datum of the local curvature and of the mod- 
uli. But at this heuristic level, the parallel between the geometric setup of uniformization 
theory and the reconstruction of dynamical triangulations from the curvature assignements 
is appropriate. 

A measure on the moduli space. According to the above remarks, in order to 
complete our gauge-fixing procedure, we have to introduce a measure on the moduli space 
Ai{Ta) of a given dynamical triangulation. This procedure is related to the characterization 
of measures on the representation variety -^°™(MA'-?"),g) ^ dimension two, the study of such 
measures has been considerably developed by exploiting its connection with the semiclas- 
sical limit of Yang-Mills theory on Riemann surfaces p9[. For dynamical triangulations in 



dimension greater than two, the connection with measure theory on -'^°™-(""^(^'^)''^) j^^g been 
stressed in ^ as the origin of the possible critical behavior in higher dimensional simplicial 
quantum gravity. 

Let us consider the representation, 6, of 7ri(Af) on the Lie algebra g generated by com- 
posing the given representation 6:'7ri(M) G with the adjoint action of G on g, viz., 
6: 7ri(M) G End{Q), (henceforth we will always refer to this representation ). The 
tangent space to ^°"^('^^^)''^) corresponding to the conjugacy class of representations [9] is 
provided by H^{tti{M), g), the first cohomology group of vri(M) with values in the Lie alge- 
bra g. It is a known fact that this cohomology group is isomorphic to the first cohomology 
group of the dynamically triangulated manifold M = \Ta\, with coefficient in the adjoint 
bundle ad{6) defined by the representation 6, i.e., H^{M,ad{6)), (for more details see the 
Appendix). Notice that H^{M,ad{9)) is only a formal tangent space to -^"™(""^(^^)''^) ^ since 
the representation variety has singularities in correspondence to reducible representations 
and/or linearization unstable representations. 
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The choice of a well-defined holonomy representation amounts to the choice of a map 
(basically a section) M^liliiAMlj^fl _^ Hom{Tii{M),G) which associates to a conjugacy 
class, [9], of representations, the holonomy representation 9: 7Ti{M,pq) g, based at the 
chosen basepoint. 

Let SG denote the generic (infinitesimal) deformations, at 9, of the chosen slice S. Note 
that 6G G Z^(7ri(M); g), the space of 1-cocycle on 7ri(M) with coefficients in the Lie alge- 
bra 0. This is the (formal) tangent space to Hom{Tii{M),G) at the chosen representation 
9. A component of the deformation 5G may be trivial in the sense that it may lie in the 
tangent space to the Arf-orbit containing 9. This part of the deformation maps, by Ad- 
conjugation, S{9) to a deformed representation in a nearby slice S{9) with [9] = [9]. The 
non-trivial part of the deformation lies in the tangent space to the chosen slice S, viz., 
TS [ ^""^i"^^^) ■'^) ] ~ H^(7Ti{M), q). This latter part of the deformation maps 6* in a neigh- 
boring representation which cannot be obtained by Ad-conjugation from 9. As recalled, 
deformations in y5'|^ -^°'^(^^(^^)''^) j qq^j^ interpreted as inequivalent Romer-Zahringer defor- 
mations of the given dynamical triangulation. 



A natural volume element on the space of such deformations is the product 



of a volume form on the Lie algebra g and of a volume on j'g^M^lllil^jMlj^'j^ order to 
define the former, consider an Haar measure on the group G, with total volume Vol{G), 
and let yi,y2, ■ ■ ■ ,yt be Euclidean coordinates on the Lie algebra g. We choose the yh in 
such a way that the measure Hi dyt on q is the chosen Haar measure at the identity element 
of G. Analogously, we can introduce a measure, d^, on the tangent space to the chosen 
representative of the (smooth component of the) variety g^E^l^lilAMli^'j q Hom(iTi{M), G). 
Recall that jg analytic subvariety of G"", m being the number of generators 

of TTi{M). The fixed Ad-invariant inner product on q induces a Riemannian structure on 
G™, and hence on g^M^^EtiEk^Ml^'j^ (^a^ least on its smooth component, to which we are 
restricting our considerations). The corresponding inner product between tangent vectors 
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in H^{M, ad{9)), will be denoted by << •, • >>, and dji can be thought of as the associated 
volume measure. We shall denote by 

V^voli'^^^^^^^Mm^ (142) 

the volume of the moduli space in this measure. We cannot use dii as such since 

we want to localize, as the representation variety becomes smaller and smaller, the choice of 
the representations 9 around the trivial representation. This is dictated by the requirement 
that as V ^ 0, {e.g., for simply connected manifolds), the measures should converge to 
the atom concentrated on the representation [a;,po] ^ Ida, where [a;,po] denotes the trivial 
loops based at the chosen basepoint po- 

To this end, we normalize the product measure d/i x on s x rg[ g»'»(MM),G) j^ 

the heat kernel on the analytic manifold Hom{7Ti{M), G) C G"^, and set 

d^ = K (iV„vo/((T")V, e, Id) dii X IJ^^ (143) 

Vol{G) 

where 6 is the given representation in Hom{TTi{M),G), Id is the trivial representation, 
vol{a"') is the volume of the elementary simplex cr", and 

K{t,x,y):R+ x Hom{7ri{M),G) x Hom{7ri{M),G) ^ M. (144) 

is the fundamental solution of the heat equation on Hom{'Ki{M),G), (since Hom{'Ki{M), G) 
is an analytic subvariety of G x G x . . . x G, K{t, x, y) can be constructed with the heat 
kernel on G). Notice that NnVol{a"')V can be interpreted as the volume of the dynamically 
triangulated manifold M as expressed in terms of the volume of the representation 
variety Hom{'Ki{M),G)/G. 

Thus, the natural normalizing volume element we consider in the procedure of gauge 
fixing the Romer-Zahringer deformations of a dynamical triangulation is 

Vol[M{T^)] = j^^^^u,,,,^^ K {N^vol{a-)V, 9, Id) d/. x (145) 



79 



Moduli asymptotics. If we exploit the standard properties of the heat kernel, we can 
get an explicit asymptotic expression for Vol[/iA{Ta)] in the limit of small volume for the 
representation variety. 

Theorem 11 For a given value of Nn and for Vvol{a"') near zero, the leading contribution 
to Vol[Ai{Ta)] is given by 



Vol[M{Ta)] 



^VlI0i{(T")<<l" 



A{V) ■ {2'KNnVvol{a''))-^/\l + 0(iV„Vt;o/(a"))) (146) 
where D = =dim[H om{TTi{M) , G)], and A{V) is a function ofV 
(The explicit expression of A{V) is provided in the proof). 

Proof. We shall consider Hom{TTi{M), G) C G"" as a compact manifold with a Riemannian 



structure induced by the Cartan- Killing metric on g. According to standard estimates ||31 
on the heat kernel on a compact Riemannian manifold, there is a neighborhood U of the 
diagonal in Hom{7Ti{M),G) x Hom{7Ti{M),G) such that for any pair of representations 
{61,62) £ U and Vol[Hom{TTi{M),G)/G] near 0, we can write 

KiNnVvolia''), 61,62) = 

{27rN^Vvol{a-))-^/'exp i^- J^'J^^l^^ ) (1 + 0{N^Vvol{a-))) (147) 

where |^i — ^2! denotes the riemannian distance, (in the induced Cartan- Killing norm), on 
Hom{TTi{M),G) C G™. 

Set t=NnVvol{(T'^), and let us consider ||32|, ||26| the representations 6, for t near 0, as 
one-parameter family of representations 6^ 

6' = exp[tu{a) + 0{t'^)]6{a) (148) 
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where a G 7ri(M), and where u: tti{M) g. Given a and h in 7ri(M), if we differentiate the 
homomorphism condition 6^{ab) = 9^{a)9^{h), we get that the u are one-cocycles of 7ri(M) 
with coefficients in the 7ri(M)-module qq, 

u{ah) = u{a) + [Ad{9{a))]u{b) (149) 

namely, u G Z^{ni{M), q), the tangent space to Hom{TXi{M),G). Among the vectors u G 
Z^{t[i{M),q), those tangent to the Arf-orbit are of the form u{a) = [Ad{6{a))]h — h, for 
some G g, these are 1- coboundaries on tti{M) with coefficients in the representation, i.e., 
elements of B\ni{M),s)- We can decompose m in a component Phu ^ H^{tti{M), g) = 
Z^{'7ii{M), q)/ B^{'Ki{M), q) and in a component PgU G g, Ph and Pg denoting the respective 
projection operators. Thus 

e' - Id = t{PHU + Pqu) + 0{t^) (150) 

Introducing this in (|147|) we get 

K{NnVvol{a''),e,Id) = 



(2vriV„V.o/(."))--/--9/. exp j (i + o(t)) (151) 

where Dh is the dimension of the smooth component of ^°'^(^^^)'^) g^Yid Dq = diniG. 
Integrating ( |151| ) over g x g^MElniE^Ml:^'^ and evaluating the gaussian integral over g we get 
to leading order: 

(2.N„V.oK.''))-^~'^-V llMP, . **P-l-'^exp (-^J^) (162) 
which, upon setting 

^(V) = llietP, . ietP-r^'^ exp (-^^ffg;^) (163) 



yields the stated result. □ 



Note that the computation of the integral ^(V), extended over the moduli space 
HomiiTi{M),G) ^ related to the theory of Reidemeister torsion , ||3^ . 

Moduli asymptotics for DT-surfaces. Both the presence of the volume of the represen- 
tation variety and of its dimension, are an important aspect of ( |146| ) . In dimension n = 2, 



we can be more explicit, and if we gather information on the dimension of -^°'"(""^(^)''^) pg 
33| , then in the case of 2-dimensional surfaces theorem |Il] specializes to the following 



Lemma 5 For n = 2 and for VfoZ(cr") near zero, the leading contribution to Vol[Ai{Ta)] 
is given as a function of surface topology by: 



(i) The sphere S^; 



Vol[M{Ta)] -^v.oKa")«i^ 1 (154) 



(a) The torus T^; 

Vol[M{Ta)] -^vvoKa^)«i^ A{V) ■ (27riV„Vt;o/(a"))-^/2 (155) 

(Hi) Orientable surfaces S/j of genus h > 1, and for representations o/ 7ri(S/j) —>■ G 
whose center has a constant dimension dimZ : 

Vol[M{Ta)] ^vvoK.-)«i^ A{V) ■ (27riV„Vt;o/(a-))-[(2^-^)'='^-^+'='^-^l/2 (156) 



(If we have a non-trivial centralizer for the representation, then we have to slightly alter the 
argument in theorem |Tl|) . 
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Proof. For a sphere dim[H^ {7ri{M) , g)] = 0, and trivially V = 0. As V — » 0, the measure 
weakly converges to 6{Id), where 6{Id) is the Dirac distribution concentrated on the trivial 
representation. This immediately yields the stated result. 

The torus T^, can be actually considered as a limiting case of (iii). From it follows 
that if z{9) denotes the centralizer of ^(7ri(S^)) in G, then the dimension of Hom{ni{T,h), G) 
is given by 

{2h - l)dim{G) + dim{z{e)) (157) 

By taking G = SL{2, m), and setting h = 1 we immediately get the stated result for T^. The 
general result for /i > 1 is a trivial consequence of the dimensionality formula ( [157| ). □. 

We wish to stress that these results (in particular the asymptotics of the volume of moduli 
space) refer only to the smooth component of the representation space Hom{7Ti{M), G)/G. 
Singularities are present in correspondence with reducible representations and around rep- 
resentations which are linearization unstable [^. Also, Hom{TTi{M)^ G)/G may have many 



distinct connected components. For example, in the relevant case of G = PSL{2,R), 
Hom{-Ki{M), PSL{2, m)) has 2'^^+^ + 2h-3 connected components, p6[. In 2d- Yang- Mills, it 



is a known fact that, at least in the semi classical limit, the singular component of the repre- 
sentation variety does not contribute to the volume of the moduli space , [^] , . Since 



there are many points of contact between our approach in evaluating the asymptotics of the 
volume of the moduli space of Romer-Zahringer deformations of a dynamical triangulation, 
and 2d- Yang-Mills, we may reasonably assume that the singular part of the representation 
variety is not playing any relevant role in our case, at least in dimension n = 2. For dimen- 
sion n > 2 this is an open issue under current investigation. As we shall see in due course, 
the asymptotic of the volume of the moduli space of Romer-Zahringer deformations plays a 
basic role in simplicial quantum gravity being related to the entropy exponent of the theory, 
thus these matters are quite relevant in addressing on analytical grounds the existence of 
the continuum limit of the theory. 



83 



A compact formula for surfaces. We conclude this section by providing a compact 
expression for the the asymptotics of the volume of the moduli space of Romer-Zahringer 
deformations of a dynamical triangulated surface. In the two-dimensional case, according 
to lemma ^ the small volume behavior of Vol[Ai{Ta)] is 

A(V) ■ iV«(^'^) (158) 

where the constant ift(V)=A(V)(27rVi;o/(a"))-[(2'*-i)'^*'^^+'^*'"^]/^ (for genus h>l, whereas 
^h(V)=A(V) for the sphere), and ^(S/^) is an exponent depending from the genus, and the 
group G. This exponent can be given a more compact expression if we introduce an effective 
dimension, De/f, of Hom{7ri{M),G) according to 

Deff=i2h - l)adim{G) + adim{z{9)) + (3 (159) 



where the parameters a and (3 are determined by fitting ( |159|) with the dimensions of the 
singular case provided by the sphere and the limiting case of the torus T^. Explicitly, 
according to lemma ^ we have 



e(S,) = -{2h - l)'^dim{G) - ^dim{z{e)) - | + 1 



(160) 



with 



^dim{G) - ^dtm{z{9)) - f + 1 = ^(s') 
.^dtmiG) - ^dim{z{e)) - f + 1 = e(T^) 



1 - 



;i6i) 



In this way we get 



a 



/3 



2dimG 



dim{z{9)) , 
diniG ' 



(162) 



and 



(163) 
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where x(S/i)=2 — 2h is the Euler-Poincare characteristic of the surface E/j. These observa- 
tions estabhsh the following 

Lemma 6 The small volume limit asymptotics of the volume of the moduli space of Romer- 
Zdhringer deformations of a dynamical triangulated surface is given by 

Vol[M{T,)] ^v,„;(,„^«,^ ^(V) • ivl^^'^''^"^ (164) 



85 



VII. CURVATURE ASSIGNEMENTS FOR DYNAMICAL TRIANGULATIONS 



According to the analysis of the previous sections, the topological information associated 
with a dynamical triangulation is encoded in the moduli space of its Romer-Zahringer de- 
formations. Here, our purpose is to characterize the metric properties of the triangulations 
in terms of the (restricted) holonomy of its natural Levi-Civita connections, and reduce the 
enumeration of distinct dynamical triangulations to the enumeration of distinct curvature 
assignements. 

Simplicial Lassos. We start by introducing the following 

Definition 11 A simplicial loop v, based at a simplex ctq, is a (small) lasso with nose at 
the hone B, if it can be decomposed into three simplicial curves ^{B) = (B) ■ T, where 

T is a simplicial path from ctq to one of the simplices a^{B) containing the bone B, r^^ is 
the same curve going backward, and uj{B) is the unique loop winding around the bone B. 

Since the open stars of the bones provide an open covering of M = \Ta\, any simplicial 
loop in M, if it is homotopic to zero, is a product of lassos whose nose is always contained 
in the open star of some bone B, (this is the Lasso lemma I^S])- Thus, in order to evaluate 
the holonomy around a generic contractible simplicial loop based at a marked simplex 
we can proceed as follows. 

Let us fix a holonomy representation of the fundamental group of M, 6: 7Ti{M) ^ G — *• 
End{s). Recall that the choice of such 9 calls for the marking of a bone o"o~^, (whose 
barycenter provides the basepoint po), of a simplex o"g D ctq"^, and of the holonomy around 
the marked bone. As the notation suggests, we naturally identify the marked simplex cTq 
with the simplex on which the simplicial lassos are based, also, for notational simplicity, we 
let A + 1 = Nn-2{T) denote the number of bones in the dynamical triangulation T^. Let us 
consider the set of simplicial lassos, {z^(a;), ctq } based at the marked Ug and with their noses 
corresponding to the A + 1 bones Bo,B{l), . . . ,B{X). The effect of parallel transporting, 
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according to the Levi- Civita connection associated with T^, along any of these lassos is a 
rotation in a two- dimensional plane, (orthogonal to the bone considered). If uj is the generic 
(contractible) simplicial loop based at ctq , then the holonomy along the generic loop uj based 
at (Jo, can be written as 

A 

^-(^^o) = X{{A{a)-^R{m{a)ct){a))A{a)} (165) 

where m{a) G Z are the integers providing the winding numbers of uo around the various 
bones, R[m{a)(f){ay\ denote the rotations in the planes orthogonal to the bones -B(a), and 
finally, A(0), A{a) denote the orthogonal matrices describing the parallel transport along 
the ropes, Tq, of the lassos, (^4(0) is associated with the trivial path ctq)- 

On rather general grounds, it is known ||36| that through holonomy we can reconstruct 
the underlying connection, (up to conjugation). In our case, as already stressed, this recon- 
struction procedure is particularly relevant since distinct connections are in correspondence 
with distinct triangulations. It follows that the conjugacy class of the (restricted) holonomy 
Ruj{o'o) is G-invariant and fully describes the underlying dynamical triangulation in the given 
representation Q:tti{M,Pq) —>■ G. The natural class functions of the holonomy (and hence 
of the triangulation) are linear combination of the group characters, namely of the traces 
of -Ra;(o"o), (these latters being thought of as elements of the fundamental representation of 
the restricted holonomy group). These traces depend only from the given set of incidence 
numbers {g(a)} of the triangulation, and a question of relevance is to what extent such 
traces, and hence the curvature assignements, characterize the triangulation itself. In order 
to answer such a question we exploit the connection between curvature assignements an the 
theory of partitions of integers. 

A. Partitions of integers and curvature assignements 

According to (^2|) the curvature assignements {g(a)} of a dynamical triangulation Ta 
cannot be arbitrary, since they are constrained by the total number Nn{T^'^^) of ra-simplices 
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in Ta- We have 

A 1 

q{o) + q{a) = -n{n + l)iV„(T«) (166) 

a=l ^ 

where the summation extends over the bones of Tj*\ (notice that we have included in the 
summation the marked bone, by considering also the incidence number q{0)). Another 
natural constraint comes about from the fact that in order to generate an n-dimensional 
polyhedron there is a minimum number q{n) of simplices that must join together at a 
bone, (according to the introductory remarks, in dimension n = 2, n = 3, n = 4 we have 
q{n) = 3, otherwise we have no polyhedral manifold. Rather than use this explicit value, we 
keep in using the general notation q{n)). If we define 

Q(0)=g(0) - q{n) 

Q{a)=q{a) - q{n) (167) 
then we can put these two constraints together by writing 

q(o)+x:q(«) = 

(recall that Ar„_2(T») = A + 1, also recall that 6(n,n-2) = in(n + l)(A/'„(T(*))/A/'„_2(rW)) 
is the average number of simplices a", in T^^\ incident on a bone). 

Thus, we have 

A 

Q(0) + J2 Qi») = Hn, n-2)- qin)] (A + 1) (169) 

a=l 

Recall that a partition of a positive integer n is a representation of n as the sum n = 
1 ■ rji + 2 ■ r]2 + . . . + n ■ rjn, where r]i,r]2, . . . ,i]n are non-negative integers, and i]i+r]2 + . . . + r]n 
determines the number of summands of the partition. 

This remark suggests that for each given value of Q{0), b{n,n — 2), and A, the set of 
curvature assignements q{a) = Q{a) + q satisfying the constraint (|169|) may be related 
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-n(n + 1) 




q{n) 



(A + 1) 



(168) 



to (some of) the possible partitions of the integer [b{n,n — 2) — q{n)] (A + 1). Among the 
attractive features of such an interpretation is the possibihty of exploiting the large body 
of knowledge we have on partition theory, a well-established chapter of combinatorics and 
analytic number theory. In order to exploit these techniques, we need to modify a little 
( |169| ) so that it resembles more directly a partition. 

First of all note that it is natural to assume that the marked q{o), (or equivalently (5(0)), 
is the largest curvature assignement in the dynamical triangulation considered Tq, i.e., q{0) > 
q{a), ( thus Q(0) is a sort of root blob of curvature). As we have seen in the introductory 
chapter, (see (|65D), the largest curvature assignements can well be of the order of A, and 
the largest curvature assignements which simultaneously minimize the remaining curvature 
assignements on the other bones are realized on particular triangulations, (see the discussion 
following (|65|)). In such triangulations, most of the top-dimensional simplices are incident 



on n distinguished bones, (bounding an embedded §"~^). Notice that, when generated with 
equilateral simplices, these particular triangulations have a discrete simmetry associated 
with the interchange of the n distinguished bones carrying most of the curvature. Since we 
have to consider ( |169| ) as g(0) varies from q to qmax, we have to factor out these residual 
symmetries, and we do so by marking not just one bone, (the one providing the basepoint 
Po), but rather n bones, (t"'~^(0), (j"~^(1), . . . , (T"~^(n— 1), which support, in the extreme case 
described above, the largest curvatures (with g(0) = g(l) = . . . = q{n — 1)) while minimising 
the incidence on the remaining bones {o""~^(a)}Q,=n,...,A- The original marked bone ctq"^ 
is one of those n distinguished bones, and can be selected by the required marking of an 
ra-simplex o"q associated with the gauge fixing procedure for the moduli space of dynamical 
triangulations. Thus, we start by replacing (|169|) with 



w 



Q(0) + = IKn, n-2)- q{n)] (A + 1) (170) 



Since nQ{{)) > Q{a), we can refer the curvature assignements Q{a) to the chosen value 
for the curvature assignement nQ{0) around the marked bones. To this end, let us add and 
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subtract to ( |170| ) the expression nQ{0){\ + 1), which allows us to eventually rewrite (|169|) 
as 

A 

<3(«) = [nQ{0) - b{n, n - 2) + g] (A + 1) - n{n - 1)Q{0) (171) 

a=n 

where 

Qia)=nQiO) - Qia) (172) 

This form ( |171D of ( |169| ) can be interpreted as providing a partition of the integer [nQ{0) — 
b + q]{\ + 1) — n{n — 1)Q{Q) in at most A + 1 — n integers Q(a), (some of the Q{a) can be 
zero), each bounded above by {h — q){\+ 1). 

Notice that is not a priori true that each partition of the integer [nQ(0) — 6 + g] (A + 1) — 
n{n— 1)(5(0) can be generated by curvature assignements through the Q{a). Thus, we shall 
denote by 

(hQ(0) - 6 + g](A + 1) - n{n - l)g(O)) (173) 

the number of partitions of [nQ{Q) — b{n,n — 2) + g](A + 1) — n{n — 1)Q{0) into at most 
X + l — n parts, (each < (6(n, n — 2) — g)(A + 1)), arising from actual curvature assignements. 

A notational remark. Notice that in comparing with standard formulae for partition 
theory , a more standard notation would have been 



pZi-u {[nQ{0) - 6 + g](A + 1) - n(n - l)Q(O)) (174) 

however, in order not to burden the notation further, we shall adopt the above simplified 
version; in some cases, if there is no danger of confusion, we shall simply write p'^^^^ for the 
above expression and for other similar expressions involving partitions. 
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B. Curvature assignements of distinct dynamical triangulations 



For small A's the actual distribution of partitions corresponding to curvature assigne- 
ments is rather accidental, however for A >> 1 we have that a large fraction of partitions 
corresponds to actual curvature assignements. We start by proving that generically distinct 
curvature assignements correspond to distinct dynamical triangulations. We can prove this 
by the elementary probabilistic methods which are typical in discrete mathematics when- 
ever attention is on asymptotic properties. We also bound the size of the non-generic set of 
distict dynamical triangulations associated to the same set of curvature assignements. The 
tails associated both to partitions of integers which are not curvature assignements and to 
distinct dynamical triangulations having the same curvature assignements provide correc- 
tions to the naive counting through unrestricted partitions of integers. Such corrections are 
difficult to control in details, however they have the structure of a rescahng that will allow 
us to implement a very straightforward counting strategy which will prove very effective. 

Let us start by noticing that it is obviously false that every possible set of distinct 
curvature assignement corresponds to triangulations having distinct incidence matrices. One 
may easily construct particular examples violating such statement. Consider for instance 
two dynamically triangulated (exagonal) flat tori. By flipping links, we may generate two 
curvature bumps on each torus, (on each torus, the bumps may correspond to distinct 
curvature assignements) . By inserting in each torus a copy of these curvature bumps in 
such a way that their distance is different in the two tori, we get distinct triangulations with 
the same curvature assignements. As another more general example, (again in dimension 2), 
consider a two-dimensional triangulation with a large number of vertices, let abc, acd, and 
efg, egh four triangles pairwise sharing a common edge (ac for the former pair, and eg for the 
latter), but otherwise largely separated, (so that their incidence numbers on the respective 
vertices are uncorrelated) . Let us assume that the corresponding incidence number arc given 
by: q{a) — a, q{h) — 7, g(c) — j3, q{d) — S, for the first pair of triangles, while for the second 
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we set: g(e) = 7 + 1, q{f) = a — 1, q{g) = 5 + 1, q{h) = [3 — \. It is immediate to check that 
a flip move will interchange the curvature assignements of the first pair of adjacent triangles 
with those of the second pair, thus changing the incidence relations in the triangulation. 
Nonetheless, the sequence of curvature assignements remains unchanged, and again we have 
two distinct triangulations with the same sequence of curvature assignements. 

Clearly, these counterexamples work either because we carefully chose a particularly 
symmetric triangulation or because we adjusted the curvature assignements to the particular 
action of the flip move. However, when moving from a dynamical triangulation to another, 
with a set of ergodic moves (Pachner ( |^)), the curvature assignements are not fixed, and 
in this sense the above counterexamples are not generic, and if we denote by {T^^^jcuru the 
set of distinct triangulations sharing a common set of curvature assignements, we can prove 
that as A >> 1, {Tj^^^^v is in a suitable sense a small subset in {Tj*^}. We can formalize 
this remark according to the 



Lemma 7 For a PL-manifold M, let 0:7ri(Af;po) ^ G, be a given holonomy representa- 
tion, and let {Tj^'>}\ denote the set of distinct dynamical triangulations of M = \Tj^^^\ with 
A^„_2(T^*^) = A + 1 bones, and let {rj*''}AeZ associated inductive limit space generated 

as A grows. Let Vt = ({^i*^}AeZ) ^t) be the probability space endowed with the normalized 
counting measure nix on every non-empty open set B C {Tj~^^}xelj- '^'^V ^ {^a*''}Aj 
{g(a)} describe the sequence of integers providing the incidence of the simplices {cr"}j^\ " 
over the A + 1 bones {a1^'^}^l'{^'^'^'^\ Then, as X 00, the curvature assignements {q{a)} 
characterize with probability one the dynamical triangulation Ta[q{a)]. 



Proof. We start by creating a probability space Vt whose elements are the distinct triangu- 
lations in {T'j*)};^^^, and whose Borel probability measure, rrix, is the normalized counting 
measure on every non-empty open set B C {7"i*^}AeZ- Note that we have to consider trian- 



gulation with variable A, since the known set of ergodic moves |^ mapping a triangulation 
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T^(*)(1)a into a distinct triangulation T^')(2)a, (with the same A), go through intermediate 
steps which do not preserve volume (hence the number of bones A at a given b{n,n — 2)). 
If, for a given value of A, {Tj^'^^curv is the subset of distinct triangulations with a same set 
of curvature assignements, we can write (A fixed) 

"r[{ri >}„„]= (175) 



Recall |^8|] that on {VT,B,mT), the Pachner moves characterize a continuous and ergodic 



dynamical system with respect to m^. If 



lim > o (i76) 



then by Birkhoff ergodic theorem we would get 



[{T!,%ur.] > (177) 



and the Pachner moves would fix mr-almost everywhere the set {T^^^}curv thus contradicting 
their ergodicity. Hence we must have 



, Carc/{T«} 
hm ' 

Card{Tr} 



lim rz::ii2iL^i:!fiii = (178) 



which reflects the fact that Pachner moves can fix the curvature assignements only under 
special circumstances. □ 

Even if this formal result may at first sight appear very reassuring we still have to 
face the fact that distinct dynamical triangulations associated to a same set of curvature 
assignements may give rise to significant subleading corrections to the counting of distinct 
triangulations of a given PL-manifold. Thus we need a way for controlling how large the 
set {T^^^curv can be. This control is provided by the following result 



Lemma 8 There are constants, n{n) > 0, and T[n) > 0, depending only on the dimension 
n, such that 
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(179) 



Proof. For a given set of curvature assignements over the A + 1 bones of a dynamically tri- 
angulated PL-manifold M, let {Tj^^^ }curv denote the set of distinct dynamical triangulations 
(with A + 1 bones), having the same curvature assignements. Consider the corresponding 
dual complexes, viz., for a given G {Tj^'^^curv, consider its barycentric subdivision, then 
the closed stars (in the subdivision) of the vertices of Ta form the collection of n-cells of the 
complex dual to the given T^. Distinct triangulations in {Tj^^^}curv correspond to distinct 
dual complexes. The use of the dual complex rather than of the triangulation itself is related 
to the fact that the incidence matrix of the dual complex is directly related to the curvature 
assignements of the original triangulation. This is most easily seen in dimension n = 2, 
whereas in dimension n > 2 the curvature assignements of a dynamical triangulation Ta, 
over the respective bones, can be read from the incidence matrix A = {a'ij}o<i,j<x associated 
with the 2- skeleton of the dual complex, {e.g., in dimension n = 3 think of a froth pattern 
in space). Explicitly, one has that the line sums of the 2- skeleton incidence matrix A are 
such that 



where q{k) is the number of simplices in {T^^^curv which are incident on the bone a"'~^{k). A 
similar relation holds for the column sums. In order to estimate Card{Tj^^'^}curv for a given 
set of curvature assignements {q{k)}o<k<x we have to enumerate the set of such distinct 
incidence matrices. The generic matrix A of such a kind, is a symmetric (0, 1) matrix 



constraints, associated with the Dehn-Sommerville relations, have to be satisfied). In order 
to enumerate such matrices we can proceed as follows. We have 



A 



Y^akj = q{k) 

j=0 



(180) 



with vanishing diagonal entries satisfying, for each line, at least the relation (|180|) , (other 



A- 1 



(181) 
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choices for the first row {aij} of A, and in general, since A is a symmetric traceless matrix, 
once chosen the k — 1-th row, we have 



X-k 



182) 



choices for the fc-th row. 



Card{Tjf>}curv is bounded above by the number of distinct (0, 1) matrices satisfying the 
constraint ( |18CI|) , thus we get 



fc=0 y 



) 



(183) 



By taking logarithms and on applying Jensen's inequality, > J^k^^^k), we get 



/ 



ln[Cart/{T«}e,,J <ln^ 

k=0 



X-k 



q{k)-{T.Uakj) 



(184) 



We set 



k-l 



T{n) = max[q{k) - ^ a^j] 

k . 



and bound the binomial coefficient in (|184| ) according to 



X-k 

^k-l 



< 



{X-k)\ 



qik) - (E;=! ak,) ) 



{X-k-T{n))\ 



185) 



(186) 



On applying Stirling's formula A! = (A/e)^y27rAe"^ , where 1/(12A + 1) < ax < 1/(12A), 
we easily get from ( |184| ) that for A >> 1 



Card{T^^'']curv < ^-in) ■ X 



T{n) 



187) 



as stated. □ 
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The above bound on CardlTj^^^curv is quite rough but sufficient for our purposes, since 
it imphes that Card{Tj^^^(.urv can, at worst, give rise to a subleading polynomial correction 
to the counting associated with the curvature assignements. 

This remark implies that the knowledge of the asymptotics of the number of partitions 
pcurv q£ ^YiQ integer [nQ{0) — b + q]{X + 1) — n{n — 1)Q{0) which are curvature assignements 
provides a way of counting distinct triangulations of a given PL-manifold, at least when A 
is large. 

Trivially p™'"'' < p^^ and a ffist estimate of the asymptotic behaviour of p™'"'' can be 
obtained by the asymptotics of px: 

Lemma 9 Let VTn{a,b, X) denote the set of distinct dynamical triangulations, of an n- 
dimensional PL manifold M, with Nn-2{Ta) = A + 1 hones, and average incidence b{n,n — 
2) = b. Then, for large A the number of distinct partitions of the integer [nQ{0) — b + q]{X + 
1) — n{n — 1)Q{0) (into A + 1 — n parts) coming from curvature assignements of a dynamical 
triangulation Ta G VTn{a,b, X) is exponentially bounded above according to 



Proof. The asymptotic behavior of Px{[nQ{Q) — 6 + g](A + 1) — n(n — l)Q{Qi)) follows by 



noticing that (see e.g pp. 32) p\{[nQ{<S) — 6 + g](A + 1) — n{n — 1)(5(0)) depends only 
on the value of [nQifS) — 6 + g](A + 1) — n{n — l)Q(O) modulo (A + 1 — n)!, and that it is a 
polynomial of degree X — n whose leading term is 

{[nQ(0) - 6 + g](A + 1) - n(n - l)Q(O)}^-" 



{X + l-n)\{X-n) 
On applying Stirhng's formula, we get for A >> n 



:i89) 



_ , H- ,1(A H- 1) - - 1X2(0)) . g ,.|„'^g'"_'-„'';!'y,,i, (190) 
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As Q{0) varies, p\{[nQ{0) —b + q]{X + l) —n{n — l)Q{0)) is bounded above by the right side of 
the above expression evaluated for max[nQ{0)] = (b — q)\, which yields px ~ e'^'^^'^^^^^^^'^A^^. 
Since p'<^^^ < px, the result follows. □ 

The bound (|188|) may be consistent as well with a subexponential growth for p^"™, a 
subexponential growth that, if present, may make the actual control of p'^x^^ quite unassail- 
able. However, more can be said on the actual asymptotics of the curvature assignements in 
T>Tn{a, b, A). Since some of the Q{a) may be zero, the number of distinct partitions of the 
integer [nQ{0) — 6 + g](A + 1) — n{n — 1)Q{0) into A + 1 — n parts coming from curvature 
assignements is approximated with greater and greater accuracy, for A large, by the unre- 
stricted partition of [nQ{0) — 6 + g](A + 1) — n{n — 1)Q{0) into curvature assignements, {viz, 
we no longer require to have a partition into X + 1 — n parts) . This latter partition function 
is in general much more difficult to estimate than the original one, however its asymptotics 
can be qualitatively controlled through known techniques of analytic number theory, and 
together with the elementary bound ( |188| ) we can readily establish an exact exponential 
asymptotics for p'j^^^. 

For a given value of b, and A >> 1, let C{b,X,Ta), (C for short), be the set of positive 
integers occurring as curvature assignements of some rooted triangulation of bounded 
geometry in X>T„(a, b, A). Set k = [nQ{0) — 6 + g](A + 1) — n{n — l)Q(O), and for t denoting 
a generic indeterminate, let 

/cn™(t) = EP"™(^)^' (191) 

fc>0 

be the generating function for the set of unrestricted partitions, p'^^^'"{k), of the integer k 
into curvature assignements of some Ta G T>Tn{a, b, A). Then, by exploiting the formula for 
the sum of a finite geometric series, one can easily show that for |t| < 1 

/cn™(t) = n(l-^')^' (192) 

kec 

where A; G C is a shorthand notation for denoting a partition of A; = {ki, k2, ■ ■ .) whose parts 
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belong to C, (recall that for large A distinct partitions generically correspond to distinct 
triangulations) , and where 

J] (1 - t'')-^ = (1 - t'^i)-! X (1 - X (1 - &)-^ X . . . (193) 
keC 

If we denote the partition /c G C in the form 

k = (r^2''^3'''' . . .) (194) 

where exactly a^of the {kj} are equal to i, we can eventually rewrite fcurvif) as the uncon- 
strained product 

oo 

/cn™(t) = n(l-^')~"^- (195) 



At this stage, we can exploit a general result due to G. Meinardus providing the 
structural asymptotics of partition functions which are associated with infinite product gen- 
erating functions. Such a technique is neatly described in Andrews' monography, and we 
shall follow closely his exposition, (see Chap. VI). 

Let us consider the Dirichlet series associated with the set of multiplicities {aj} associated 
with a curvature assignement partition /c G C, viz., 

D{s) = f:^ (196) 
j=iJ 



with s = u + iv, and let a G be its abscissa of convergence, {i.e., ( |196D is assumed to 



converge for s > a). Let us further assume that D{s) admits an analytic continuation in the 
strip u > —Co with < Co < 1, and that in this region D{s) is analytic with the exception 
of a pole of order one located at s = a with residue Res[D{s)]s=a = Res{a). Under a 
number of further technical assumptions which are not relevant for our purposes, we get the 
following result which is a straightforward application of Meinardus' theorem: 

Lemma 10 As k ^ oo, the leading asymptotics of p'^^'^'" {k) is given by 
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pcurv^f^^ ~ C/t^exp + ^)[Res{a)r{a + l)C{a + (1 + . . .) (197) 



where C is a suitable constant, (see e.g., Theorem 6.2 in W\] ), (1 + • • •) stands for higher 
order terms, r(a+l), C{a+1) respectively denote the Euler gamma function and the Riemann 
zeta function, and finally 

D(0) - 1 - , , 

1 + a 



In our case, k stands for [nQ{0) — b + q]{X + 1) — n{n — 1)Q{0) and for nQ{0) = {b — q)X, 
we get from ( |197D 

pcurv ^ ^[(^ _ g)A2]'7exp |[(6 - g)A2]"/(°+i)(l + ^)[Res{a)r{a + l)C{a + l)]V("+i)| (1 + . . .) (199) 



This asymptotics is consistent with the bound ( [L88| ) only if the abscissa of convergence of 
the series -D(s) is a = 1, and if its residue Res{a) and the value -D(O) are bounded according 
to 



/ 6 2 + ln(6-g) , , 



D(0) < -\ (201) 



C. The counting principle 



The bound (|188|) and the asymptotics of lemma [T0| establish that p^"^'' has an exponential 



growth with A, and that the coefficients a„ characterizing the Dirichlet series D{s), (and 
the partitions in curvature assignements), have to satisfy some growth condition (related to 
the above bounds). Clearly, more geometric information would be necessary for getting a 
complete control of p^"™, however what really matters for our purposes is the observation 
stemming from the above analysis that the ratio 
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^ = c„ (202) 
Px 

where < c„ < 1, either is a constant or an exponentially decreasing function of A, and 
b{n,n — 2) viz., Cn = exp[(ai — a2b)X + ao] for some suitable constants ao, ai, a2- This 
functional dependence of c„ on A can be conveniently recast in the form 

^ = c„ = Coe°"^"(^») (203) 
Px 

where the constants < Co < 1, and a„ < depend only on the dimension n. 

Since c„ has exactly the structure of the exponential of the volume part of the dynamical 
triangulation action, and since p^"™ generically enumerates distinct dynamical triangula- 
tions, we can as well as use for such a purpose the much more manageable counting function 
Px provided that in the final results we allow for a renormalization of the coupling A;„, 
according to 

K^K + 6n (204) 

with 5„ < 0, a constant depending only on the dimension n. 

In other words, by using p\ as an enumerator, we are using a discrete measure on a 
space (the space of distinct partitions of the integer [nQ{0) — b + q]{X + 1) — n{n — l)Q{Qi)) 
which contains the relevant space of curvature assignements, but which is definitely larger. 
The discrete measure associated with px induces a corresponding measure, p^"^^ on the 
space of distinct dynamical triangulations, and the two measure scales according to (P03D , 



(which appears as a the discrete Radon-Nykodim derivative of the measure p^"^^ with respect 
P\). The nature of this scaling allows us to use directly the overcounting measure p\ since 
in seeking the infinite volume limit and a possible scaling limit for the partition function 
(p6|) we have to vary both couplings Roughly speaking, in using directly p\ in 

place of p^"'"^ it is more or less like including tadpoles and self-energies in the triangulation 
counting. The shift (p03|) and ( p04| ) allows us to remove, in the final results, the effect of such 
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inclusions. We formalize the results of the lemmas § ^, ^ and of the associated remarks, 
in the following lemma providing the rationale of our counting strategy, 



Lemma 11 There exists constants, < Cq < 1, a;„ < 0, and T{n) > depending only on 
the dimension n, such that the number of distinct curvature assignements of rooted dynamical 
triangulations Ta G VTn{a,b, X) is enumerated, up to the scaling 

Card{T^'^ 



Px 



CoAr^(''^e""^"('^''\ (205) 



by the set of distinct partitions px of the integer [nQ{0) — b + q]{X + 1) — n{n — 1)(5(0) into 
at most X + 1 — n parts, ( each < (b — q){\ + 1)). 

D. A geometrical remark 

According to the above results, curvature assignements generically characterize dynami- 
cal triangulations. In the two-dimensional case this can be related to the fact that curvature 
assignements provide the rotation matrices defining the holonomy, and thus the triangula- 
tion. Let us verify that geometrically also for n = 3, and n = 4 curvature assignements 
provide most of the information characterizing the holonomy for dynamical triangulations. 
We work out in details the three-dimensional case by using an SU{2) connection. The case 
for dimension n = 4 can be easily discussed along the same lines by exploiting the isomor- 
phism 5*0(4) ^ sui2)xsui2) ^ fjj^jg analysis is due to J. Lewandowsky, to whom we wish to 
express special gratitude. 

Consider the marked simplex ctq, (in this case a tetrahedron), and fix an orthonormal 
basis {e^^y^z} in ctq- Among the six bones (edges) in the boundary of o"q, we can choose 
^n{n — 1) = 3 bones, Tj, i = 0,1,2, one of which is the marked one, tq, and all sharing 



a common vertex, (recall the need of marking n bones in interpreting ( [L69| ) in terms of 
partition of integers!). These bones are defined by the vectors E{Ti) G M^, (with components 
referred to the fixed orthonormal basis {cx^y^z})- Set 
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n = + + (206) 

\E{n)\ '\E{r^\ '\E{n)\ 

with ax, (Ty, <Jz denoting the Pauh matrices. The holonomy matrix describing the 2- 
dimensional plane rotations generated by winding around the bones Tj can be written as 
the SU{2) matrix 

U{Ti) = I cos((/)(i)/2) + Ti sin(0(i)/2) (207) 

where / is the identity operator, and 0(i) = q{i) cos~^ | denotes the rotation angle around the 
bone Ti. In general, we can reconstruct the holonomy matrix Uiji) for every distinguished 
bone Ti in the marked base simplex ctq, by giving the traces of U{Ti), and the traces of 
U{Ti)U{Tj), i 7^ j. The traces tr[U{Ti)] are proportional to the incidence numbers q{i), 
i = 0,1,2, while tr[U{Ti)U{Tj)] are proportional to the scalar product in SU{2) between 
the spinors fi and fj, respectively associated to the bones r^, Tj. These latter data are 
actually already known, since the base simplex is equilateral. Thus the SU (2)-holonomies 
around the base simplex ctq are trivially determined once we give the corresponding curvature 
assignements. The reconstruction of the SU{2)- holonomies around the remaining A + 1 — 
|n(n + 1) = A — 5 bones = {a'^~^,0 < a < A, cr^"^ ^ ctq}, is similar. Formally, 
for such a reconstruction, we would need the assignement of the traces x{c() = trlU^Ta)], 
and of x(Q;,j) = tr[U{Ta)U{Tj)], with j = 0,1,2. Indeed, with the characters x(q;) — 
tr[U (to)] — 2 cos we can reconstruct the rotation angles (j){a), while with the characters 
X(«,i) = tr[U{Ta)U{Tj)], viz., 

X(a,j) = 2 {cos ^ CO. ^_ < ?„, q > ,i„ m (208) 
we can reconstruct the SU (2) inner products 

<Ta,Tj>^-^Tr[fa-fj\ (209) 

and hence obtain the spinors through which the 5'C/(2)-holonomies [/(tq,) = 
/cos(0(q;)/2) + faSm{4>{a)/2) are defined. 
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The characters x(a) = tr[U{Ta)] with a — 0, 1, . . . , A are given (up to a constant fac- 
tor) by the incidence numbers {q{a)}'^^Q, and thus are known. Since the simphces in the 
triangulations are all equilateral, the bones fj are not arbitrarily distributed with respect 
to the reference bones Ta, but are in correspondance with the 12 points on the unit sphere 
§^ marked by the vertices of the inscribed Icosahedron, (arrange the icosahedron in space 
so that the uppermost vertex coincides with (0,0,1)). Hence, the SU{2) inner products 
< Ta-iTj > are known, and the characters x(q;,j) = tr\U {Ta)U {tj)] are given in terms of 
the {q{a)}a=Q- It follows that the curvature assignements {g(a)}^^Q provide grosso modo the 
SU (2) holonomy matrices, (up to conjugation), what is missing at this detailed level, (which 
is independent from any asymptotics) , is the combinatorial data of how to put together the 
local holonomies so as to reconstruct the underlying dynamical triangulation T^. 
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VIII. ENTROPY ESTIMATES 



According to the results of the previous section, dynamical triangulations are generically 
characterized by the corresponding curvature assignements {g(a)}, and the simple counting 
of the possible curvature assignements provides a very accurate estimate of the number 
of distinct dynamical triangulations a PL manifolds (with given fundamental group) can 
support. The previous results have already evidentiated the close connection between the 
theory of partitions of integers and the curvature assignements, our strategy is to further 
exploit this connection in order to estimate the number of distinct dynamical triangulations. 



A. The generating functions for the enumeration of dynamical triangulation 

Since the partitions of [nQ{{])—h{n,n — 2) + q\{\ + l)—n{n — l)Q{{]) into at most A + 1 — n 
parts, (each < {b{n,n — 2) — g)(A + 1)), can be thought of as enumerating, in the large A 
limit, distinct curvature assignements, theorem ^ implies that they also enumerate distinct 
rooted dynamical triangulations. This observation establishes the following 

Theorem 12 Let 0:7ri(M;po) ^ G, be a given holonomy representation of an n- 
dimensional PL-manifold M , n>2. The number W{Q, X,b,Q{0)), of distinct rooted dynam- 
ical triangulations on M, with A + 1 bones and with a given average number, b=b{n,n — 2), 
of n-simplices incident on a bone, is given, in the large A limit, by 

W{e, A, 6, Q(0)) = c„ ■ Px{[nQ{0) - b + q]{X + 1) - n{n - l)g(O)) (210) 



where Cn is the rescaling factor of lemma 11 



Remark. Note that selecting the value (5(0) of the curvature at the n marked bones 
is equivalent to considering dynamical triangulations with a boundary, (the links of the 
marked bone), of variable volume, (the length being proportional to the chosen Q(0)). Thus 
we have the 
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Lemma 12 The function W{e, A, b, Q(0)) = c„ ■ px{[nQ{0) - 6 + g](A + 1) - n(n - l)Q(O)) 

provides the number of distinct dynamical triangulations with boundary dTa consisting of 
the disjoint union of n = spheres of dimension n — 1, EP"^^ , each one of volume Vol[dTa\ = 
g(0)vo/[(T"-i]. 

Proof. This trivially follows from the above theorem by noticing that by removing from 
Ta the open star of the marked bones, (on each of which q{o) simplices {o""} are incident), 
we get dynamical triangulations with a boundary dTa — U^""^ of volume Vol[dTa] = 
nq{0)vol[a''-^]. □ 



B. Gauss polynomials and dynamical triangulations 



By exploiting the properties of the partitions Px{[nQ{0) — b + q]{X + l) — n{n — l)Q{0)) we 
can actually characterize the generating function for the number of distinct rooted dynamical 
triangulations. 

Theorem 13 Let < t < 1 be a generic indeterminate, and let px{h) denote the number 
of partitions of the generic integer h into (at most) X + 1 — n parts, each < {b{n, n — 2) — 
g)(A + 1). In a given holonomy represention 0:7ri(M;crQ) G, and for a given value of 
the parameter b = b{n,n — 2), the generating function for the number of distinct rooted 
dynamical triangulations with A^„_2(Tj*^) = A + 1 bones and given number 



n(A + 2 — n) 

of n- dimensional simplices incident on the n marked bones is given by 

{b{n, n - 2) - g)(A + 1) + (A + 1) - n 



g[W{Q,X,b;t)]=Cn-Y.Pxih)t'' 

h>0 



(A + 1) -n 



(211) 



(212) 



where < c„ < 1 zs an exponential rescaling factor, (see lemma |7i| ), and 



n 
m 



it\ 



(i-r)(i-t"~i)...(i-t™+^) 

(l-t"-™)(l-t"— l)...(l-t) ~ {t)n-Ut)r 



(213) 



105 



is the Gauss polynomial in the variable t. 
Note that the Gauss polynomial — — — 



lim 



n 
m 



j-j^ is a polynomial in t of degree {n — m)m, and that 

/ \ 

(214) 



m\(n — my. 



see 



p7| for an extensive review of the properties of Gaussian polynomials). 



Proof. The proof is a straightforward consequence of the fact that the Gauss 
polynomials are the generating function of the partitions p\{h). □. 

From these observations it follows that the generating function for dynamical triangu- 
lations with (spherical) boundaries is provided, up to an exponential rescaling, (see lemma 
TID , by a Gauss polynomial of degree (6 — g)(A + 1)(A + 1 — n), and we expect that their 
properties bear relevance to simplicial quantum gravity. It is interesting to remark that 
these polynomials already play a basic role in the celebrated solution of the Hard Hexagon 
Model by R.J. Baxter |4T|]. A first important property, related to the theory of Gaussian 
polynomials, that we can exploit for understanding the geometry of dynamical triangu- 
lations is related to the unimodality of the Gaussian polynomials. Recall that a polyno- 
mial P{t) = Oq + ait + . . . + a„t" is said to be unimodal if there exists m such that 
o-o ^ c^i < '^2 < • • • < flm > > O'm+2 > • • • > ^n- A simple applications of the unimodal 
properties to our case yields the following 



Theorem 14 For all {b{n,n- 2) - g)(A + 1), \, h>0 



w{e, A, b, h) = w{e, {b -q){\ + 1), b, h) 



(215) 



where W {Q, {b — q){\ + l),b,h) = Cn ■ P{b-q)(x+i)ih) , with p(^h_^)(^x+i){h) denoting the number 
of partitions of the integer h into at most [b — q){\ + 1) parts, each < A + 1 — n. Moreover 



W{e, A, 6, h) = W{e, A, b, (6 - g)(A + 1)(A + 1 - n) - h) 



(216) 



an 



d,forO<h< (b-mx+m+i-n) ^ 
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w{e,x,b,h) -w{e,\,b,h-i) >o (217) 



Proof. The proof is obtained by the known proof of unimodahty for px{h), by 
replacing px{h) with W{Q, \,b,h) = c ■ px{h). □ 

Since W{Q, A, b, h) enumerates the distinct rooted dynamical triangulations with a 
marked curvature assignement related to the value of the parameter h by 

h = [nQ(0) - (6 - g)](A + 1) - n{n - l)g(O) (218) 

the above unimodality properties tells us that grosso modo the largest entropic contribution 
to the dynamical triangulation generating function comes from triangulations such that 

,(0) . Mi±i) pig, 

(this bound being attained on the n marked bones). 

C. Asymptotics and Entropy estimates 

If in the generating function Q[W{Q, X,b;t)] we let t ^ 1 we get the enumerator of 
the distinct dynamical triangulations, for a given representation of the fundamental group. 
Roughly speaking, for t — > 1, ^[Vr(6, A, b; t)] reduces to the sum over all possible curvature 
assignements g(0) in the set of distinct dynamical triangulations in the given holonomy 
representation. We have the following 

Theorem 15 In a given represention 0: 7ri(M ; ctq ) — > G, and for a given value of the 
parameter b = b{n, n — 2), the number W{Q, A, b) of distinct dynamical triangulations with 
A^„_2(Tj*^) = A + 1 bones, n of which are marked, is given by 



w{eA,h) = cn- 



' {b{n, n - 2) - g)(A + 1) + (A + 1) - n ^ 
^ (A + l)-n 



(220) 
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Proof. The sum \imt->iQ[W{Q, X,b;t)] involving the partions px{h) can be evaluated by 



exploiting the property ( pi4| ) of the Gaussian polynomials, viz., 
limg[W{Q,X,b-t)]=Cn- E 

^ {b{n, n - 2) - g)(A + 1) + (A + 1) - n ^ 



1^ (A + l)-n 

estabhshing the desired result. □. 



(221) 



On applying Stirling's formula, we easily get that the above result yields, for large A, the 
asymptotics 



(6- g + 1)^-^+1^^+' 



(b-q) 



b-q 



A-2 (i + 0(A-5)) (222) 



In order to obtain from ( |222| )the asymptotics of the number of distinct dynamical tri- 
angulations T^, with A + 1 bones, we have to factor out the particular rooting we exploited 
for our enumerative purposes. We have first of all to factor out the marking of the n 
bones (j"^^(0), cr"~^(l), . . . , cT"~^(n — 1), and the marking of the base simplex cTq . To fac- 
tor out the marking of the n bones, we have to divide (|222|) by A". There can be at most 
j- (fc(n,n-2)-g) j^^_l_j^^ simplices cr" sharing one of the marked bones {o""~^(z)}j=o,...,n-i5 (actually 
the actual count is [ '■''^"'"~^-'~'^-' ](A + 1) where q can be larger than g; we use q for simplicity). 
Since the marked simplex o"g can be incident on any of the n marked bones, we still have to 
divide by n. For A >> 1, we get in this way a normalization factor [6(n, n — 2) — g]A"^^. 

We have also marked an orthogonal representation Oq of 7ri(M;po) corresponding to the 
marked base simplex o"o , and as discussed in section IV, besides the above purely combina- 
torial factors, we have to divide ( p22| ) also by the volume of the moduli space, V ol[M.{Ta)\, 
associated with the possible non-trivial deformations of the class of dynamical triangulations 



considered. According to theorem |Tl| this volume is provided by 



4 ' 



Vol[MiTa)] ~ Aiy) (^47rVt;o/(a") ^^^^^^ Aj (223) 
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where ];=Vol[ """''^^'J^^'^^ ] and D=dim[Hom{7ri{M), G)]. Hence, in order to get an asymp 
totic estimate for the number of distinct dynamical triangulations, we have to divide 
by 



A(V)[b(n,n- 2) - g]A"+^ ( AnVvoUa'')—^ — - 

Y n{n + 1) 



-D/2 



A 



(224) 



The Entropy Function. These observations estabhsh the following 

Theorem 16 The number of distinct dynamical triangulations, with A + 1 bones, and with 
an average number, b = b{n,n — 2), of n-simplices incident on a bone, on an n- dimensional, 
(n <4), PL- manifold M of given fundamental group tti{M) , can be asymptotically estimated 
according to 

W{X,b) ~ 



(6 -g + 1)1-2" 



2^\ {b-qy 



(6 -g + 1)^-5+1 
{b - qf-^ 



A+l 



(%^A)^/^A-^ (225) 
^ n{n + 1) 



where < c„ < 1 zs the exponential rescaling factor of lemma 11, and where we have 
introduced the topology dependent parameter 



W^=A-\V)[AixVvol{a'')f'^ 



(226) 



□ . 



As we shall see momentarily, the estimate ( p25|) fits extremely well with the data coming from 
numerical simulations in dimension n > 3, and in dimension n = 2, (where 6(n, n — 2) = 6), 
it is in remarkable agreement with the known analytical estimates Notice in particular 
that for a two-dimensional sphere the integration over the representation variety drops, and 
7251) yields 

C2 



24V67r 



44 
33 



_ 7 
V 2 



(227) 
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where we have denoted hj v = X + 1 the number of vertices of the triangulation. This 
estimate should be compared with the known resuh provided long ago by Tutte M, viz., 



1 



64V67r 



33 



_ 7 
V 2 



(228) 



Thus, if the exponential rescaling < C2 < 1, characterized by lemma [TT|, just reduces 
to a constant rescaling (c2 — 24/64), the agreement between the two asymptotic analysis 
would be excellent. 



D. The 2-dimensional case 



If we exploit the compact formula ( 164 ) of lemma ^, providing the asymptotics of the 



volume of Romer-Zahringer deformations for surfaces, then the agreement between our anal- 
ysis and the known results in dimension n = 2 is complete. For a surface T^h of genus h, 
and for representations of '7ri(S/j) in a general semi-simple Lie group G, the properties of the 

representation variety |2^ exploited in lemma ^ imply that ( p25| ) takes the form 

-44- 



A+l 



24V67r 



^-i-|x(s.)(^^ (229) 



where xi^h) is the Euler characteristic of the surface S/^, and . . . stands for terms of higher 
order. 



If we rewrite ( p29| ) as 



■ (^)"A^-=^ (1 + . . .) (230) 



where 7,; is the critical exponent, then we get 

7s(M) = 2 - ^x(M) (231) 

The exponent jsiM) defines the universal critical exponent of pure 2D-gravity. The ex- 
pression ( f^31| ) coincides with the well-known expression obtained by the use of Liouville 
theory. 
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E. The n > 3-dimensional case 



In order to check our analytical results against the MonteCarlo simulations which up to 
now are the only available results in dimension n > 3, let us rewrite the partition function 
for dynamical triangulations (^) explicitly in terms of the entropy function W{X,b), as 

Z[K.2, K] = E m\ N = ^^^^^^h e'-A e-'^"^" (232) 

where a sum over all possible A appears. This sum is needed since triangulations with 
distinct values of A contribute to the number of distinct triangulations with a given volume 
A^. Notice that in dimension n = 2 this further summation is absent since, in that case, the 
average value of 2-simplices incident on a bone (vertex) is always 6 = 6. 

The presence of this summation over A shows that in dimension greater than 2, what 
entropically characterizes the infinite volume limit of a dynamically triangulated model as 
a statistical system is not just the entropy function W{X, b) , (which can be identified with 
the entropy in the micro-canonical ensamble) but rather the effective entropy provided 
by the canonical partition function 

W{N, K-2)eff = E W{X, iV)e'="-^^ (233) 

In order to evaluate the canonical partition function and compare the results with the 
know numerical simulations we can limit ourselves to considering dynamical triangulations 
of simply connected manifolds, since numerical data are available just for the 3-sphere §^ 
and the 4-sphere S^. Already in this case computations are far from being trivial, and much 
care is needed in handling the asymptotic estimation of the sum ( p33| ) . 



According to proposition |T3|, the entropy for simply-connected PL-manifolds is 

W{X,b) 



^ (6(n, n - 2) - g)(A + 1) + (A + 1) - n ^ 



[(6-g)(A + l)«+i] 



(A + 1) -n 



(234) 
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where the factor [{b — q){X + l)"""*"^] comes from removing the labelhngs associated with the 
rootings, (see the discussion following theorem ^). For discussing the asymptotics of ( ^341 ) , 



we need Stirling's formula in the form k\ ^ \phx{k + l)'^+2e ^ ^, which is accurate also for 
small k. 

With these caveats along the way, and setting explicitly g = 3, a long but elementary 
computation provides the asymptotics 



^2vr iV^(iV - 3A + l)(iV - 2A + 1 

where N = ^n{n + 1)N, and where we have set 

/i(iV,A) = -AlogA + log^^^ V : (236) 

Thus, ( |233| ) reduces to 

W{N, k^.,)^jf = ^ . A-V^(iV-2A + l)-" ^237) 
v2vr ivy'(Ar _ 3A + l)(Ar - 2A + 1) 



where for n = 3 



and for n = 4 



2 

— = Amm < A < Amax = (238) 

6 9 



— ^min ^ A < ^max — (239) 



We can estimate the sum (|237|) by noticing that the function appearing in the exponent, 
viz., 

f{N,X)=h{N,X) + kn-2X = 

-A log A + (iV - 2A + 1) log(iV - 2A + 1) - (iV - 3A + 1) log(iV - 3A + 1) + A;„_2A (240) 

has a sharp maximum in correspondence of the solution, A = A*(fc„_2), of the equation 
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[N-3\ + If 

log = —kn-2 

^A[iV-2A + l]2 



(241) 



A straightforward computation provides 



3 " -4(fc„_2) 
where for notational convenience we have set 



(242) 



—e'-' + 1 + 
2 



27 

[ — ekn-2 + 1)2-1 

2 



1/3 



+ 



— e'"-^ + 1 
2 



27 

r_efcn-2 + 1)2 _ 1 
^ 2 ^ 



1/3 



- 1 (243) 



The structure of (|242|) suggests the change of variable 



A = (iV + l)ri 



(244) 



(we wish to thank G. Gionti for this remark), and by replacing the sum (|237|) with an 
integration, we get 



W{N, K.2)eff 



r/-V^(l-2r^)-» j^^. 



2vr A™„ ^/(l-3r]){l-2r]) 



(245) 



where 



f{ri) = -7] log 7] + {1- 2ri) log(l - 2ri) - {1 - 3r/) log(l - 37^) + kn-2V (246) 

and ?7™„ = ri{Xmin), Vmax = r]{\max)- For n = 4, we get ?7™„ ~ 1/5, r^^^x ^ 1/4, whereas 
for n = 3, r^min ~ 1/6 and r/^ax - 2/9. 



The obvious strategy is to estimate (|245|) with Laplace method, however attention must 
be paid to the possibility that, as kn-2 varies, the maximum 



3' A{kn-2)' 



(247) 
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crosses the integration limits rjjnin and rjmax- It will become clear in a moment that the 
quantities telling us when we are nearby these particular regions are 

ll^min{kn-2) = -Sgn{r}* - Vmin)\/'^[f(V*) " fiVmin)] (248) 

and 

i^max{kn-2)^sgn{r]* - r]max)\l'^[f{'n*) " f {rimax)] (249) 

This suggests that in order to control the large N behavior of the effective entropy as kn-2 
varies, we have to use uniform Laplace estimation in terms of ■?/'■ The following general 
theorem provides such uniform asymptotics. We state it first for the case involving the 
crossing of the lower integration limit r]min- A completely analogous result holds, mutatis 
mutandis, for the upper crossing r]max-i and we state it as an obvious corollary . 



Theorem 17 Let us consider the set of all simply- connected n-dimensional, (n — 3, 4j, 
dynamically triangulated manifolds. Let denote the unique solution of the equation 

^(l-777^)=ta (250) 

3 A[kn-2) 

Let < e < 1 small enough, then for all values of the inverse gravitational coupling kn-2 
such that 

C-4 - e < K-2 < C-4 + e, (251) 

the large N-behavior of the canonical partition function W{N, kn-2)e}f is given by the uni- 
form asymptotics 



^_n-l/2g(iV+l)/(r;^i„) 



27r 



Woilpminikn-2)y/^) + ^W-i{ljjminikn-2)\f^) 



N N 



(252) 



where Wr{z) = r(l — r)e^^^^Dr-i{z), (r < 1), Dr-i{z) and r(l — r) respectively denote the 
parabolic cylinder functions and the Gamma function, and where the constants ao and cti 
are given by 
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Vmin (1 "^Vn 



n-2 



(253) 



n-2 



^^-1/2 (1 _2r/*)- 



1 - 3r7*)(l - 2r/*) 



i'min{kn-2)\J -frir){ll* 



(254) 



Notice that the above expression, (in which fr=df /drj, and f^j^=(Pf/d7f), provides the 
leading asymptotics. The full asymptotics is discussed during the proof of the theorem. 
Notice also that in some circumstances, (notably in the 3-dimensional case when employing 
as variables Nq, and A^^a), there can be an equally important constribution to W{N^ kn-2)eff 
due to the upper crossing rjmax- In that case W{N, kn-2)eff is, at leading order, characterized 
by the sum of the uniform asymptotics around both f{rimin) and f {rjmax)- 

Proof. We have to provide a uniform asymptotic estimation for large of the integral 
appearing in ( |245| ) , viz., 



Am«n ^'(i _ 3^)(i _ 2r^) 

The uniformity requirement stated here refers to the existence of a (5 G M"*" such that for 
all \ip{kn-2)\ ^ ^ the error after a finite number of terms in the asymptotic expansion should 



be smaller in than the last term which is kept |^. As discussed above, the integral / 
has various distinct asymptotic regimes, according to the relative location of the maximum 
T]* with respect to the lower integration limit rjmin or to the upper limit rjmax- We explicitly 
discuss for simplicity the case where rj* may approach rjmin- The case in which rj* rjmax, 
can be dealt with similarly, and we shall indicate the necessary modifications in the proof 
and state the final result. 



In order to carry out the required asymptotic estimation, we first transform the expo- 
nential in the integral / by introducing the variable p = p{rj) according to 
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/('?) - fijlmin) = - + 'ipmin{kn~2)pj (256) 

where for rj = rjmin we assume p = and {drj / dp)\p=o > 0, so that the orientation of the 
path of integration remains unchanged. Differentiating this expression we get 

T~ = ~(P + ^min)^ (257) 

dr} drj 

In order to have dr]{p)/dp finite and non-zero everywhere we require that p — >■ —ipmin as 
rj ^ 1]*. Thus, by L'Hopital rule 



hm ^ = \ (258) 
P^-^dp V fmiv*) 

To determine the expression of the parameter ipmin{kn~2) , (see ( ^48| )), controlhng the uniform 
expansion, we evaluate /(?]) - f{r]min) at r] = r]*, {viz., for p = -ipmin)- 

fivl - fiVnun) = l^Ln (259) 

whence 



i/j^in = ±\/2[fiv*) - fiVrmn)] (260) 

The branch of the square root is selected by examining -^[fii]) " fiVmin)] at 77 = r]min, {i-e., 
for p = 0). We get 

df df d p"^ . 

— |p=0 = — |p=0 = "^("2" ~'~ ^™"^''lp=0 ~ ~Wmin (261) 

Since {drj/dp)\p=o > by hypothesis, and {df /dT])\r^^,^ > 0, for r]* > r]min, {{df /dr])]^^^^ < 0, 
for r]* < r]min), we get that 

sgnipmin = -sgn{r]* - r]min)- (262) 

Thus 

V'mm(/cn-2) = -sgn{ri* - rimin)\/'2[f{ri*) - f{r]min)] (263) 
(note that ipmin{kn-2) is a monotonically decreasing function of kn-2)- 
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With these remarks out the way, we can write 



Jo 



(264) 



where we have extended the upper hmit of integration to oo, (on introducing a Heaviside 
function), and where we have set 



G{p) = 



^-1/2(1 _2r^)-" 



df] 
dp 



(265) 



If p = —ipmin £ [0, +C)o), then we can expand the function G{p) in a neighborhood of 
p = —ipmin, (corresponding to r] = r]*). Otherwise, if p = —ipmin ^ [0, +00) we need to 



expand around p = 0, (corresponding to rj = rjmin)- According to a standard procedure |^ 
we can take care of both cases by setting 



G{p) = ao + aip + p{p + ip)Gi{p) 



(266) 



The integral / now becomes 



I = I^ + h = e(^+^)^(^— ) / («o + «ip)e-(^+i)(5P'+'^™^"P)c;. 



P 



+ e(^+^)^(^™") / pip + ^min)GAp)e-^''^^^^-^P'+^--P^dp 
Jo 



(267) 



Recall that for r < 1 



Wr[Z) 



Jo 



(268) 



where r(l — r) is the Euler Gamma function and Dr-i{z) is the parabolic cylinder function 
Thus we get 



In = e 



Wo{ll)min{kn-2)\l N) + (^mm (^n-2) V A^) 



(269) 



where the constants ao and ai are determined by 
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ao = G(0) 



(270) 



and 



Namely: 



G{0) - G{-^„ 



Vmin '^Vrr 



"^minij^n- 



^J{1- 3r7™„)(l - 2ri^in) fviVmin) 



(271) 



(272) 



n-2 



r] (1 _ 2r/)-" 



;i - 3r7*)(l - 2t]*] 



(273) 



The remaining integral /i can be reduced to the same structure of / and an inductive 



argument shows that the asymptotics of / is provided by 



Wo{i)min{kn~2)\J N + 1) 



+ 1 



+ e 



(iV+l)/(r,_„) \ W.i{i)min{kn-2)\l N + I) 



N +1 



(274) 



where the constants and 0:21+1 are determined recursively from oq and ai according to 



a2i = GiiO) 



(275) 



_ G,{0) - Ci-tPrmn) , dG, . 
"2^+1 - + ^'^=-^ 

and where the functions Gi{p) satisfy the equation 

Gi{p) = a2i+pa2i+l +p{p + 1pmin)Gi+i{p). 

dp 



(276) 



(277) 



Thus, the leading uniform asymptotics is given by 
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which provides the stated result. 



A similar analysis can be carried out when r}* approaches rjmax — \- In this case, it is 
convenient to make the preliminary change of variable y — 1 — 3r] thus reducing the integral 
I to 

" 1 ^Vmin 



If we set 



-'oo — "J 



(279) 



'^maxikn-2) = -sgn{y {?]*)) J 2[f {t]*) - 



(280) 



and introduce the variable p — p{y) according to 



f[v{y)]-f[vm 



^-J + 'ipmax{kn-2)P 



(281) 



we get again an integral whose asymptotic estimation can be reduced to the integral represen- 
tation of the parabolic cylinder function and to a set of parameters (32i and (32i+i recursively 
determined from the function H(p) defined by 



Hip)p 



-1/2 



(l_^(p))-3/2 

^y{p){l + 2y{p)) 



dyjp) 

dp 



(282) 



Explicitly we obtain 

T _ p(JV+l)/(r,„ 



^) j Wo{llJrnax{kn-2) \ N + 1) 

(TV + 1)1/2 
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_^ g(iV + l)/(»?maa.) 



(7V+1) 



(283) 



where 

6 = i^(0) (284) 

and 

5. = ' (285) 

Ymax 

The constants and ^2i+i are determined recursively from and according to 

^2^ = H,{0) (286) 

and the functions Hi{p) satisfy the equation 

Hi{p) + P— ^ = + Pe2i+1 + P{P + ^max)Hi+i{p) (288) 

Prom these expressions it is straightforward to obtain the leading 7V-asymptotics for I as 
kn-2 is such that 77* crosses 77,^03: ■ ^■ 

Summing up we have the following 

Lemma 13 Let us consider the set of all simply- connected n- dimensional, (n = 3, dy- 
namically triangulated manifolds. Let denote the unique solution of the equation 

Let < e < 1 small enough, then for all values of the inverse gravitational coupling kn-2 
such that 

k^l'^ - e < kn-2 < +00, (290) 
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the large N -behavior of the canonical partition function W{N^ kn-2)eff is given by the uni- 
form asymptotics 



27r 



JY-n-l/2g(iV+l)/(7?,„a,) 



(291) 



The notation characterizing the value of kn-2 for which we have the upper hmit crossing, 
(as compared with the notation appearing in theorem is justified by the observation, 
discussed in the next paragraph, that corresponding to this values of kn-2 we have the phase 
transitions observed in numerical simulations of simplicial quantum gravity. 

F. Distinct asymptotic regimes 

As kn-2 varies around k^nl2 and k'^%, the two parameters ipmin{kn-2) and ipmax{kn-2) 
vary monotonically from positive values to negative values. This monotonicity allows us to 
describe in details the asymptotic behavior of / in terms of kn-2- Strictly speaking we should 
consider separately both a neighborhood of k^n-2 well as a neighborhood of k^^\. However, 
corresponding to the former, available numerical simulations are so polluted by finite size 
effects that the final confrontation between the analytical results and the numerical results 
would be quite unreliable. Thus we confine ourself to a detailed analysis of what happens 
nearby kn-2^ ^^is corresponds to the value of the inverse gravitational couplings probed, 
reliably, by the MonteCarlo simulations. In any case, it is trivial to extend the conclusions 
we reach to k^f^-2- Thus, limiting ourself to a neighborhood of ^^^"2 we can distinguish the 
following cases: 

(i) Strong coupling: kn-2 < kn-2 ^ ^5 where e is implicitly characterized by the condition 
In dimension n = 4, this implies 

fc2< 1.387 -e (292) 
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(notice, that 1.387, is a theoretical bound-see the discussion in the next paragraph). Since 
fc„_2 is proportional to the inverse gravitational coupling, this asymptotic behavior describes 
a strong coupling regime for simplicial quantum gravity. For such a range of values of kn-2, 
V* < Vmax and we get 



i'max{K-2) = Sgn{T]* - r]^ax) \J '2[f (v*) - fiVmax)] = ^2[/(?7*) - /(//max)] < (293) 

In such a case the appropriate asymptotic expansion for the parabolic cylinder function is 

/ ^ (r + l)(r + 2) ^ (r + l)(r + 2)(r + 3)(r + 4) ^ ^ ^ \ ^^g^^ 



r(-r) V 2z^ 2-4^4 

(which holds for all complex z such that \z\ >> 1, \z\ » r, and j < arg{z) < ^n). Setting 



z = i^maxy N + 1 and r = —1, —2 in (|278D we get that the leading asymptotics for / is 



^ ^ r^ *-y^l~2r^*r ^ l_2n^^ ^ ^)-i/2g(^+i)/(,*) (295) 
y^(l -3r/*)(l-2r7*) V fvvW) 

(this, as expected, is twice the expression obtained in the limiting case kn-2 = k'^-2-^Ge 
below, (a), for the discussion of this limiting case). The corresponding expression for the 
effective entropy in this range of kn-2 is given by 



W{N, kn-2)eff = C, 



( {A{kn-2)+2) y\^_^_, 
H ^MK-2) ) 



■ exp 



.1 . A{kn-2) + 2 , 

[-n{n + 1) In 



{1 + 0{N~^/^)) (296) 



We have a iV " ^ subleading behavior which changes as kn-2 



(ii) Critical coupling: A;^"2 — e < kn-2 < kn-2 + *=• 

The above asymptotic is no longer appropriate when, for a given (large) value of N, ip 
is so small that 
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2i;\k'^i\±e)N « 1 (297) 
In such a regime, the parabohc cyhnder function is more correctly described by 



Dr{^) ^ 2-^/^r((l-r)/2) "^P[~V-^^~^/^" + ^^] ^^^^^ 



where 



vi ~ -sgnz — ' , ' ... (299) 

V-^ - 1/2 (2v'-r - 1/2)2 



(this holds for —r — ^ >> z^^ see pp.689). Setting r = — 1 and z = '4^rnax\l N + 1, we 



get for the canonical partition function the expression 



W{N, kn-2)eff = ^(V^„„.\/i^7I) 
V ZTT 



, ^max(l '2'^max) , 

In — 7z ^ r:^ kn-2 



max I 



1 3?7}Tta3;)(l '2l]max^ 



where . . . stands for higher order correction terms (which can be obtained by the appropriate 



substitutions from (|299|) ). 



Notice that in this latter regime, the subleading asymptotics of ( pOOj ) is no longer poly- 
nomial, but rather exponential in ipmaxVN. This transient in the subleading asymptotics is 
necessary for taking care of the fact that as kn-2 ^n-2 term (—1//^) tends to diverge. 
This divergence is compensated for by the vanishing of iprnax\J iV + 1, and indeed, the sin- 
gularity does not manifest itself in the exponential leading behavior of the effective entropy 
W{N^ kn-2)eff which goes sufficiently smootly from e^^^^^-'^'^^™^) into e'^^"'"-^^'^^''*^ but rather 
in a change of the polynomial part N~^~^l'^ of the subleading asymptotics, which changes 
by a factor iV^^^. This transient exponential subleading asymptotics is of relevance to the 
MonteCarlo simulations since it appears consistent with the numerical data for kn-2 nearby 



h,crit 
'^n-2- 
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Notice that the range of kn-2, for which the exponential subleading asymptotics is active 
depends on the actual value of A^, and as cxd, the range k^l\ — e < kn-2 < k^l\ + e 

becomes smaller and smaller since e = e(A^) is characterized by the condition 2^/'^(fc^f ± 
e)A^ ^ 1 which implies that lim7v-»oo ^{N) = 0. In this case we have 



i>max{kn-2) = Sgn{r]* - r]max)\/'2[fir]*) - fiVmin)] = (301) 



and by setting ipmax = in (|278|) we get that the leading asymptotics for / is 



y^(l-3r7*)(l-2r/*)V -^fmiv*) 
The corresponding expression for the canonical partition function is easily seen to be 



exp 



1 A{K-2) + 2 ' 
\-n[n + 1) In \N 



(1 + 0(iV-3/2)) ^3Q3) 



(Hi) Weak coupling: k'^^2 + e < < +oo. 

This is a weak coupling regime in a rather loose sense, since the asymptotics we obtain 
holds both for fc„_2's in a small right- neighborhood of k'^^2 in a large left-neighborhood 

of kn-2 = +00. 

For such range of values of /c„_2, f]* > rjmax and 

'4^max{kn-2) = Sgn{r]* - r]niax)\J'2[fiV*) " f iVniax)] = \/2[f{r]*) - f {r]max)] > (304) 

In such a case, and as long as ijj'^N >> 1, the appropriate asymptotic expansion for the 
parabolic cylinder function is 

BA.) . e-"".' (l - ^ + - . . .) (305) 

(which holds for all complex z such that \z\ » 1, \z\ » r, and \arg{z)\ < |7r). Setting 



z = i^maxy N + 1 and r = — 1, —2 in (|278|) we get that the leading asymptotics for I is 
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{N + i)-ie(*+^)^(^™-) (306) 



which provides the following expression for the canonical partition function in this regime 



277 



2 n-1 



1 r/max(l-277 

'max I 7 

m fcn-2 

(1 — 67], ''^ 



max ) 



(1 - 3r7„ax)(l - 2?7^ 



which clearly exhibits a different subleading asymptotic behavior, iV~"~^/^, with respect to 
the previous one discussed in the strong coupling case (i). 

Obviously, these results could have been obtained by standard Laplace approximation 
working in the appropriate range of values of the inverse gravitational coupling, but the 



expression (|291|) , (and its related counterpart (p52|) ), has the advantage of describing in 
a uniform way all possible regimes for fc„_2 in an open neighborhood of Roughly 



speaking, the asymptotics ( 296 ) corresponds to the location of maximum of f{r]) well'mside 
the integration interval rj^nin < rj < rj^ax, whereas the asymptotics ( |307| ) occurs when rj* is 
larger than rjmax, (but sufficiently near to it). Notice also that A{kn-2) > 1, V/c„„2 ^ ^■ 

We wish to stress the fact that the transition between the two regimes ( p96|) towards 
( p07|) manifests itself in a change of the subleading asymptotics: the growth N~^~^ which 
characterizes the expression ( P96| ) for W{N^kn_2)eff changes into the N~^~^^'^ subleading 



behavior characterizing ( |307| ), and this change is mediated by an exponential subleading 



asymptotics in the range /c^_2 ~ ^ ^ ^n-2 < ^^-2 + ^• 

The infinite volume limit. From the above asymptotics we can also read the critical 
value of the coupling /c^*"** corresponding to which the infinite-volume limit may be taken. 
In general, i.e., for all values of the inverse gravitational coupling kn-2, this critical value is 
provided by 
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(308) 



where W{N, kn-2)eff is given by ( P91|) . By considering separately the various ranges of /c„_2 
we exphcitly get 



\ A{k2)+2 Incn 
in h hm 

3 N^oo N 



(309) 



which, for < e(iV) < 1 characterized by the condition 2?/'^(A;^"2 i — 1; holds for all 
kn^2 < ~ ^- Whcrcas, for k'^^^ + ^ < kn-2 < +00 we get 



fcr(^n~2) = 2^^ + !) In 



Ijmax \i- Climax )^ ' 



+ kn-2'nmax + hm \ (310) 



Notice that if one works at finite N, (as is the case in MonteCarlo simulations), care must 
be taken in the transition region k'^^i — e < kn-2 < + ^ since in that case a non-trivial 
exponential subleading asymptotics is present, (see (ii))- 



With these elements in our hands, we can compare our analytical estimates with some 
of the known numerical data coming from MonteCarlo simulations. We start with the 4- 
dimensional case. 



126 



IX. ANALYTICAL VS. NUMERICAL DATA 



The interplay between the asymptotic estimates arising from ( 296 ) and the results com- 
ing from Montecarlo simulations requires some care, since there are manifolds which are 
algorithmically unrecognizable. Denote such a manifold by Mq, assume that n = 4, and 
let Mq be finitely described by a dynamical triangulation T{Mq). Monte Carlo simulations 
make use of two kinds of moves altering the initial triangulation T{Mq): a finite set of local 
moves [0 , and global baby universe surgery moves. The local moves are ergodic, since given 
any two distinct triangulations of the manifold it is possible in a finite number of moves, 
carried out successively, to change on triangulation into the other. The algorithmic unrec- 
ognizability of Mq means that there exists no algorithm which allows us to decide wether 
another manifold M, again finitely described by a triangulation T(M), is combinatorially 
equivalent (in the PL-sense), to Mq. It has been proved |^2| that the number of ergodic 
^ needed to connect two triangulations of Mq, T and T*, with A^4(T) = A^4(T*), 



moves 



cannot be bounded by any recursive function of N^. This result implies that there can be 
very large barriers between some classes of triangulations of Mq and that there would be 
triangulations which can never be reached in any reasonable number of steps, to the effect 
that the local moves, although ergodic, can be computationally non-ergodic. The manifold 
used in actual Monte Carlo simulations of four- dimensional simplicial gravity is the four- 
sphere which is not known to be algorithmically recognizable or not. Up to now, there 
is no numerical evidence settling this issue, even using the non-local moves associated with 



baby- universes surgery |^^, |^5|. This implies that either the class of configurations sepa- 
rated by high barriers contains only very special configurations which are of measure zero 
and not important for numerical simulations, or the barriers expected are so high that they 
have not be encountered yet in the actual simulations. It is important to have these remarks 
clearly in mind when discussing the comparision with the analytical data. 
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A. The four-dimensional case 



In order to discuss how the asymptotics associated with ( |296D compares with MonteCarlo 
simulations, let us rewrite W{N, k2)eff as 



W{N, k2)eff = lOcn ■ ^e^4"*{'=2)A^Ar^4'=2)-3(^ ^ _ _ 



(311) 



where according to theorem 17 and lemma |T^, the critical ^4"* is provided by: 



kr{kn-2) = 10 In 



A{k2) + 2 



lnc„ 

+ lim 10—^ 



(312) 



for ko < fcn"* — e, and 



A;r(^2) = 10 log 



:i - 2r]r, 



\(l~2r?,„aa;) 



77^^(1 -3?7^„,)(l-3^—) 



+ k2'nmax + lim 



In c„ 



(313) 



for /c2 > k2"^ + e and r]max = |, with 

"27 , r27 



A(A;2) = 



.k2 + 1 



gfe + 1)2-1 



1/3 



+ 



27 



e^^^ + 1 



,27 



efe + 1)2-1 



1/3 



1 (314) 



and where the factor c„ is the normalizing factor of lemma p 



The critical exponent 75(^:2) is also provided by lemma p91| , (see the discussion on the 
distinct asymptotic regimes), according to 



7^(/c2) = 2-n 



(315) 



for k2 < k^2^ — e, 



ls{k2 



-00 



(316) 



for k^2^ — e < k2 < k^2^ + e, ( recall that in this range, for large but finite A^, ( pOO| ) 
provides a non-polynomial subleading asymptotics which formally corresponds to such a 
critical exponent), and finally 
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Isih) = 2 ~ ^ 



(317) 



for k2 + e < ^2- 

Notice that these critical exponents can get an additive normahzation by the term r(n), 
(depending only on the dimension n), associated with the rescaling factor c„, (see ( |205| )). 
Thus, the actual critical exponents are provided by 

7s(fc2) = 7.(^2) + r(n). (318) 



This remark reminds us that we still have to fix the scale factor c„, affecting the canonical 
partition function W{M; N, k2)eff- According to (|204|) and lemma |ll| this factor normal- 



izes the counting of the partitions px to the partitions p^"™ coming from actual curvature 
assignements. In dimension n = 4, a simple exponential (down)shift times a constant poly- 
nomial subleading correction is all is needed, (later on we see that this latter correction to 
the subleading asymptotics is related to a polymerization mechanism), and if we set 



-ION 



^" ^2.066' 



iV^/^ (319) 



then we get an extremely good agreement with MonteCarlo simulations. It is interesting to 
note that the exponential scaling factor {j^)"^^^ ^ very close approximation to 

[cos-i(-)U=4]-'°^ (320) 
n 

suggesting that the ratio p'j^^'"/px is, as expected, of geometrical origin. We wish to stress 
that we do not have yet a convincing geometrical proof that the scale factor c„ is actually 
given by (|32CI|) , and leave this as a conjecture to be proven by the interested reader. 



The above choice (|319|) for the rescaling factor c„ provides 
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with /c2 < ^2"''^ — e, and 

~kZ'\k2) = 10 (log 



(1 - 2r]„ 



\{l-2r)m.ax) 



Vmax , 

I /max \ 



l-3r/„,,)(i-3w.) 



log 



2.066 



(322) 



(with k2 > + e). 



(Note that the subleading correction N^/"^ is irrelevant for the characterization of ( p21| ) and 

(HI)). 



The rescaled critical value of so obtained gives the critical /C4 for any simply-connected 
closed 4-manifold, in particular for the 4-sphere S^. Indeed, it provides an extremely good 
fit with the numerical data relative to the characterization of the cxo-volume line in the 
(^4, ^2)- couplings space. Explicitly, if we denote by Kf"^* {Mont. Car) the critical value of /C4 



obtained from Montecarlo simulations for dynamically triangulated four-spheres, then 
the following table demonstrates how well this agreement works out 



k2 


= 0.0. 


.K2''''{Mont.Car) 


= 0.98. 




= 1.07 


k2 


= 0.2. 


.K2"'{Mont.Car) 


= 1.45. 




= 1.47 


k2 


= 0.3. 


.K2"'{Mont.Car) 


= 1.65. 




= 1.67 


k2 


= 0.4. 


{Mont. Car) 


= 1.89. 




= 1.87 


k2 


= 0.5. 


.K2"\Mont.Car) 


= 2.07. 


. kz'' 


= 2.07 


k2 


= 0.6. 


{Mont. Car) 


= 2.30. 


. kz'' 


= 2.33 


k2 


= 0.7. 


.K2"\Mont.Car) 


= 2.54. 


. kz'' 


= 2.55 


k2 


= 0.8. 


.K2"\Mont.Car) 


= 2.78. 


. kz'' 


= 2.78 


k2 


= 0.9. 


.K2"\Mont.Car) 


= 3.02. 


. kz'' 


= 3.01 


k2 


= 1.0. 


.K2"\Mont.Car) 


= 3.25. 


. kz'' 


= 3.25 


k2 


= 1.1. 


.K2"\Mont.Car) 


= 3.49 . 


. kz'' 


= 3.49 


k2 


= 1.3. 


{Mont. Car) 


= 3.97. 


. kz'' 


= 3.97 
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k2 = l.5... K2"\Mont.Car) = 4.47 . . . i^^* = 4.46 
A;2 = 1.8 . . . K2"\Mont.Car) = 5.20 . . . K^"^ = 5.19 

(323) 

(note that the numerical data are affected by an error of the order ±0.02, and the anahtical 
data are rounded to the largest second significant figure). According to the distinct asymp- 
totic regimes provided by theorem and lemma |291| , we have used ( ^22| ), (evaluated for 



Vmin = 1/5), in the range for < A;2 < 0.5 , (however we stress that in such region, and in 
particular near ^2 — 0, the MonteCarlo data are rather unreliable being strongly polluted 
by finite size effects). For 0.6 < k2 < A;2^'** — 1.38 we have used ( |321| ). Incidentally, it is not 
difficult to check that the use of ( P21| ) works quite well also for < A;2 < 0.5, with just a 
sligth discrepancy between the MonteCarlo simulations and the analytical data. Finally, for 
/c2 > kf^^ we have used (|322|) . The numerical data in (323) come from simulations [4^ which 



are already a few years old, however we wish to stress that the agreement with the analytical 
estimates holds beautifully also when using more recent and accurate MonteCarlo data, a 
fact that confirms that our simple counting strategy based on the use of the generating 
function of partitions of integers px has a sound foundation in the geometry of dynamical 
triangulations. 

B. Polymerization 

A well known consequence of the numerical simulations in 4- dimensional simplicial 



quantum gravity is the onset of a branched polymer phase [TT[] in the weak coupling regime, 
i.e., for k2 > k^^*. This polymerization is very effectively described by the asymptotics of 



theorem 17 and lemma P91|, and it is generated again by a mechanism of entropic origin. 



Let us start by noticing that if in the expression yielding for the canonical partition 



function (|245|) we replace the upper integration limit rjmax = 1/4 with the limit rj = 1/3, 
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associated with the algebraic singularity of the integrand at 77 = 1/3 (and corresponding to 
the unphysical average incidence b{n,n — 2) = 3), then we would get a canonical partition 
function W{N, kn-2)poiy whose large N limit is provided by the asymptotics ( p96| ), viz., 



T—n—l 



n-2 



■ exp 



rl f ^.x, A{kn-2) + 2 ^ 

[-n[n + 1) In 



(1 + 0(Ar-3/2)) ^324) 



for all ^2"'^ + 6 < k2 < +00, where k^^^ corresponds to the lower crossing rj* — >■ rjmin = V^) 
and (5 G is a suitably small number. 

In the range k^^-^ + 6 < k2 < kf^^ — e this is exactly the expression for the asymptotics 
of the canonical partition function W {N , kn-2) ef f , but this identification breaks down as 
soon as k2 approaches kf^^. The geometrical reason for such a behavior is very simple: 
as ^2 —>■ ^2*^**, the possible maximum value of the (positive) average curvature for a given 
volume, {vol{M) oc N), is saturated: no more positive average curvature is allowed for that 
volume, (think of the average curvature on a sphere of given volume: it is maximized for the 
round metric corresponding to that volume). On the other hand, ( ^24| ) as it stands, would 
allow, as k2 grows, for more and more (positive) average curvature at fixed volume, and in 
order to comply with the constraint relating volume to positive curvature we have to turn to 
the correct weak coupling description of the partition function, viz., ( pOTp . However, (|324|) 



can still be interpreted as the partition function of a large positive curvature PL manifold 
of given volume, but now the manifold in question is disconnected. Such a disconnected 
manifold is generated by means of a polymerization mechanism, namely it results as the 
(Gromov-HausdorfF) limit of a network of many smaller manifolds carrying large curvature. 
Such curvature blobs are connected to each other by thin tubes of negligible volume so that 
the sum of the volumes of the costituents manifolds adds up to the given total volume, (the 
standard example is afforded by a collection of highly curved small spheres connected by 
tubules of negligible area, vs. the standard round sphere of the some total area). The basic 
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observation is that for /cg*"** + e < ^2 < +00, the effective entropy ( P24| ) is larger than (pOTj) : 



thus in the weak couphng region, for a dynamical triangulation is energetically favourable to 
polymerize rather than to comply to the curvature- volume constraint with a large connected 
manifold. 

The above remarks can be easily formalized. According to (^), the average curvature 
(around the collection of bones {B}) on a dynamically triangulated manifold is proportional 
to 



1 1 JV„_2 

<KiB) >=-E^15M^)«7v E 

Namely 



27r - q{B) cos"^ - 
n. 



(325) 



, , 27r-6(r2,n-2)cos-i(i) , , 

< K{B) >oc — — (326) 



The canonical average, << K{B) >>, of < K{B) > over the ensemble of dynamical 
triangulation considered is, up to an inessential positive constant, provided by (see ( P45| )), 

« KiB) »= rf,^=;i!ii=BtL[2.,- - cos-(i)U-4]e(--^i)/(^) x 

JVrmn ^ ( 1 - 3?/) ( 1 " 2?/) ^ 



X 



^(1 - 3r7)(l - 2r]) 



(327) 



For k2 = k2"^ — e, with e G sufficiently small, the maximum for f{r]), (see (|245| ) ) , is 
attained for an rj* which uniformly approaches 77 = | as e — ^ 0. Correspondingly, the inte- 
grands in ( P27| ) are peaked around 77 ^ 1/4 and << K{B) >> is positive, (explicit expression 
can be easily obtained from the uniform asymptotics of the previous paragraphs). Let us 
now extend ( p27D to rjmax = I, namely consider the ensemble average, << K{B) »poiy 
of < K{B) > over W{N, kn~~2)poiy It is trivially checked that << K{B) »poiy formally 
remains positive for all k2"^ < k2 < +00, (corresponding to the maximum 77* ranging from 
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i < Tj* < |). But according to lemma |l|, there cannot be connected four-manifolds (re- 
gardless of topology) with 6 < 4, (i.e., > |), because for any dynamically triangulated 
4-manifold 

hNi > 2N2 - 5x (328) 

An inequality which implies that as 77 — |, the average curvature < K{B) > attains its 
maximum oc [27r — 4 cos^"'^(|)]/4 compatibly with the given volume, (oc N). 

There is a natural way of circumventing the implications of the constraint ( |32^ ): we have 
to remove the connectivity requirement for M. 

Let us assume that we have a collection {M^} of dynamically triangulated manifolds 
of fixed volume oc N^lMa) = N, and where each Ma is generated by a large number m 
of dynamically triangulated 4-manifolds {S(z)}j=i m(«) C Mq,, connected to each other 
by tubes of negligible volume. Such a sequence of manifolds naturally converges, in the 
Gromov-Hausdorff topology, to the disjoint union of a set of m 4-manifolds, viz., 

m 

{MJ^M = I]H(0 (329) 



such that A^4(M) = A^4(Mq,) = A^. The constraint ( p28|) , which holds for each connected 
manifold Mq, of the sequence, carries over to the disconnected limit space M and implies 
that 

m m m 

5Y.N,{E{^))>2Y^N,{E{^))-5Y: xi^) (330) 



namely. 



Thus, if each is topologically a 4-sphere, S"^, and A'2(2(i)) is 0(1), we can make Bm take 
any value in the range 3 < 6 < 4. In particular, if 
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-1 m 

^2(S(.)) = -EiV2(H(z)) = 20 (332) 
m 

and 

-| m 

Mm) = -T.N,{m) = Q (333) 

i 

then for such M = ]JS(i), we get 

6m ^ 3 (334) 

Such an M is the disjoint union of m ~ ^"^g*^^ small 4-dimensional spheres, (recall that the 
standard triangulation of a 4-sphere obtained as the boundary of the standard 5-simplex 
is characterized by the /-vector (A^'o, A^i, N2, N3, N4) = (6, 15, 20, 15, 6)). 



It is trivially checked that for k!^"^ + e < k2 < +00, the effective entropy ( |324D is larger 



than (|307|) . Thus in the weak coupling region, it is energetically favourable to generate 
dynamical triangulations associated with a disconnected manifold M = IJS(i), rather than 
to comply with the curvature-volume constraint on a large connected manifold. From the 
asymptotics ( |3071 ) and ( p96| ) it follows that the critical exponent associated with such a 



polymer phase, ■jpoiy, is related to 7s (^^2) = 2 — n, (see (|315|) ), by the relation 



Ipoly = ls{k2) + ^ (335) 



This jump of | has been observed in the numerical simulation |TT[| exactly around fcg*^** ~ 1.3. 
In such simulations there is also some indication that 7s(A;2 — /c"'**) ^ 0, so that jpoiy is 
exactly the critical exponent of branched polymers, but this evidence is not yet conclusive 
due to critical slowing down near the transition point. 

According to our counting strategy, we still have the freedom to add to the critical 
exponents the term r(n) associated with the count of distinct dynamical triangulations with 
the same curvature assignements. The correct critical exponents to be confronted with the 
MonteCarlo data are then provided by 
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Isih) = Isih) + Tin) 



(336) 



and 



'fpoly 'fpoly ~l~ T{ 



n] 



(337) 



Since the critical exponent of a branched polymer is exactly ""ypoiy = |, and as we have 
shown ( |324| ) seems to geometrically describe a polymer phase, we can use this result to fix 
the value of rin). Indeed, if at least provisionally, we accept this heuristic argument, from 



( pl5|) we find that the critical exponent scaling factor T{n) must satisfy the relation 



2-n + r(n)|„=4 = ^ (338) 



namely, r(?7,)|„=4 = |, a result already anticipated in the expression ( |319| ) of the normalizing 



factor c„. Not surprisingly, this position provides also the correct critical exponent near 
criticality, viz., 

7sih) = (339) 



These latter results, although not fully rigorous, seem to be conforted by the numerical 
simulations. 

We wish to stress that the polimerization phase is not properly speaking described by the 
dominance of a disconnected phase [i.e., by a gas of weakly interacting small 4-spheres), but 
rather by the dominance, in the statistical sum, of the dynamically triangulated 4-manifolds 
Ma generated by a large number of curvature blobs {H(i)}j=i^...^m(«), connected to each 
other by tubes of negligible volume. 

Few remarks are also in order for what concerns the transition between the various phases. 
The transition threshold towards a polymer phase occurs, in our modelling, somewhere 
around the value of k2 such that r/* = rjmax- Explicitly, for a ^2 — kf^^ such that 
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This condition provides ^2"* — 1.387 which appears quite in a good agreement with the 



value indicated by the numerical simulations. For instance in the critical point is 



located around kf"^^ {Mont. Car I) ~ 1.336(±0.006). More recent simulations tend to 
move this point towards kf^^ ~ 1.22. The explanation of this discrepancy is related to the 
observation that rj* = rj^ax is a theoretical lower bound, and that kf^^ changes very rapidly 
with small variations of the ratio Thus for instance, if we set r]max = 1/4.1 we would 

get ^2"* — 1-2. As is clearly indicated by our analysis, this is simply a question related to 
the fact that it is very delicate to exactly locate numerically the true critical point. 



C. Summing over simply-connected 4-dimensional manifolds 



It is important to stress that the asymptotic estimates of theorem |T7| refers to all simply- 
connected 4-manifolds, not just to the 4-sphere case, the only for which numerical data 
are currently available. 

Besides S^, well-known example of closed simply-connected topological 4-manifolds are 
provided by CP^ ^ §2^ the Kummer surface K3 = {[Zq, Zi, Z2, Z3] E CP'^\Z^ + Zf + 
Z2 + Z^ = 0}, and more in general by all nonsingular algebraic surfaces of degree n in CP^. 
The basic topological invariant for closed simply-connected 4- manifolds is the intersection 
form [|l^]. This is a natural symmetric bilinear form, u!m{,), on the second homology group 
H2{M, Z) of the manifold: for a, /3 G H2{M, Z), let a* be the Poincare dual of a G H2{M, Z), 
then c^Af (a, P) =< a*, P >k, where <, >k is the Kronecker pairing between the cohomology 
group and the homology group. The case of 2-dimensional (connected) surfaces S provides a 
direct illustration of what is meant by this construction. Here the relevant homology group 
is Hi{J], Z2), every homology class in Z2) can be represented by a simple closed curve 

C C S, and if Ci, C2 are any two such curves, their intersection number CrC2 = C2-Ci G Z2 
is defined. Note that the self-intersection number C ■ C is zero as soon as the curve Ci has 
a neighborhood which is orientable, (by deforming the curve C — C, while remaining in 
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the same homology class, we can make C disjoint from C), while if such a neighborhood is 
non-orientable, the self-intersection number will be odd, (think of a deformation of the curve 
on a Mobius band). Actually, by Poincare duality we can make Z2) an inner product 

space over Z2 if we use the intersection number Ci-C2 as inner product. It can be proven that 
any two closed connected surfaces Si and S2 are homeomorphic if and only if the associated 
inner product spaces iJi(Si,Z2) and ifi(S2,Z2) are isomorphic (see Milnor-Husemoller's 
book [1^), and that every inner product space over Z2 is isomorphic to /7i(S,Z2) for some 
closed connected surface S. 

In a similar fashion, the intersection form ujm{c(,I3) on H2{M,'Z) may be thought of as 
representing the transversal intersection of two oriented surfaces a and f3 in the compact, 
simply-connected 4-manifold M, a and f3 being cycles representing elements in H2{M,Z). 
Again, we can think of H2{M, Z) as an inner product space over Z if we use the intersection 
number ujm{c(, as inner product. Moreover, it can be proven that any two closed oriented 
simply- connected 4-dimensional manifolds Mi and M2 have the same homotopy type, [i.e., 
there exists an orientation preserving homotopy equivalence Mi —>■ M2), if and only if 
the associated symmetric inner product spaces H2{Mi,Z) and H2{M2,Z) are isomorphic, 



(Whitehead's theorem [|14|). Actually, as in the case of surfaces more can be stated in 
terms of this isomorphism, since the set of simply-connected oriented closed 4-manifolds 
is in correspondence [|1^] with the set of all unimodular, (any matrix representing uj{,) 
has determinant ±1), bilinear symmetric form, (the manifold is unique when the form is 
even, and there are exactly two distinct manifolds when the form is odd), (Freedman's 
theorem [Q). However, if the 4-manifolds in question carry different iable structures, the 
correspondence is with positive-definite intersection forms, (Donaldson's theorem ||14|| ). 

Let us denote by W{M; N, k2)eff the canonical partition function enumerating the dis- 
tinct triangulations on a 4-dimensional simply-connected manifold M of given topology and 
given volume (N). The set of such M is formally labelled by the set of all positive-definite, 
(we are dealing with PL 4-manifods), unimodular bilinear symmetric forms ujm- Accord- 



138 



ing to theorem the asymptotics (|252|) refers to all simply-connected 4-dimensional PL- 



manifolds of fixed volume, thus we can write 



eff 



= iy(§^ iV, k2)eff + W{f X iV, k2)eff + W{K^- N, k2)eff + ■■■ (341) 

For large but fixed, the dynamical triangulations entering ( |341| ) are of bounded geometry, 
thus according to theorem ^, there are only a finite number of distinct homotopy types of 4- 
manifolds contributing to the above sum, and Whitehead's theorem mentioned above implies 
that likewise finite is the number of distinct types of (isomorphism classes of) intersection 
forms u entering ( p41|) . 



Thus from ( |341| ) we get 
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W{S'; N, k2)eff + W^(§' X iV, k2)eff + W{K^- N, k2)eff + ■ ■ ■ (342) 



where 75(^2) is the scaled critical exponent ( |337D . The identity ( p42| ) makes clear that 
the expression ( p291| ) for the canonical partition function being common to all dynamically 
triangulated 4-manifold (closed) with trivial fundamental group, tends to overestimate the 
canonical partition function of any 4-manifold M in (|342|) . However, guided by the anal- 



ogy with the 2-dimensional case, and by the agreement between our analytical results and 
the MonteCarlo simulations for the 4-sphere, we may tentatively assume that all simply- 
connected DT manifolds contribute in (|342| ) to the total entropy W{N, fc2)e// with the some 
exponential behavior and that the topological differences come about only in the subleading 
asymptotics. Since ( P41|) is a finite sum, ( |342| ) implies that the subleading asymptotics in 
each entropy function W{M; N, k2)eff is of polynomial type, (at least in the range of k2 
considered). Thus, we are led to the following ansatz: 

W{M; N, k2)eff = 10c„ ■ ^e^^"'('=2)^Ar^^("'^2)-3(i ^ _ _ ^343) 

V 27r 
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where 73(0;, ^2) is the critical exponent for the topological 4-manifold characterized by the 
intersection form uj. Since by hypothesis the exponential leading term does not depend on 



topology, the above ansatz allows us to rewrite ( p42| ) as 

^ jY73(aj,A:2)-3 _ jY74fc2)-3 (344) 

The above hypothesis is fully consistent with the finiteness result for the topology of 
dynamical triangulations of bounded geometry discussed in Section V, (theorem accord- 
ing to which the number of distinct homotopy types for such triangulations grows at most 
polynomially with N . 

In order to be more specific, let us recall that the rank r{uj) of a unimodular sym- 
metric bilinear form is the dimension of the space on which uo is defined, namely r{uo) = 
dim[H2{M, Z)]. The signature s{uj) of u is defined by choosing an orthogonal basis bi, . . . ,bn 



for the (rational) inner product space Q (g) H2{M, Z), then JT^ s{iu) = s{Q (g) H2{M, Z)) G Z 
is the difference between the number of basis elements hi with uj{bi,bi) > and the number 
of basis elements bj with uj{bj,bj) < 0. Roughly speaking, s{ijj) is the number of positive 
eigenvalues minus the number of negative eigenvalues of uj (thought of as a real form). Fi- 
nally, the type of u is defined by stating that u is even (or of type II) if u{a, a) G 2Z for 
all a G Il2{M,Z), [i.e., the diagonal entries of u are even), otherwise u is odd, (or of type 
/ [0). The relevance of such definitions lies in the fact that any two (indefinite) inner 
product spaces {Il2{Mi,Z),uJi) and {Il2{M2,Z),uj2) are isomorphic if Ui and UJ2 have the 
same rank, type, and signature [|lj]. According to this observation, we can replace the sum 
over the possible intersection forms in ( p44|) with a sum over the possible ranks, signatures 
and types. To this end, let us introduce the functions A± (r, s) enumerating the distinct 
(isomorphism classes of) intersection forms uj of rank r{uj) = r, signature s(co') = s of even 
type, (v4+(r, s)), and of odd type, (v4_(r, s)) on a 4-manifold of bounded geometry. Then, 
we can write (13441) as 
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f s 

^ ^[A+(r, s) + A_{r, s)]Ar^^('^'^2) _ ^7.(^2) ^345) 

where f < +00 and s < +00 are the largest value for the rank and signature allowed on a 
4-manifold of bounded geometry. In analogy with the 2-dimensional case, (where the critical 
exponent depend linearly on topology), let us make the further natural ansatz 

%{uj,k2)=ar + ps + S (346) 



where a, (3, 6 are constants possibly depending on the inverse gravitational coupling ^2- 
Notice that since H2{S^, Z) = 0, 6 is the critical exponent for the 4-sphere. 

The ansatz concerning the structure of the effective entropy of each simply- connected 
4-manifold M, implies that 

^ ^[A+(r, s) + A_{r, s)]iV"'^+/3^+5 = AT^^te) (347) 

which is a very strong bound on the set of possible intersection form for 4-manifolds of 
bounded geometry. Indeed, without the control on geometry, the number of inequivalent 
intersection forms, say of given rank r, grows with r at an extremely fast rate. To make 
a rather well-known example, (see Milnor-Husemoller book quoted in [|1^]), consider the 
subcase of positive-definite intersection forms of odd type. The number of such forms is 
known to grow with the rank r at least as 



\27^e^/e J 

Thus, the interplay between our analytical approach and MonteCarlo simulations seems to 
indicate results of great potential interests in the borderline between 4Z)-simplicial quantum 
gravity and the topology of 4-manifolds. 
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D. Higher order phase transitions for simply connected 4- manifolds 



To conclude, we will directly prove that for a dynamically triangulated 4- dimensional 
simply-connected manifold M, {e.g., the sphere S'^), the jump in the critical exponents 
7(M, A;2), as the inverse bare gravitational coupling k2 approaches the value k2 = k^'^^, 
corresponds to a higher order phase transition of our statistical system. This result settles 
down in the affermative a long standing debate about the real nature of the critical point k^'^* 
characterized by the MonteCarlo simulations. As we shall see, some of the contraddictory 
results of the numerical simulations find their rationale in our analysis. We wish to stress 
that most of the material in this section is based on joint work with D. Gabrielli and G. 
Gionti, and we wish to thank them for granting us permission to present part of their results 
here. 

The optimal way for discussing the onset of a phase transition for simplicial quantum 
gravity on a 4-manifold M, would require the characterization of the two-point function 
associated with the set of all dynamical triangulations of M. The evaluation of such two- 
point function requires the enumeration of all distinct dynamical triangulations on M with 
two marked simplex whose distance is declared to be some running parameter r. As already 
remarked, this counting problem is, at the moment of writing these notes, not yet solved. 
However, a reliable indication that our system exhibits a phase transition may came also 
from exhamining the simplicial version of the curvature susceptibility. This is defined by 

^ 1 d'^\nW{E.^;N,k2)eff _ ]_ 

where N = ^n{n + 1)A'"|„=4 = 10 A^. Since we work at fixed A" = A4, {i.e., at fixed volume), 
such a xr can be equivalently written as 

XR = {{1?) - (350) 

where as usual b{n,n — 2) = ^n'{n + 1)t^- Thus xr characterizes the fluctuations in the 
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<X^>-<X>' 



(349) 



average incidence b on the set of all possible dynamical triangulations (at fixed volume) on 
For this reason, it is not difficult to show on heuristic grounds that such a xr can be actually 
considered as the dynamical triangulation counterpart of a connected curvature-curvature 
correlation function: 

connected 

(351) 

It is worth noticing that such a xr has been considered for describing the onset of phase 
transitions also in computer-assisted simulations. However, its direct use in such a setting 
is hampered by the fact that in such computer simulations one is forced to allow for volume 
fluctuations. Moreover, it has been observed that the peak of such xr is rather asymmetric 
and difficult to locate. Thus, in MonteCarlo simulations it is difficult to provide a clear cut 
interpretation of the behavior of xr- 

Such difficulties in using the curvature susceptibility have a clear origin in the analytical 
computation of xr- Indeed, contrary to what one would like to think, at the critical point xr 
does not diverge with A^, it rather exhibits a finite jump. According to standard statistical 
mechanics, such a jump may or may not be associated with the onset of long range order, 
and not surprisingly, for characterizing the order of the transition at ^2^**, we are naturally 
led to a finite scaling argument a ' la Fisher not much dissimilar from the kind of argument 
adopted in the actual MonteCarlo simulations. 



Adopting the notation iof theorem |1^ and lemma |T3|, the possible existence of long range 
order at the critical point k2 = k'^^^ rests on the following 

Theorem 18 As \k2—kf^^\ and N varies there exist pseudo critical points k2 = kf^^{N) 
characterized by the condition 

ifjma.{kr\N))\f^ = 5 (352) 



143 



where < S <1, is a fixed ( suitably small) constant (indipendent from N-see ( ^91\ )), such 
that 



(353) 



N 



and 



lim XRih = kr) = ^ 



fvviv* 



(354) 



This indicates that at the critical point ^2 = ^2 — 1-387 the dynamically triangulated 4- 
manifold M undergoes a higher order phase transition possibly associated with the onset of 
long range order. 



Proof. According to the uniform asymptotic estimation of the canonical partition function 



W{N, k2)eff provided by theorem |T7| and its corollaries, the approach to the critical point 



^cr«t ^ 1.337 is characterized by the parameter (see (|249|) ) , 



i^max{k2)=Sgn{T]* - Vmax) \/ ^[f [t]*) - f{r]max)] 



(355) 



through the condition 



i^maxik2)\/ N ^ 



(356) 



Whereas, the separation between the weak (strong) coupling regions and the critical 
region corresponds to values of the inverse gravitational coupling k2 such that 

^|Jmaxik2)^/^ ^ I- (357) 



It follows that if we choose a 6, with < 5 < 1, sufficiently small, then through the 
condition (|352|) , we select the rounding interval feg*"** — k'^'^*{N) < k2 < k'^^^ + k'^^^{N) 
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corresponding to which the critical behavior occurs. As the notation suggests, this condition 
characterizes a value of the coupling /c2 = k^ '^^^N) according to 

ipmaAkT''{N))\/^= 6 (358) 



which localizes the maximum rj* of the function f{7], ^2) to an exponentialy small neigh- 
borhood of the transition point 77 = rj"^"-^. For a given A^, this k'^^*{N) maximizes the 
corresponding value of the curvature susceptibility xr{N), and as such it defines a pseudo 
critical point. Explicitly, from the definition of xr and from the proof of theorem we 
compute 

1 92 



Xr 



Ndk 



2 



In / G(p)e~(^+^)(5P +^ma.p)^p ^359) 



By exploiting the uniform asymptotics estimation of the integral appearing in the right 
member of ( P59| ), we can write, up to order 0{N~"^'^), 

Ndk,' \{Ny/^t,{Ny {N) toiNyi ' 



where z = ipmax{kn-2)\l N, (henceforth we will omit the annoying 1 whenever appears). 
If we introduce explicitly the expression of the Whittaker functions wq{z) and W-i{z) in 
terms of the parabolic cylinder function Dr{z), viz., Wr{z) = r(r — l)e^^^'^Dr-i{z), then a 
long but straightforward computation provides 



+ 2^-- dkl + 



iAln(^E^ + ^E^| (361) 



In the critical region characterized by ipmax{k2)\J N « 1, the asymptotics ( |298D of the 
parabolic cylinder functions shows that, to leading order in A^, we have 
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" 2 ( ^"^""^^ 'V 2^"^- dkl 



In order to compute the above derivatives, let us consider the following limit 



lim Ml^!^= lin, 



By I'Hopital rule one easily obtains 



crit\ 2 



and thus, not unexpectedly, we obtain 



(363) 



. r I = 7r\'^2=kr^ = ^V*a - 3r7*)(l - 2r^*). (364) 
But this latter is exactly the expression of —f^^ijf^-, thus one computes 



lim If-f = (365) 



This latter result implies that for |/c2 — k^ \ small enough we may write 



which proves that at criticality the susceptibility does not diverge with A^, but rather it has 
a finite jump. 

We can rewrite ( |367| ) as 

= 2(^2 -Ir'y ^^^^^ 

Thus, if k^°'^{N^ is defined as the value of k2 for which 

i'maAk2 = k2{N))4^ = b (369) 

we obtain 
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XR{S,N)=^jj (fcrW-^n^ (370) 

which shows that at fixed N, and for a given 6, (fixing the rounding interval), xr attains its 
maximum at k'^'^^{N). This indicates that k'^^^{N) is a pseudo critical point for the given 
rounding interval and ( |366| ) provides the way such pseudo critical points approach, as 
increases, the true critical point /C2"*. Indeed, from ( |366D and ( ^69| ) we obtain 



-f (r)*) 

\kT\N) - k^"'\ ~;v»i !1 (371) 



N 



as stated. 

As is well known, the scale fixing the rounding of a first- order transition takes place on a 
scale proportional to the inverse of the total volume of the system. In our case, this rounding 
takes place on a scale proportional to the square root of the volume of the manifold. Thus 
the phase transition at kf^^ is not a first order transition and it is possibly associated with 
the onset of long range order in our statistical system. □ 

This result strongly suggests that the ensemble of 4-manifolds will develop a genuine 
extension when k2 k2^"^ ^ and make this transition point a potential candidate for a non- 
perturbative theory of quantum gravity. As already emphasized, in order to carry out this 
program (and eventually remove all the possihly's in the above statements concerning the 
existence of long range order), we have to compute the two point function and characterize its 
scaling properties at criticality. The almost spectacular agreement between the analytical 
approach and the MonteCarlo data indicates that the picture of 4-dimensional simplicial 
quantum gravity which would emerge from this program should be very similar to that one 
already elaborated on the basis of computer simulations. However, a delicate point emerges 
here for what concerns the assessment of the nature of the transition at k'^^*. Some recent 
MonteCarlo simulations have cast doubts to the standard wisdom concerning its second 



order nature. It has been observed a double peak in the distribution of order of vertices, 
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and this has been interpreted of some indication of a first order nature of the transition. On 
the other hand, finite size scahng results, as above, point to a higher order transition. Prom 
the numerical simulations point of view the situation is not yet settled. In fact the situation 
might be quite similar to the one we encounter for compact QED in four dimensions: we 
have a strong coupling phase were the string tension is different from zero (to be compared 
with the crumpled phase in 4-dimensional gravity). The string tension scales to zero at 
the critical point (in 4D-gravity this would correspond to the mass parameter scaling to 
zero at the critical point). At this critical point it is unclear if there is a two- phase signal, 
and if the transition is first or higher order. A two-phase signal has first been observed, 
then it has disappeared in a sort of see-saw game, and in a way akin to the situation we 
are experiencing in 4D-gravity. This similarity persists in the weak coupling phase: for 
small coupling constant we have the coulomb phase in QED, viz., for all coupling constants 
smaller than the critical one we have a continuum massless phase, and indeed in 4D-gravity 
for all (bare) gravitational coupling constants smaller than the critical one we have branched 
polymers, which indeed admit a continuum limit. Since, as we have shown, there is more 
than a superficial agreement between MonteCarlo simulations and our analytical approach, 
one wonders if there is a way of understanding the origin of a possible two-phase signal and 
how this may be compatible with the results above concerning the higher order nature of 
the transition. 

The explanation is that the uniform asymptotics controlling the leading behavior of the 
canonical partition function may generate a bistable behaviour: we have discussed the ap- 
proach to criticality as if the only leading contribution to W{k2, N)eff comes from developing 
in a neighborhood of ipmax{k2)y/N — > 0. However, as it is easily checked, f{rjmin) is not much 
different from f{rjmax), and W{k2, N)cff may receive an important contribution from devel- 
oping around ipmin{k2), too. It is this further contribution that may generate a double peak 
and similar hysteresis effects, in a way completely analogous to what happens in the three- 
dimensional case and which we discuss in details in the next section. These effects tend to 
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obliterate the nature of the transition. They are more manifest and more important in the 
three-dimensional case essentially because the difference \f{rimax) — f{Vmin) \ is smaller than 
in the four- dimensional case. 



E. Hysteresis effects in the three-dimensional case 

The confrontation between MonteCarlo simulations and the analytical estimates is quite 
satisfactory also in the three-dimensional case. 



Numerical data in dimension n = 3 often employ |43| the pair of variables A^s and A^o 
rather than N3 and A^^i, thus we have to adapt our asymptotic estimates to this new pair of 
variables before comparing with the simulations. This can be trivially done by noticing that 
the Dehn-Sommerville relations Nq- Ni+ N2- N3 = and N2 = 2N3 imply A^i = A^g + A^o, 
so that the partition function can be written as 

Z[h{ko), fcs] = E f E W{X = N + No, A^)e^^(^=^+^«) I e'^^^^ 

N \X=N+No J 

= E( E W{X = N + No,N)e'''''Ae^^'-'''^'' (372) 

N \\=N+No J 

where N=N3, and A = A^i. 

The delicate point to stress is that in dimension = 3, we cannot any longer use the 



results of theorem and lemma |T3| as they stand. The basic observation is that the function 
f{r]), (now expressed in terms of the variable A^o), is, in the interval 

concave function symmetric with respect to the point 77* where it attains its maximum, viz., 

fiVmin) — fiVmax) (373) 



Since the parameters driving the uniform asymptotics of the saddle-point evaluation of the 
canonical partition function are, (see ( p48| ) and (|249|) ) , 



^mm(fcn-2) = -Sgn{r]* - r]min)\l2[f{7]*) - f\rjrnin)] 



(374) 
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and 



^max{kn-2) = Sgn{r]* - Vmax) \/ 2[f {t]*) - /(?7„ax)] (375) 

then ( P73| ) implies that the correct asymptotics of the canonical partition function, in the 
range rjmm < < Vmax is provided by a linear combination of the two uniform asymptotics 
obtained by developing the integral (^55|) in terms of both ipminiki) and ipmax{ki), viz., 



Theorem 19 Let us consider the set of all simply- connected 3 -dimensional dynamically 
triangulated manifolds. Let , and /cg^** respectively denote the unique solutions of the 
equations 



3^ Aiko)' ' 6 



1 _ _ 1 

77 Vmin 



6 



(376) 



3*-"^ A{ko)' ' 6 



1_ _2 
y 



(377) 



Let < e < 1 small enough, then for all values of the inverse gravitational coupling such 
that 



4"^ + e<ko< k^"' - e, 



(378) 



the large N -behavior of the canonical partition function W{N, ko)eff for 3-dimensional sim- 
plicial quantum gravity on a simply connected manifold is given by the uniform asymptotics 



^7V-7/2g(JV+l)/(r,„,„) 



--WQ{'ll)min{kQ)\l N) + -^W-iitpminiko)'^ N) 



N 



+ 



+ _r!!_]V-7/2g(^+l)/(r?ma.) 



Wo{'l/jrnaxiko)\f^) + (V^max (/^o) V A^) 



N 



(379) 



where the constants ao, ai, and ^0, C,i are given by 
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tto 



Vmin (1 '^Vn 



\J{1~ 3r7mi„)(l - 2rijnin) fviVmin) 



(380) 



^^-1/2 (1-27^*)- 



^mm(/^o) ^(1 - 3r7*)(l - 27]*) 



'ipmin{ko)\/-fm{V*] 



(381) 



\J{1- 3?7max)(l - 2?7mai') MVmax) 



(382) 



^0 r7*-V2(i_2r^*)- 



4^max{ko) ^{l - 3^*)(l _ 2r/*) 



i^max{,kii)^-fr,r,{ri*) 



(383) 



From this result it follows that k^"^{kQ) can be formally obtained from an expression 
similar to (|309|) by the substitution ki ako and kf"^^{ki) — > k^"^{ki — > afco) — a^o, where 
< a < 1 is a suitable constant, viz., 

In c„ 



where 



^3 [ko) = 6 In ^ ako + 6 Injv^oo 



y4(aA;o 



27 / 27 

ye°'^" + l + W(ye"^o + i)2_i 



1/3 



+ 



(384) 



+ 



27 

_J_ ako _|_ 1 

2 ^ 



27 

r_e"fco + 1)2-1 
^ 2 ^ 



1/3 



(385) 



If we fix the rescaling factor c„ according to 

,1 



Cos-l(-)]-('5/5)iV 

3 



(386) 



and set the constant a ~ 0.16, then it is not difficult to check that the rescaled kf^^{ko) 
has a distribution of values, as ko increases, which is in very good agreement with the 
corresponding values determined by MonteCarlo simulations. Explicitly we get 
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fco = O...K^'''\Mont.Car) 
ko = l...Kl"\Mont.Car) 
ko = 2...Kl"\Mont.Car) 
ko = ?,...Kl"\Mont.Car) 
ka = A...Kl"\Mont.Car) 



= 2.050 . . . 7^3^"* = 2.04 
= 2.055 . . . isTg"'** = 2.07 
= 2.070 ... iTg""* = 2.11 
= 2.100...ir3""* = 2.15 
= 2.220 . . . Kl"^ = 2.20 



The critical point (corresponding to the crossing of the upper hmit rjmax = |) occurs for 
f^cnt ^ 3.845 which agrees well with the location of the critical point kQ^^ {Mont. Car) ~ 3.9 
obtained by the MonteCarlo simulations. 

Some aspects of the nature of the transition at the above critical point k^^^^ are also 
rather easily characterized. It is standard wisdom among groups performing MonteCarlo 
simulations that this transition should be a first order transition, a credence fostered, among 
other things, also by the fact that in three-dimensional gravity we expect no dynamical 
graviton to be present, (but this is a very misleading concept in a non-perturbative approach 
to quantum gravity!). More seriously, the assumed first order nature of the transition is 
associated to quite pronounced hysteresis effects, the same hysteresis effects that (even if 
at a much less pronounced level) seem to be manifest also in some recent simulation of 
4-dimensional simplicial quantum gravity. 



From theorem 19 is clear that as ko — > ko^^^ the system shows a cross-over behavior 



with the phase connected with small values of ko, since as ko approaches /cg"* both ipmax 
and ipmin tend to go to and the system shows a bistable behavior yielding for an hysteresis 
associated with the distinct phases described by ipmax and ipmin- This bistability, if strong 
enough, can be consistent with a first order nature of the transition at k^"*. 

It is very important to develop the above analysis in order to study in detail the neighborhood 
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of the transition point either in 3-dimensions as well as in 4-dimensions. The analytical ap- 
proach developed here already provides useful indications, but it can be certainly improved, 
hopefully to bring to full power dynamical triangulations as a most effective approach to 
deal with non- perturbative issues in euclidean quantum gravity. 
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X. APPENDIX: THE TANGENT SPACE TO THE REPRESENTATION 

VARIETY 



Given the representation 6:7ri(M) — > G let us consider the associated representation 
on the Lie algebra fi generated by composing 6 with the adjoint action of G on g, viz., 
9:7ri{M) — >e G -^Ad End{Q), (henceforth we will always refer to this representation ). The 
representation 6 defines a fiat bundle, flg, over the underlying space M of the dynamical 
triangulation considered. This bundle is costructed by exploiting the adjoint representation 
of G on its Lie algebra g, i.e., Ad: G — > End{Q), and by considering the action of 7ri(M) on 
Q generated by composing the adjoint action and the representation ©: 

fle = M X fl/TTi ® [Ad{Q{-))]-^ (388) 

where tti ® [Ad{Q{-))]~^ acts, through 7ri(M), by deck transformations on M and by 
[yl(i(0(-))]~^ on the Lie algebra q. More explicitly, if (3"i, a2 are simplices in M, and gi, 
g2 are elements of g, then &i,gi) ~ (a"2, g'2) if and only if (T2 = aia, and g2 — [Ad{Q{a))]~^gi, 
for some a e 7ri(M). 

Given a dynamical triangulation and a fiat bundle qq over M = \Ta\, we can define 
homology and cohomology groups with value in q^, (namely homology and cohomology with 
values in the sheave of Killing vectors of Mriem)- On Ta we define a p- chain with values 
in to be a formal sum Qi'^^ where erf e Ta is an ordered p-simplex and gi is an 
element of the fiber of qq over the first vertex of erf. Namely, first we consider chains with 
coefficients in the Lie algebra g, viz., gj&j with gj e g, and then quotient the resulting 
chain complex C(M) C?) s by the action of tti ® [Ad{9{-))]~^. We denote the group of such 
chains by Cp{M,Qg). There is an action of tti{M) on the above chains expressed by 

ai^gj^]) - EiiMOiaWgM^j) (389) 
j j 

for any a G 7ri(M), {i.e., we are considering g as a 7ri(M)-module). This action commutes 
with the boundary operator, and we may define homology groups H^{M,Qg) with local 
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coefficients in the flat bundle 29- By dualizing one defines the cohomology H*{M, Qg), which 
enjoys the usual properties of a cohomology theory. 



These homology groups are naturally related to the geometry of the infinitesimal de- 
formations of representations of the fundamental group of M. Given the representation 
6: 7ri(M) G let us consider the associated representation on the Lie algebra g generated 
by composing 6 with the adjoint action of G on g, viz., 6: 7ri(M) ^ G ^ End^g), (henceforth 
we will always refer to this representation). Then a 1-cocycle on 7ii{M) with coefficients in 
6' is a map c: 7Ti{M) q such that for 7, 5 G vri(M) we have 

c{jS) = c(7) + 0(7) ■ c{S) (390) 

The space of such cocycles (often called crossed homomorphisms) is denoted by Z^{tti{M), g). 
We denote by B^{ni{M),g) C Z^{'Ki{M), q), the subspace generated by the cocycles which 
are 1-coboundaries, namely the set of c G Z^(7ri(M), g) for which there exists a (7 G g such 
that 

c(7) = e{'y)g - g (391) 
for all 7 G 7ri(M). The first cohomology group of vri(M) with values in g is then defined by 

The definition of the higher cohomology groups is somewhat more complicated, and the 
interested reader may consult In order to relate H^{7Ti{M), g) to the infinitesimal 

deformations of the representations of 7ri(M), let -^°'"(""^(^^)''^) denote the set of all conjugacy 
classes of representations of the fundamental group 7ri(M) into the Lie group G. Notice 
that if 6 and F6F~^ are two conjugate representations of vri(M) in G, then through the 
map Ad{F):g ^ g we get a natural isomorphism between the groups ifj(7ri(M), g) and 
/7,(7ri(M),FflF-i). 

We are interested in understanding the structure of -^°™-(""^(^^)''^) -when deforming a par- 
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ticular representation ^:7ri(M) G through a differentiable one-parameter family of rep- 
resentations 6t with 9q = 9 which are not tangent to the G-orbit of G i/om(7ri(M), G). 
To this end, let us rewrite, for t near 0, the given one-parameter family of representations 



Ot as pg, PI 



dt = exp[tu{a) + 0{t'^)]9{a) (393) 

where a G tti{M), and where M:7ri(M) q. In particular, given a and b in 7ri(M), if we 
differentiate the homomorphism condition 9t{ab) = 9t{a)9t{b), we get that u actually is a 
one-cocycle of vri(M) with coefficients in the 7ri(M)-module gg, viz., 

u{ab) = u{a) + [Ad{9{a))]u{b) (394) 

Moreover, any u verifying the above cocycle condition leads to a map 9t: tti{M) — > G which, 
to first order in t, satisfies the homomorphism condition. This remark implies that the 
(Zariski) tangent space to Hom{iTi{M),G) at 9, can be identified with Zi(7ri(M),0). For 
this reason, the space Z^{tti{M), g) is called the space of infinitesimal deformation of the 
representation 9. 

In a similar way, it can be shown that the tangent space to the Arf-orbit through 9 
is (tti^M) , g) . Thus, the (Zariski) tangent space to ^^"^i'^^^^)''^) corresponding to the 
conjugacy class of representations [9] is //^(7ri(M), g) which, therefore, can be identified 
with the formal tangent space to the representation space. 

It is not difficult to establish an explicit isomorphism between the tangent space 
H^{tti{M), g) and the first cohomology group with coefficient in the fiat bundle gg, H^{M, gg), 
associated with the underlying space M of the dynamical triangulation T^. We choose a 
base Xq G M, (say a vertex of the base simplex ctq), and a corresponding basepoint xq G M, 
the universal cover of M, such that vr(a;o) = Xq, tc: M —>■ M being the covering map. If Qq is 
the fiber of gg over xq, then go is also the fiber of 7r*gg over a;o. Since parallel transport (for 
the fiat connection) is independent from the path on M, we can parallel transport the fibers 
of TT*gQ to go, and this generates a map e: C°(M, 7r*gg) — > go, where C°(M, 7r*gg) is the space 
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of functions assigning to each vertex of M an element of the fiber of 7r*g0 over the vertex 
considered. By exploiting this map, we can define, for every 1-cocycle a on M with values 
in Qg, a function 4>a- 7ri(M, xq) go by setting 

0,(7) = e(7r*a(7)) (395) 

where 7 e 7ri(M, xq), and 7 is a simplicial lift of 7 to M (starting at Xq). Such 0^ is a 
1-cocycle on 7ri(M, xq) with values in q, and the mapping a — > 0^ establishes the seeked 
isomorphism from the cohomology group H^{M, qq) to the cohomology group H^{7ii{M), g). 
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